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Part I: Introduction and Discussion 


‘The reception of radio, radar, and television signals is always accomplished in the 
presence of noise, which exclusive of outside sources—either man-made or natural— 
arises because of the random motion of electrons in the various circuit elements of the 
receiver, principally in the stages preceding the detecting devices. The noise is ac- 
cordingly assumed to be of the fluctuation type; impulse noise is not considered in this 
paper. Knowledge of the spectral distribution and the power output following demodu- 
lation of a frequency-modulated carrier in such noise is essential to any critical theory 
of receiver behavior. Accordingly, the principal purpose of the present work is to in- 
vestigate the spectral aspects of the problem. One of the chief reasons for our interest 
in the spectrum lies in the fact that the signal-to-noise ratio at the output as a function 
of the ratio at the input depends markedly on the spectral distribution, in the so-called 
case of broad-band F M (where the maximum change in carrier frequency may be several 
times the highest modulating components). The dependence in the instance of narrow- 
band FM (for which the maximum deviation in the carrier frequency is comparable with 
or less than the highest significant audio modulating frequency) is not nearly so critical. 

Previous investigations, of which the work of Carson and Fry,' Crosby,’ Blachman,”* 
and Rice* are particularly to be noted here,** have dealt chiefly with the strong-carrier 
case, in which the noise is largely suppressed and appears linearly with the signal in 

*Received May 27, 1948. The research reported in this document was made possible through support 
extended Cruft Laboratory, Harvard University, jointly by the Navy Department (Office of Naval 
Research) and the Signal Corps, U. 8. Army, under contract N5ORI-76, T.O.1. The author wishes to 
thank Miss Marilyn Lang of the Electronics Research Laboratory, who performed most of the calcula- 
trons In this paper 
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the output. Rice’s work* is an exception, in that he determines the spectrum for the 
case of arbitrary input carrier-to-noise power ratio (p), but only for the idealized ex- 
ample of extreme limiting (as do likewise Carson and Fry’ and Crosby’). Further, the 
carrier is not modulated. Blachman’s interesting treatment,’ on the other hand, is more 
comprehensive, since he includes modulation and a discussion of the case of no limiting. 
Similarly, the present paper is more general than Rice’s work in that modulation and 
no limiting are also considered, and is an extension of Blachman’s efforts” since spectra 
and their associated correlation functions) are obtained for all values of the ratio p. 
We consider the case of no limiting because it gives us one extreme of receiver per- 
formance as a function of limiting; superlimiting gives us the other. Since the passage 
from the former to the latter is monotonic, we may be able roughly to estimate the 
behavior for intermediate degrees of limiting. The ability to handle all values of input 
carrier-to-noise power ratios is particularly useful when we desire the spectrum for 
threshold signals (p — 1 or less), an important region when contrasting the maximum 
sensitivities of AM/ and FM reception. Specific calculations of spectra are made, which 
utilize a Gaussian input noise distribution. This spectral shape is usually a somewhat 
better fit to actual responses than the rectangular spectrum. The results of the analysis 
re given in detail in Part II, while Part III contains a treatment of important special 
eases: (1) noise alone, (2) unmodulated carrier and noise, and (3) a sinusoidally modu- 


lated carrier in noise 

The actual nonlinear or demodulating elements in an F'Al receiver are too involved 
to vield satisfactorily to direct treatment, but a possible model of these essential circuit 
r—may be constructed if we assume 


elements the limiter, followed by the diseriminat 
that the physical discriminator is replaced by an “ideal”? one which responds every- 
where linearly with frequency (the ‘‘quasi-stationary”’ hypothesis of Carson and Fry 

Thus, the output current (or voltage) is directly proportional to the instantaneous 


difference frequency between the vave and the central or resonant frequency of the 
[F, limiter, and discriminator bands (all tuned to the same frequency and assumed to 


be svmmetrical). The filter characteristic of the limiter is taken to be wide enough to 


pass the /F portion of the limited signal and noise without distortion due to frequency 
selection. In practice, of course, the real discriminator saturates for sufficiently large 
excursions about resonance, but such behavior IS rare if only an insignificant portion 
of the wave’s energy is redistributed into higher harmonics due to the “spreading”’ 


This is assumed to be true here even for the extremes of super- 


il 


‘tion of the limiter. 


limiting, in which the incoming disturbance is heavily clipped at top and bottom 
Blachman has shown (ref. 3, ch. IV) that when there is no limiting, the characteristic function 

method used by Rice and others® in the solution of A.W receiver problems may be adapted for all carrie: 

the problem in which the idealized discriminator is replaced by one closer to actual eleetronie 


practice name . two diode rectifiers m7 radians out of phase with one another Unfortunately, the 


strengths t 
method breaks down when there is limiting, since it does not scem practicable to determine, in any 
useful form, the distribution functions W; and W.2 of the clipped, filtered noise a 
of the above type. Kae and Sicgert (J. Appl. Phys. 18, 


1d signal WAVES le aving 


the limiter to enter a second nonlinear device 
383 (1947)) have calculated the first-order distribution W, , but it is the second-order density Wy» that 


lysis of the present pape! has been extended by the author to 





is needed for the spectrum. (The 


inelude arbitrary limiting, with the idealized discriminator and similar assumptions as to the width 


of the limite filter response. See Technical Report No. 62, Cruft Laboratory, Harvard University, 
Nov. 24, 1948 
61). Middleton, The response of biased aturated linear and quadratic recti fie to random noise, 


J. Appl. Phys. 17, 778 (1946). 
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The wave leaving the JF and entering the limiter is narrow-band, that is to say, it 
may be represented analytically by 

V(t) = R(t) cos [wot + O(d)], (1.1) 

where f(=)/2r) is the central frequency of the JF, limiter, and discriminator, RP is 

the envelope and @ a phase angle; both R and @ are slowly varying functions of time 


compared with wol. If we let f(7z) be the Fourier transform of the limiter’s dynamic 


characteristic g(V), 1.e. 


fiz) = | e*¥g(V) aV, (1.2) 
the output of the limiter may be written’ 
: gr a 
V,(0) : | [(az) exp |2z R cos (wot + 6@)| dz. (1.3) 
2r Je 
The contour C extends along the real axis from — © to + © and is indented downward 


about a possible singularity at the origin. With the help of the expansion 


exp (7a cos d) = > ie, J,(a) cos no, é = I, én, = 2, (1.4) 
we obtain from (1.3) 
V (2) 2 B(R) cos n(wpt + 8), (1.5) 
v} ere 
ee es a es : ine 
BR) = =" | f(i2)J,(Re) dz. (1.6) 
2r Je’ 
The quantities B,(R) are the envelopes of the n (=0, 1, 2, ---) spectral zones that are 
produced in the limiter by its nonlinear action;"'” 7@ is the phase relative to nwot, the 
th harmonic of the original central (angular) frequency. Only the band centered 
1) is passed to enter the discriminator. The input wave is therefore 
V; B,(R) cos (wot + 4). (49) 
ector language, this input is represented by V, = Ve'*, where real quantities 
orrespond to the 2-components, and imaginary quantities to the y-components. Here 
\ B,(R) and ® = wt + 6. Then we have 
} Vicos ® + zsin &) + Vo(—sin ® + 7 cos &) rV + 0V®, 
since the polar unit vectors (r, 6) are related to the rectangular ones by the transforma- 
tol 
r cos ® + 7sin ®, ft) sin & + 7 cos ®; 


V is the time rate-of-change of the modulus, and V@® is the time-derivative of the 


S. O. Rice, Mathematical analysis of random noise, Bell Syst. Tech. J. 24, 46 (1945). 
8—D. Middleton, Some general results in the theory of noise through nonlinear devices, Quart. Appl. 
Math. 5, 445 (1948). (D. M. I. 








132 DAVID MIDDLETON [Vol. VII, No. 2 


phase, or, what is the same thing, the instantaneous frequency.’ The output /% of our 
idealized discriminator is precisely this second term, namely 


Ev(t K Vb = «B,(R)(we + @), 


where « is a constant of proportionality with the dimensions (sees). In practice the 
high frequency term, represented above as the coefficient of w) , is not passed by the 


discriminator filter, so that the final discriminator output is 
; ; iOK F ; 
E(t) = «xB (R)6 | ((iz) J (Fz) dz. (1.9) 
T JC ; 
Notice that the effect of the limiter is contained in the term B,(/2), the exact form of 


which depends on our choice of g(V), ef. (1.2). We assume that the limiter amplifies 
linearly (for both positive and negative amplitudes) up to a certain level, denoted by 
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Fig. 1. A typical limiter characteristic. 
R, , beyond which point saturation occurs abruptly. The dynamic characteristic is 
shown schematically in Fig. 1. From (1.2) we find that here 
f(iz) = 2B(1 — exp [—7R,z])/(zz)’, (1.10) 


where 6 is the dynamic transfer constant of the limiter and the factor 2 is present be- 


°B. van der Pol, The fundamental principles of frequency modulation, Jour. I. k. E. 93, Part ITI, 


53 (1946) 
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cause both the positive and negative portions contribute to the envelope. Accordingly, 
B,(R) in (1.9) is” 


B(R) e ***") J, (Re) — = BR, O< R<R,, (1.1 1a) 


rT JC rs 


202 i dz 
(] 


* aR, F,(—1/2, 1/2; 3/2: R;/R’), 
(1.11b) 


> 

S 
lA 
> 


JR, | | 
ad ) ( Be) Te Fi 


bo 


The two cases of chief interest are (a), no limiting, and (b), superlimiting. From (1.11a) 
and (1.11b) it is evident that B,(?) = BR (no limiting); and B,(R) = (4/1)BR, (super- 
limiting) in which the rms saturation level for the envelope is much less than the input 
noise level (b,)'’”, i.e. Re <K Qhy . 

\lthough the spectral behavior of the discriminator output depends on a considerable 
number of parameters, chief among them (1) carrier strength, (2) degree of limiting, 
3) power and spectral distribution of the noise entering the demodulating elements, (4) 
deviation from resonance, (5) intensity and wave-form of the modulation, etc., certain 
general observations can be made with regard to these factors and their effect on the 
power output and spectrum. We summarize the more significant features below, referring 
to the appropriate figures and equations, the analytical details of whose derivation are 
available in Parts IT and ITI. 

1. Signal Output (Power and Spectrum). Spectrally speaking, the modulation is 
reproduced without distortion (2.32), since our model of the discriminator responds 
linearly to all frequencies about resonance. In practice, this is only approximate, be- 
cause of unavoidable nonlinearities which result in harmonies of the modulation in the 
low-frequency output. The effect is small, however, if the maximum deviation remains 


mostly the linear portion of the discriminator characteristic. 
Figure 2 illustrates the output signal power P,(0) [Eq. (2.32)] as a function of 


the input carrier-to-noise power ratio p. A significant difference between the limiting 
and non-limiting conditions of receiver operation is immediately evident. For no limiting 
\ 1) and sufficiently large values of p (p* > 1) the output signal power is directly 
proportional to p, while for extreme limiting (A = 2) saturation is observed: the output 
signal strength is independent (~p") of the carrier. The amount of noise is assumed to 
be constant in both instances. 

On the other hand, the threshold behavior for weak signals (p < 1) shows that the 
signal output of the discriminator is proportional to p*, whether or not there is limiting. 
Thus, if the carrier is sufficiently weak, the limiter-discriminator combination behaves 
like any (half-wave) second detector acting on an A. wave,"' as far as the dependence 
on carricr strength is concerned. Similar remarks apply for the d-c output (Fig. 2), 


The integration is readily accomplished by representing exp iRoz) as the sum of two Bessel 
functions and evaluating the resulting Weber-Schafheitlin integrals with the help of Eq. (2), p. 401, 
Watson, Theory of Bessel functions (Macmillan, 1945). See also D. M. I., Appendix ITI. 

1D). Middleton, Rectification of a sinusoidally modulated carrier in the presence of noise, Proc. I. R. E., 
36, 1467 (1948 
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provided that we replace the mean-square modulation by the square of the mean. 
Accordingly, when the modulation possesses no d-c component, there will be no d-c 
output, .[ef. (2.36), (2.37)]. From these equations and from physical considerations it 
is also obvious that if the modulation is slow enough, n = 0, there will be no observable 


signal. 
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Fic. 2. Output signal power for no limiting 1) and extreme limiting (A 2). 
2. Noise Output Power). The low-frequency noise power Trom. th discriminator 


is independent of the carrier when there is no limiting, provided the carrier is unmodu- 
lated [(E« 2.34 It is also true that the outpul loes not depend on the modulating 
verified from (2.34) and (2.35). However, unlike 


frequencies in either case, as can be 
the analogou situation in AJ ect ption (mentioned above in connection with the 
sional) the magnitude of the noise power emanating from the discriminator dov depend 


on the spectral distribution of the noise, i.e., on thy /F-limiter-discriminator frequency 
response. The dependence is not very critical, covering a range of about 20-50 per cent 


lar and “optical,” or single-tuned responses ol the 
i i 


between the extremes of rectangulat 
same energy (see Appendix V, and ‘Table 1 of ref. 12 
Modulation suppression, ra phenomenon whereby the signal uppresses the noise 


’ ceouay tans 
when the signal is large enough, occurs in FM as well as in AM reception. For no 


limiting \ | the low-frequency noise output is independent of carrier str ngth 

when p” > | and is still random noise, though mixed with components of the modulation 

Kq. 2.40) et seq.| A like behavior occurs in extreme limiting (A 2), except that 
-D. Middleton, Spurious signals caused by noise in trigge ed circuits, J. Appl Phys. 19, 817 (1948). 

sual, aural, and mete C- 


133, H. Van Vleck and D. Middleton, A theoretical comparison of the vi 
lsed a in. the presence of noise 5 Appl Phys. 17, 940 (1946). 
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now the noise is absolutely, rather than relatively, suppressed, proportional to p' 
The reverse situation arises when the carrier is weak relative to the noise. The noise 
is then dominant, and it is the signal which appears in the output as a perturbation of 
the noise. (See See. 1 above.) 

3. Noise Output (Spectrum). It should be emphasized that discrimination, with or 
without limiting, is noé a linear process, even though the noise output following it may 
be the linear sum of the effects of the transient pulses (and their derivatives) which 
make up the incoming noise wave. In general, from the spectral point of view the noise 
output after the nonlinear operations of limiting and discrimination will consist of three 
types of component, just as in conventional AM cases. These modulation products are 
produced by (a) beats between noise and noise (nxn), (b) beats between the signal 
and noise (sxn), and finally, (¢) a signal contribution (sxs), due to products generated 
among the signal harmonics. For weak carriers the output spectrum is primarily de- 
termined by the distribution of the noise before discrimination (on the assumption of 
in ideal discriminator), while for strong carriers the spectrum of the discriminator 
output is not so heavily dependent on the spectral shape of the input, but is modified 
ina way that denends markedly on the degree of limiting. Limiting spreads the spectrum 
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3. Low-frequency noise spectrum following discrimination, no carrier 


IG. 3 


buting the energy into the higher harmonics of the distorted wave obtained 


i 


1 the clipping process 
The significant. feature of no limiting (strong carrier) is that the spectral intensity 
[t has been implicitly assumed throughout that the limiter-discriminator filters are wide enough 


limited wave without appreciable distortion due to frequency selection. Thus, to consider 


peetra in the neighborhood of Q 5, for example, is to assume these filters are at least five times wider 
than the JF band. The region of principal interest, however, is in most cases Q < 2, so that when there 
limiting our model of a limiter-discriminator response several times as wide as the /F is not too tar 


with actual practice For no limiting this response need not be wider than the /F. 
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near the low-frequency end of the output does not as a rule vanish. For extreme limiting, 
however, the intensity is always very small near f = 0, unlike the analogous situation 
in AM reception. Moreover, the spectrum near f = 0 is proportional to f° and is roughy 
triangular, becoming more nearly rectangular as the carrier diminishes in power relative 
to the noise, in accordance with Crosby’s results.” Consequently, in broad-band FM 
the amount of noise passed in the video or audio stages is noticeably reduced when 
there is heavy limiting as compared with the instance of no limiting (Part [[1). Compare 
the spectra in Figs. 4, 5 and 9, for example. This shows the necessity for limiting if there 
is to be an improvement (over AM, and FM without limiting) in the signal-to-noise 
ratio when the carrier is strong.”’* Limiting, however, is not mandatory for weak signals, 
to which FM receivers respond like AM receivers in their dependence on carrier strength 
(see Sec. 1 above). For narrow-band F./ the dependence on spectral distribution is 
not critical, since it is the entire low-frequency spectrum, rather than a fraction of it, 
which is a measure of the interfering noise. Spectral shape is then unimportant. Figures 


6-11 illustrate various other special Cases. 
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Fig. 5.* The same as Fig. 4, except now the limiting is extreme. 


*Correction added in proof: The ordinate scale factor in Fig. 5 is 27'/2b9/1.85w, = 1.08 m'/2ho/wp . 
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Fia. 6.* Figure 5, for large carrier amplitudes. 
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*Correctton added in proof: For f < 0.3, the curve for p = 10 in Fig. 6 should lie between the curves 
for p = 5and p = 20. 
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Part II: Correlation Function of the Output (General Theory) 


The correlation function of the low-frequency output of the discriminator is 








Ro(t) = (E(to)E(to + O]av- = e(B,(R,)B,(R2) 61 9e]ev- ° (2.1) 
The bar indicates the statistical average over the random variables and [ ],,. denotes 
the average over the phases of the modulation, if any."° The mean power spectrum is 
given by the well-known Wiener’*-Khintchine’’ relation 
Wif)=4] RW coswtdt, w= 2nf, (2.2) 
“0 
with the transform 
R(t) = | W(f) cos wt df. (2.2a) 
Jo 
The input voltage to our limiter-discriminator combination may be written 
V(t) = V, + Vw = Ay cos @t + ¥) + >) (a, cosw, + b,sinw,t), (2.3) 
n=l 
where f, is the carrier frequency (in the 7F range) and 


— [ D,(t) dt, with D,(#) the modulation. (2.3a) 


The quantities a, and b, are random parameters with the properties 








a = bn = 0; Aubn = 0; Gidm = Dndbn = Wf, ASS, ; (2.4) 


w(f,) Af is the mean power dissipated in a unit resistance by the n-th component f, ,* 
and 6%, is the familiar Kronecker delta, such that 5), = 0, n # m, 6; = 1. Here Ao isa 
constant representing the peak amplitude of the frequency-modulated carrier. Since 
the disturbance is assumed to be narrow-band, let 


O=w tw and w,=w+t+w, (2.5) 


where w) is some constant (angular) frequency, in this case the resonant or central 
frequency of the JF-limiter-discriminator elements; w, is then a frequency relative to 
w, , and w, is simply the (fixed) deviation of the carrier from exact tuning (w = @»). 
We have from (2.3) and (2.5) 


Vy = > 2 fa, COS (wo + wi)t + b, sin (wo + w/)t] = V. cos wot — V, sin wot, (2.6) 


1’The subscripts 1 and 2 on current or voltage amplitudes throughout refer to these amplitudes at 
the initial time 4; = t and a later time t2 = t + t( t > 0), respectively. 

16N, Wiener, Acta Math. 55, 117 (1930). 

7A. Khintchine, Math. Ann. 109, 604 (1934). 

188, O. Rice, Bell Syst. Tech. J. 23, 282 (1944) and M. C. Wang and G. E. Uhlenbeck, Rev. Mod. 
Phys. 17, 323 (1945). 
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so that 


bo 
~] 
— 


V. = > (a, cos wit + b, sin w/t); yi = z. (a, sin wit — b, cos wit); (2 


1 


n=1 
V. and V, are the portions of the slowly varying part of the wave that are in phase 
with cos wot and sin wot, respectively. A similar treatment for the signal yields finally 


for the complete wave 


V(t) = [V. + Ao cos (wat + V)] cos wot — [V, + Ao sin (wat + V)] sin wol 
(2.8) 
R COS (Wol aa 6), 
so that 
: we Ve. + B\ 
R= [((V.+ a) + (V, + B))]”, and 0 : tan-"( ¥ | na 
( a) + V+ os nH 
2.9) 
n= 0. 2,2 
in which 

a = A, cos (wat + VW) and 8B = Av sin (wat + W). (2.9a) 

From (2.9) we obtain 
6 = [(V. + a)(V, + 8) — (V., + BV. + a))/R’. 2.10) 
Consequently, the low-frequency output of our discriminator is [ef. (1.11)] 

E(t) = y{(V. + a)(V. + B) — (V. + BV. + @)/R’, (2.11) 
where y, = «8 for \X = 1 (no limiting), and y, = 4«8R,/m for \ = 2 (superlimiting). 
The correlation function (2.1) is therefore given by the eightfold integral 

R(jy =r} dV] at | dV.WAVa , Vin st 
(2.12) 


X [(Var + ar)(Ver + B) — (Vas + B)(Ver + )) 
xX (V2 + a)( Viz + Bo) — (V.2 + Bo)(Ve2 + a)] RiR? : 


The subscript N indicates that only the average over the random variables V,, , --- 
V.2 has been performed, and not that over the phases of the modulation. The complete 
correlation function requires the additional average. The quantity W, is the joint 
probability density of the eight random variables V,, , --- , ae 

The evaluation of (2.12) will require the Fourier transform of W, , namely, its char- 
acteristic function F.(z; , Z2, --- , 2s ; 4), which is specifically” 


19The subscripts 1, 2, --- , 8 here clearly do not have the time significance of those on Va , Vee, --* ; 


cf. footnote 15. 
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F(z, 22, °** , 23 3 t) = [exp {12,V.. + t22Veo + 123V.1 + t2,Vi2 + iz,V-; 
+ t2gVeo + t27Vir + tee V2} ]stat. av. 


bs 


| 


I 


bo 2 2 2 2 2 2 2 € 
exp - A (2; + 22 + 23 + 2%) (25 + 2 +2 + 2%) (2.13) 


\ 


bas o( t) (ZZ. ao 2x24) = o,(t) (2:26 — 22%5 


2323 — 2427) + oo(t)(Z52% + zd, 


where 


i) 


b, = | w(f)(w — wo)” df, and 410) = o w(f) cos (w — wo)t df. (2.14) 
‘ C “0 





Since w(f) is the mean input power spectrum, determined by the shape of the JF filter 
response, by is the mean input noise power and ¢(¢) is the correlation function associated 
with it, by virtue of (2.2a). Observe that ¢.(0) = bo ; then, for later use we find it con- 
venient to write ¢,(¢) in semi-normalized form, viz. 


¢,,(t)/bo = 7, (2), and Ta(b)max > Tn-1(b)max 2 **°* & Tolt)max = 1, (2.15) 
as is easily shown. The determination of W, and F, is given in Appendix I. Our result 
2.13) is somewhat specialized in that the input spectrum w(f) is assumed to be sym- 
metrical about fo . 

Now, from the definition of the characteristic function and its Fourier transform 
relation with WW, , the probability density may be expressed as 


WAV. , Ver, -**, Vie 30 = (2s) [ dz, 


. | dz. exp {—2(z, i + a - =~ 4 23 V2) } Fo(z, ofa, *** 2s ; t). (2.16) 
The correlation function (2.12) is then 
RAdy = AUS — 1 — 1 + 1), (2.17) 
where 


I = [ dV. °°: [ dV,2(2m)~* . dz, 


@ 


| dz,F (2, ge es OG (Vier ees V2) (2 18) 


x exp {—-1(2,V.. + ++: + 2eV..)}RPR2 - (j = 1, 2, 3, 4) 
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with : ; , ; 
= (Vir + a)(Vir + Bi) (Ven + a2)(Vi2 + 8:2), 


G, = (F a + a)(Vir + B)(Vie + B2)(Ver + a2), 
(2.19) 
Gs — (Vis + B,)(Ve: + a)(Veo + a2)(Vi2 + Bo), 


= (Vir + Bi)(Ver + o1)(Vi2 + B2)(Veo + on). 


At this point we follow a procedure suggested by some recent work of Rice. 
tegrate first over the random variables V,, and V,, , V.2 and V,, , taken in pairs, fol- 
lowing which we perform the integration over the V. We need the value of 


* We in- 


; oe on - 
K, = dx | dy —-——,, exp {—tzz — tb} 
je due + Y) 





and (2.20) 
4 . y 
K, = [ dx dy ———x,, exp {—taz — tyé}. 
seas I. (x? + y’) 
Transformation to polar coordinates and integration give finally 
— riz 
K, = oe 
rw? + 

and (2.21) 

— a 2> R(Qu — 2) > 1. 





P(ua**(? + By" 
For superlimiting d is 2 and yu is unity, while in the other extreme of no limiting d is 1 
and yu is 1/2. For this latter case the integrals (2.20) do not converge, but the results 
(2.21) still apply, as Appendix II demonstrates. Now, letting tc = a + Va , 
y = 8, + V.,, --- and using (2.20), we obtain the twelve-fold integral 


I” = i av. i ee ¥ dV ,o(2n)~* [ de, 


| | dz, exp {—i(esVer + o-°: + ZaVuz)}Fo(z: , -°* 283 8) (2.22) 


(—2mi)*2,2¢ exp {t(z1a1 + 222 + 28: + 2482)} (Vir + 8:)(Vi2 + Be); 
@ +2)" +a" te 
| eg and IS are found from (2. 22) if we replace 2:22(Vir + Bs Via + Bs) by 


a 2a(Var + 8 (Wer + a2), zota(Ver + a) (Vie + B.), Zsta(Ver + a)(Veo + a2) respec- 
tively. The integrations with respect to Ver, --° , Veo are effected with the help of 


~ | dY [ exp {—7zY}g(z) dz = g(0), 





xX 


and (2.23) 


a oe - = -| 2] 
= | Yay / exp {—izY}g(z) dz = 22 at 
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The first relation follows from Fourier’s integral theorem and the second also, on in- 


tegration by parts; in both instances it is assumed that g(+ ©) vanishes properly. Ap- 
plication of (2.23) in (2.22) for I{” yields the fourfold integral 


I” = —(2n)? [ re [ i. 





222 {Pieata — bo + 1Bidizs 93 1Bob ize + B:Bo} (2 24) 


x 
(i + 3)" + 4)” 


bo 2 2 2 2 * 
xX exp fo. 1 + 22 + 2s + 24) — Goler22 + 2524) + torts + ante + Bits + azo}. 


The other J{” are found from (2.24) if we replace its numerator, excluding the expo- 
nential, with 

Zrtalbizita + Biber — iaabrts + afr}, 

2223 {Pie2zs + tardi2s; — iBopiz2 + anBo}, 
and Zala{ (bio — bo) + tarde, — tardiz, + ayer}, 


respectively for j = 2, 3, 4. Returning to (2.17) with these expressions for I}, adding 
and collecting terms give us finally 


2 


Riv = ro I dz, --: [ dz, {2(t)(z:z2 + 2324) 


+ o1(t)(e:24 — Zo%s)” — ibs (t)(er24 — 202s)(Cnts + 2%, — Bit, — Boer) 


> (B:B22122 + OyOkaks = OB i224 _ a ,Bokezs) } (21 + 25)? "(é + jai (2.25) 


xX exp ‘ ; bo(zi + 23 + 23 + z4) sues do(t)(z:z2 + 2324) + i(asz; + az. + Biz; + auz)}. 


° ° : . . 20 
At this point we make the obvious transformation to polar coordinates 
2; = p, cos 0, , Ze = pz COS 6, 


23 = Pa sin 6, ; 24> po sin 02 , 





20When A = 2 (superlimiting), we can follow an alternative method suggested by Rice‘ for the case 
(8 = 0, a’ = 6° = 0) in which the denominator [(z*, + 2*s)(2*s + 2%)]?~*/* is expressed as a pair of infinite 
integrals in (say) u and v, whose integrands have the form exp [—ubo(z*, + 2*s)/2] and exp [—vbo(z*2 
+ 22,) /2]. The integration over 2 --- z4may be achieved with the help of a principal-axis transformation, 
leaving a sum of double integrals, each of which may be further reduced to a single one after suitable 
changes of variable. The final integration, however, cannot be effected in closed form, except in special 
cases. When there is no limiting (A = 1), this technique fails, as it is not possible to eliminate the radicals 


[(z*1 + 2s) (2*2 + 2%) )*/2, 
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for which the Jacobian is simply p,p. . From (2.3) and (2.3a) we may define a phase 7 by 


7 =a t+vV, 
whence 
n=o+ WV, 
and 1 > Waly + Vv, ; = wal to + t) + V, ’ te > ty — to ’ (2.26) 


D 2 a2 /2 “2 52 /2 z . . 
and since (a;,2 + §;,2)'"" = A, and (aj.. + Bi.2)'”” = m.2d0, wecan write the correlation 
function in its general polar form 


Rov = v3 | pr exp (—bop;/2) do, 


/0 


A-2 “" dd 
ps? exp (—bop/2) dp. | 


oo, 


/0 


°* dO, | 
Xx [ _ exp {—@opip2 COS (82 — 8,;) + tAop, cos (6; — m) + tAope cos (82 — n2)} 


0 


xX {do cos (@. — 6,) + 1P1p2 sin® (0; — 0:) + tAog, sin (0, — 4:)[mp1 cos (0; — m) 


+ 12P2 cos (A. — m2)] — Aomne cos (6; — mm) cos (A. — m)}. (2.27) 
The integration over 6, and 6, is accomplished in straight-forward fashion with the 
aid of the Bessel expansions of the various trigonometric exponentials [ef. (1.4)] and 
the orthonormality relations for the sine and cosine. The details are available in Ap- 
pendix III. The correlation function becomes finally 


ro) ( rae 
R,(i)x = Yr pi pw ro" {| — 6.00 + (a3 k+1 COS (k + 1)(n2 — Hay 


k=0 m=0 


2 


+ Ayn. 1e-1; cos (kK — 1)(m — »)) 


€E (1)? 1 2 . kina 
: (Fs. cos k(nz — m) — = Ae ns1.z+2 COS (Kk + 2)(n2 — m) 


+ ,( t)” > 


1 2 ne , 
— 5 Aiom+i,jr-2) 008 (k — 2)(m2 — n)) 


(n 


bo |++ 


) n>) 0) (1) ) 
+ (2bop)'’"d, (2) “hl c om Be. ap Ae cass cig Biante 


X sin (k + 1)(m — n) 


— H, (ere (ie eis - | 5g 1,1k-21) sin | k-1 | (m2 = | 
a “ave I « (0 2 
+ 2bopm Ne cos k(n. — ni] BH + 5 (Ei oa.a01 — Hy om. ie-11) } 


2 


where p = A}/2b, is the ratio of the input mean-square carrier to mean-square noise 








i) 
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voltage. The amplitude functions Hi amines 9 = kK —1i,k + I; qn = 0); ¢q = k, 
k+2,'k — 21, (n = 1) are obtained from 
_ i” k)/2 


Hi sane: X) = — ae - 
Qem-!)/20 mim + k)t]'” 


x | p™"**** exp [— bop /2)J,(Aop) dp 


_ hk +2m+n+ q+ A — D/2)lp” (2.29) 
qiby" P7228" “mim + kt? pe 


x i'i({k + 2m + n + q + A - 1} 239 a Ll: —p), 


where k, m, n, g are integral and positive. The integration is performed with the help 


of”? 


‘ , Ip ; iy 2 > 
| Jax) exp (—b'a)a* dx = ‘ v + u)/2! (#) iF iv + w)/2;v + 1; —a’/4b), 
d Dh’ Tp 20 
2b* T(v + 1) (2.30) 


Riv + “) > 0 


in which ,/’, is a confluent hypergeometric function. Recurrence relations for the various 
H’s are available in Appendix IV.*? The complete correlation function R(t) follows 
after the average over the phases#of the modulation has been taken, viz. 


r 


’ We ie 
R,(t) = T | Ro(t)w Atimoa) « (2.31) 
‘The mean power spectrum may be obtained in the usual way [ef. (2.2)], by determining 
the Fourier transform of P,(2). Specific cases are considered in Part IIT. 

\s in the analogous situation for amplitude modulation,”’ the output will consist 


of three sets of contributions: (1) (mxm) noise generated by the beating of the input 
noise components with each other, (2) (sxn) noise produced by the cross-modulation of 
the signal and the noise, and finally, (3), the discrete signal components (sxs). In our 
expression (2.28) for the correlation function, only the quantities in the first bracket 
for which g = 0 [ef. (2.29)] represent (nxn) terms; the remaining ones of the first 
and all those in the second and third brackets [| ] represent (sxm) noise, except in the 
latter when m = k = O; these are signal components: 
RAO vane Vi2bopl ni nelav-Hoo 
(2.32) 


‘ Qi bop oo { [wa + Dol to) [wa + Do(to + t) |} av- 


2G. N. Watson, Theory of bessel functions, Cambridge University Press, 1945. 


In its present form the series (2.28), or its Fourier transform, does not converge very rapidly 
when p < 1, so that it is more convenient to expand the amplitude functions H(p; A) and collect 
coefficients of p" (n = 0, 1, 2, ---) before summing. Only a few such terms are then necessary when 


p is small, say 0.3 or less; interpolation between these values and those for p > 1 yields values of the 
-orrelation function and spectrum for the intermediate signal-to-noise ratios: 0.38 < p < 1. 
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from (2.3a) and (2.28). Observe that the modulation is received without distortion since 
our idealized discriminator responds linearly to all frequencies. 

The mean low-frequency power, Wo , may also be found from the correlation function 
on setting ¢ = 0, according to a well-known theorem (D.M.I.). Now ¢,(0) vanishes, 
(ef. (2.14)], and n. — m, = 0, so that (2.28) reduces to 


W, a R,(0) = [Eo( to) Jav- — y. Zz z. {—6.(0) 2 


k=0 m=0 


x Cr + Fi on. 10-11) + Zbop{ [wa + Dg(to)]? } av- 
(2.33) 


x AE sans = | ae — Riscincel 


of which 2bopy,Hoo. {[we + Do(to)]*},,- is the signal power; the rest is attributable to 


the noise. A direct calculation of [Eo(to)*],.. = K?(B, (R)? 6],,. using the first-order 
distribution density W,(R, @) gives us*’”* for no limiting, 
1 = 


e 


(Wo)rnr = «B'S —¢,(0) + 2bop{ 1 — —— {fue + Dol(to)]*}av-?, (2.34) 
2p 





and for superlimiting 


eo 
(Wo) = 6 (2 ate +! 4,(0) + bop { [wa + Do(to)]”} av-) 
wr \2bo \ 
* 
7 me / _ Tf _ Re.) - 9 _nt2 2 
x (= 0.3444 Bi( ab,)) + > [(2n + 3)p"™"{ [ws + Dolto)]"}av- (2.35) 
— (n + 2)p.(0)p"*")/(n + 1)(n + 2h Ro K 2b, , 
where Ei(—x) = —fy e” dy/y. The exponential integral is tabulated in Jahnke and 


Emde’s Tables of Functions,” (pp. 6-9). Note from (2.34) that for no carrier or for an 
unmodulated one that is on tune, i.e., w, = 0, the mean noise power output (— — ¢,(0)) 
is independent of the carrier when there is no limiting. Then the spectral distribution 
of the noise, whatever its shape (for different values of p, w, = 0) will always enclose 
the same area. Equation (2.35) shows us that such is not true when extreme limiting 
exists (except, of course when p = 0). It is also evident that the total output power 
is independent of the modulation frequencies, since [Do(to)*],v- does not depend on 
frequency. 

The d-c output is easily found from our result (2.28) if we allow ¢ to become infinite. 
There is then no correlation between the wave at time t, = ¢ and at a later time t. 
Since ¢,(t), » > 0, vanishes in the limit t ~@, only the term m = k = 0 independent 
of ¢, and ¢, will contribute to give us the square of the d-c component. Taking the 
square-root, for no limiting we have (A = 1) 





2.), Middleton, On theoretical signal-to-noise ratios in FM receivers: a comparison with amplitude 


modulation, to appear in J. Appl. Phys. 
“4B, Jahnke and F. Emde, Tables of functions, Dover Press, New York, 1943. 








w\ 
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[Eo(to) lave = KBD? [awe + Do(to)]av-(/2)'*p,F (1/2; 2; —p), 

(2.36) 
on KBbo’ [wa + Do( to) lav-(w/2)' "pe" [To(p/2) + 1,(p/2)}.) 

in terms of modified Bessel functions (see Appendix IV). The superlimiting case (A = 2) 


yields 


Il- 


[Ko (t, = . ; KBRolwa + Do(to) lav-P iF’; 2; ——), 
vs 


| 
f (2.37) 
ae : KBRolwe + Dolto)jav-(1 — e°”). | 
T 


Both (2.36) and (2.37) are in precise agreement with previous results.*’* If 
[ws + Do(to)],y. vanishes, there is no steady output component. (Figure 2 also os 
the square of [E ‘o(to) Jeve » provided we replace {[ws + Do(to)]?}ar- by [wz + Do(to) |e -) 

We are now ready to examine the two important limiting cases which occur when 
the incoming carrier is strong, or weak, relative to the rms (input) noise level. Let us 
consider first the case when the carrier is strong, i.e., p> >> 1. The asymptotic expansion 
of ,F, (see D.M.I.) is needed here, namely, 


Fi(a;B; -—xz) = __T(6)_ aE o% ala ot + 1) 














r(B— a)” 
(2.38) 
, aa + Ie - 6+ D@-6+2 , | 
2!x” : 
1. Strong Carrier. For large ratios (p) we find accordingly that 
Ho ( ee r(k + Qm + n + q + Xr iain 5 ath eamaataaate 
k,2m+n,c\P3 A) = pint? 19/2Q 8") an Wm +4 ky? 
(2.39) 
x DL +k+2m+n+qtrA— IDK+2m+n—Ggt+rA—)/4pt-::) 
r[(q+3 —k — 2m —n — 2)/2] é 
The correlation function (2.28) becomes 
VT (A/2)?p 
R,(t) = = —¢,(t)[cos — bie 
ol) 2 TT(4 _ d)/2)° { 2( [ (n2 m)] 
(2.40) 


+ (A — 1¢,(é[(m + 1) sin (m — m))er 


+ (A — 1)*g0(t)['m2 cos (m2 — m)lav- + 2boplmmav-}, p> 1. 


The first three terms represent (sxn) noise, and the last term is the signal contribution 
(sxs); (nxn) noise is suppressed. Note that the contribution of the noise enters linearly 
in the correlation function: only ¢) and its derivatives appear. This means that the 
output noise is still random in the “pure” sense in which it entered the limiter-dis- 
criminator elements. In other words, the effects of the transients which gave rise to 
the original noise in the JF still add linearly in the output. In spite of this, discrimination 








148 DAVID MIDDLETON [Vol. VII, No, 2 
remains a nonlinear process, as the mixing of noise and signal indicates (Eq. 2.40). 
That the noise can still remain random in the present case is due to the fact that dis- 
crimination is (for our ideal discriminator) frequency selective and not amplitude 
selective: the noise loses its ‘‘randomness” (normal properties) only when subject to 
amplitude limitation, i.e., distortion and/or clipping. 


From (2.2) and with the help of (2.14) we may write the output noise spectrum 


when there is no limiting as 
p > 1, 


Wo(f)n_1 —4x"p° | ¢.(t)[cos (no — )]ay- COS wl at, p 


I§ 


) 


cos (nw, + w,)t cos wt dt | (w’ — wo)” (2.41) 


pn oD 
) “0 


= 4x’8° p> . | 


X w(f’) cos w’ — w)t df’, (f’ = w’/2nx), 


after expanding [cos (yn. — 7:)],,. in a Fourier series; here f,(=w,/27) is the fundamental 
frequency of the modulation. Integration using the Dirac delta-function” ~ gives us 


A. | (w’ — wo) w(f’)[ 5(w’ — Wo — NW, — Ws — w) 


IS 
of 
M 


+ 5’ — wo — nw, — w, + w)] df’ 


(2.42) 
= «8? Do Ajl(na. + ws + w)*wlwo + nw. + ws + w) 
n=0 
+ (nw, + w, — w) ww + nw, + w, — w)]. 
Observe that the spectral density at and near f = 0 is finite and nonvanishing, no 


matter what distribution the input spectrum w has. 
The case of superlimiting (A = 2) when the carrier is strong is particularly inter- 
tp and y = t, + #, the three noise terms in (2.40) may be written 


esting. If we let x = 

. > \ ( ( ny "2 

( cs "a \ | 

= ee - . “ 29°! bp '’ —¢,(t) | eos (wat + | D,(t) dt’) |. 
as T \2b, { os 


*y 


4 os(0{ + D(x) + Do(y)} sin (cet + | D,(t’) w).. 


+ bul + D,(x)]lwa + Doly))| cos (ut + | D(t6) w)}...} 
(2.43) 


“a oA _f ¢ 2\2 =) 
= a'( )p (- ot + 0xz/\ot + OY, G00) 


T 2bo ( 


x [ cos E + [ Dt’) at’ |az/x, 


) Pa 
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where 2%o(=2z/w,) is the period of the modulation and ¢,(t) = (0"/dt")¢o(t). Expansion 
of the integrand in a Fourier series, followed by the average over x, gives 


of 


Lo d/o 


ex 


cos E + | Dt) a | dz = ) A,, cos (nw, + wa)(y — 2), 
ea na (2.44) 


y—-xr=t, 


since the result must be even in ¢. The coefficients A, are the same as those for the 
development in (2.41). Next, we substitute (2.44) into (2.43), which shows that 


Ro(t)acisecae2) : ca Be) >> A,[—¢2(8) cos (nw, + w,)t 
” a (2.45) 


+ 2(nw, + wz)d,(t) sin (nw, + wa)t + (nw, + wi) do(t) cos (nw, + w,)é). 
Use of the integral form of ¢,(¢) gives us, finally, a relation similar to (2.41) for the 


noise spectrum: 


6 ova; ie 
W(f)rno 2 a x 3B” (& -)p ‘ is Ad cos (nw, + w,)t cos wt dt 
. T 2bo n=0 /0 


5 


4 | [(w’ — wo)” + (nw, + w,)}w(f’) cos w’ — wo)t df’ (2.46) 


“0 


a 


~—2 | (nw, + wy) Sin (Nw, + w,)t cos wt dt (w’ — w.)w(f’) sin (w’ — wo)t as’ 


“0 “0 


6 0 = , \2 ¢ 
Wolf r-2 = x K a( Ee ‘\p : : ad [” [(w’ — wo)” + (nw, + ws) — 2(nw, + w,) 
T 


X (w’ — wo) ]w(f’)[6(w’ — wo — Nw, — wa + w) + 5’ — wo — Now — ws — w)] as" 
which is 


6 “ _ , 
WoA(f)re2 = = kK “a(R. \p 1S A,w*[wlwo + nw, + ws + w) 
™ (2.47) 


+ ww + Nws + ws — w)], p > 1. 

The significant fact about (2.47) is that the spectral intensity for superlimiting always 
vanishes at f = 0 and is small in the vicinity of f = 0. Compared with our result for 
the pena d case le f. (2.42) Lt this shows why ~ a Sin ited aid carrier, any 
in a deal near zero ca ney, Ww hich is narrow valid to the w sidth (ys, b ) of the + eam 
[F-limiter-discriminator filter response, i.e., f K< f, . This result is identical with Blach- 
man’s,’ which was derived by a different method, and it agrees with the earlier work 
of Carson and Fry." 

From (2.40) we observe that for p* > 1 the noise is always suppressed, its contri- 
bution being proportional to p'* and that of the signal to p* *. Accordingly, when 
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there is no limiting, the noise power is constant and the output signal strength increases 
as p’, whereas for extreme limiting, the noise actually vanishes as p ©, and the signal 
reaches a limiting value independent of p. These comments are illustrated in Fig. 2. 
We notice also from (2.40) that when \ = 1, the noise power output is independent of 
modulation; this is not true for extreme limiting (A = 2). 

2. Weak Signals. When the signal is weak or, at most, of the same order of magnitude 
as the noise, the correlation function (2.28) is most conveniently described” by a power 
series in p; if we use (2.29) and omit terms in p"(n > 3), we obtain finally after con- 
siderable reduction the following result: 


R(t) = y2T(A/2)?2?b2{Gp(t) + pG,(d) + piG.(t) + 2{mne}ar-] + °°: } 
(2.48) 
p <1, 


where 
Go(t) = (r? — rote) oF (4/2, d/2; 2; 70), (2.49a) 


G.(d) = (ri - ray\-r oF'(A/2, 4/2 + 1; 2; 19) 


- To 7: oF (A /2 + iF /2 + ] ; 3; ro) [cos (Ne aoa nla} + {rin cos (no = m) lav: 
(2.49b) 


= ril(m + n2) sin (72 — m)Jav- — T2[CO8 (m2 — nlp Ps, d/2; 2; 70), 


and 


2 
G,(t) = “i - rane oF (A/2 + 1,4/2 + 1; 2575) + A/4(A/2 + 1)+2Fi(A/2, 4/2 + 2; 2570) 


2 
“7 + (4/2 + 1) oF,(A/2 + 1, A/2 + 2; 3; 70)[cos (n2 — m)]ev. 


»? . 
a 16 [cos 2(n2 = m:) Jav.[2 oF (A/2 > ie A/2 + ks 2; ro) aac off (A/2 + i, 
2 / 2 3X ee 
r 2 ao ] ; 3; an} a oF (X/ 2, /2 + | & 2; a{- 2 Tol m2 cos (ne _ m1) Jav. 


+ Ari[(n: + no) sin (ne a 71) lav + T2 \/2[cos (ne som 4 ndler} (2.49¢) 





+ oF,(d/2 + 1, 4/2 + 1;3; aye [(m + 12) sin 2(n2 — m)]ev- 





8 
22 ad 23 
+ _ [m2 cos 2(m2 — m)lav- — = [cos 2(m2 — nla} 


+ 2lmmlav-[2F(A/2, 4/2; 1; 73) — 1). 








Ww 
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The reduction is accomplished with the aid éf the recurrence relation 


2F,(a, b3c; z) — F(a, bse + 1; 2) = Fi(a+1,b6+ 1;¢ + 2;2). (2.50) 


cc + 1) 7 1) 
We see at once from (2.48) and (2.49) that the output noise is no longer “pure,” i.e., 
the transients producing the noise no longer add linearly in the output, unlike the case 
of a strong carrier and weak noise. 

The first term in (2.48) represents the correlation function for the noise output 
that would exist in the absence of a carrier, and 2p” [mM2lav- is the expression for the 
output signal or modulation. Jn threshold perception (p < 1) the observed signal is sup- 
pressed by the noise and the output signal-to-noise (power) ratio is therefore directly pro- 
portional to p” (for the leading term), just as in the amplitude modulated cases."""*** This 
is true whether or not there is limiting.*”° Figure 2 shows the output signal power as a 
function of the input carrier-to-noise power. 


Part III: Special Cases 


1. Noise alone. The first important example occurs when the demodulated wave is 
narrow-band noise. Here 7; = 72 = m = m2 = 0 and p = 0 also, so that the correlation 
function (2.48) is simply. 





o-Agh-3 Ye To L(m + 2/2)’ 
ToT2) do . 3 ’ p> = a + oe ’ 3 1) 


R(t) = y(ri sd 
=n(ri — Tof2) bo *2**F'(A/2)” oF'(A/2, 4/2; 2; To). 
In the instance of no limiting we obtain 
Ro(t)rer = F «'B bo(ri — rota) aFPy(1/2, 1/25 25 19), 


(3.2) 


b a 
xB 2 (ri — ror2)K(r), 


where K(r3) is an elliptic integral of the first kind, with modulus 7, . With the aid of 


. . n . LOTG -—a- fp 
2F’,(a, B; Y) 1) —_ Ty nn a)T(y es g)’ 





R(y) ¥ 0, —1, —2, Ry — a — 8) #0, -1, —2,-*-, 


we find the mean low-frequency power to be 


(Wo)rer = Ro(O)a-1 = —*°B'G.(0) = «’B" / : (w — wo)*w(f) df, (3.3) 





*See, in particular, Sec. VII, Part I and Sec. II, Part II, ref. 13. 
25See references 23 and 3 for a discussion of signal-to-noise ratios. 
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€ 


a special case of (2.33) or (2.34) when p = 0. The spectrum follows from (2.2) and (3.1); 


it is 
Wo(f)rar = thox’B”? > a | [rs()® — ro(t)ro(0)]r3"(0) cos wt dt. (3.4) 


At this point, and for all subsequent calculations, we assume a Gaussian spectral dis- 
tribution for the incoming noise, representing the composite /F-limiter-discriminator 
frequency response, viz., 


w(f) = Wo exp {—@ — wo)” /ws}, and bo = Wow,/2n'”’, (3.5) 


where W, is the maximum spectral intensity of the noise entering the limiter. From 
(2.14) and (2.15) we find that 

ro(t) = exp {—w,t /4}, r(t) = —w;tr,(t)/2, ro(t) = wi(w,t’ — 2)ro(t)/4, (3.6) 
so that (7; — ror2) = w,ro(t)”/2. Similar relations for the “‘optical’’ or single-tuned char- 
acteristic and for rectangular spectra are discussed briefly in Appendix V. Substitution 
of (3.6) into (3.4) gives 


ae a ate (1/2): ; 
Wilf). = | K Bw,bo —_,- —_—_=—5 exp {—2°/(2m + 2)}, 
— \2 2» Gay + ia Pi a 
(.4) 
where 2 = w/w, = f/f. . 
The distribution is shown in Fig. 3. 
The correlation function for extreme limiting assumes the form 
DG 26 (Bo ) ‘ 2 
Ro(br-20. = 8 K\- DT. — Toe) of (1, B32 Be 7) 
Oo} ’ ? 
T 2bo 
3.8) 
—16 . ( e (m — tt , 
= i a - 1bo§ —3—— _ ] log, (1 — 7). 
T Zbo/ “\ 1% = 
Although R,(0),,-2) is logarithmically divergent, the mean (low-frequency) output 


> 


power is actually not infinite, but zero, since (3.8) applies rigorously when t — 0 only 
if R>/2b, — 0 also; in other words, the limiting is so extreme that essentially no low- 
frequency components are produced. All the energy goes into higher harmonies.” For 
computation of spectra, however, it is sufficient that Ro/2b) < 1. 

We obtain the spectrum (A = 2) at once from (3.1) or (3.8) and (2.2) for the Gaussian 


distribution (3.5). We have 


a2 [hi exp {—°/(2m + 2)} 
‘ OL i 2 2 ) MZ é 4) ” 
Fis ele i, Fe 3.9) 
1 (27) 2b, aon (m + 1) 
Figure 3 also illustrates Wo(f),-2 as a function of 2. The spectral intensity at Q = 0 is 
Wo(O)r-2 = a i72 (5, BK webos (3/2), (3.10) 
T\ 27) \4aU00/ 


where ¢ is Riemann’s zeta-function™ and ¢(3/2) has the value 2.612. Some details of 


9 


the evaluation of (3.7) and (3.9) are available in Appendix V. From (A5-2) it also 


*6Jahnke and Emde, loc. cit., pp. 269-273. 
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follows that when Q becomes infinite, the spectral intensity (3.9) falls “o zero as 
Wol(f)x-2 = 168°x°w,Ro/e AQ > @). 

2. Tuned, unmodulated carrier. The unmodulated carrier in noise is the next more 
complicated case. Now m, = n. = 0 and m = ne = () also, provided the carrier is tuned 
to the center of the band, i.e., w, = 0. The correlation function (2.28) then becomes 


R(t) = > ea be ro(t)”” bg.) omnes SS 
(3.11) 
T :(0)"( Hi amet. — 3 are aa | call 


“a x4 ) 
with p Ao /2h, as before. The power output associated with the low-frequency terms 
is obtained at once from (2.34) and (2.35), viz: 


(W, x B’o,(0), 
(3.12) 


2 


: ‘ 16 (Ro \,0 ‘a rf Ro) _ . 
(Wo)r-2 = - ifs }3aoO i - bil 3%) — La + Die + 0) 


Observe that the mean-square low-frequency output of the ideal discriminator is inde- 
pendent of carrier strength when there is no limiting. This is what one would expect, 
since the carrier is tuned to the central frequency of the symmetrical band, and its 
d-c therefore vanishes. 

The mean power spectrum follows as before from (2.2) and (3.11). For the Gaussian 
input noise distribution it is 


‘ ‘ { ek ( ‘ 22° Je 2 
Hi y\2r wr bo§ > > (om +k+ om re 3 + 5 Ce 


0 22 0 


exp {—’/( (2m + k + 1)} 29° ) —n 
) we + b1 — ————— (3.13) 
. (2m + k + 1)” nf 2m+k+2 oan 
: Beas oa exp | —2°/(2m + k + 2) u | 
a — =H; — = Hi msi, 12-21) — aS 
(ME suot.n — 5 Mi cowes.ees — gM onesie ) ee as 


fhe spectra for no limiting are shown in Fig. 4 for various values of p, and for super- 
limiting in Figs. 5 and 6. The significance of the results is discussed briefly in Part I. 

When the carrier is large, the spectrum is found at once from (2.42) or (2.47). These 
results give the leading term, while a further use of the asymptotic development (2.39) 
gives the succeeding one. Here 


A, = & 


a —1 1/2 / 2 2 
and w(f) = 2bow, wr’ exp {—(w — wo) /ws}, so that 


(3.14) 


Wo(f)x-1 2 4 XB" bya) exp (— 2) — 5- E exp (— 0) 
<p 


an aon (1 + 9’) exp {-—2 2 | > 4s p p >i. 
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A similar procedure for superlimiting [cf. (2.47)] yields 


_ ’ 
Wo(f)a =2 _ 77 K 3° (Fos : {or exp (— 0) + ORS = [2 + ol hedind 7 (3.15) 


The mean output power corresponding to (3.12) when p* > 1 is simply —x«°8°¢.(0) 
when A = | and is 


\ 


9-7/4 
) bf r. o>). 298 
Pp \ Pp / 


Ro 


; ae ; 8 
(Wo)r-2 = | W o(f)reo df = 3 (KBr) {Bs 
T Io 


Jo 
We notice that when there is extreme limiting (Rj « 2b,) and the carrier is strong 
(p° > 1), suppression of the noise occurs, analogous to the similar situation in the re- 

. . . 11,13 
ception of amplitude-modulated waves. 
illustrated in Fig. 6. 


The spectra for the strong carrier case are 


The small-signal problem, on the other hand, gives us a different output spectrum. 
From (2.2), (2.48), and (2.49) the spectrum of the noise may be written 


Wf) = yiT(A/2)’2*""b5-* | [G,(t) + pG,(t) + pG.(t) + ---] cos wt dt, (3.17) 
which in the instance of the Gaussian distribution (3.5) becomes finally 
Wolf) = Wolf) noice?” + V2I'(A/2)*wye'/?2 703 


"4 ; 
J — (A, Pile 1), Se. {= 2 /\ (am + 2)} 
x a a” m'(m + 1)'(2m + 2)’ 


2(d/2)°(A/2 + 1), exp {—2°/(2m + 3)} 
mim + 2)'(2m + 3)" 


(n/2); [(: 20° 1 
~ mi (m + 1) 2m 4+ 1)” a A ~ Wm + 1] (2m + 1) 


— 1] exp {-2 am +1)}| 

[EEG AGEs 
F s\L4 \2 ~ 43 2/,,\2 ‘ 
1 (A\/A “(m+ 1 

+ 2 (3) (3 T i) (2 + 5) 











(3.18) 
x eB 1= 0 /Om + 2} 
m\(m + 1) "(2m + 2)'” 
9X 2)*(A/2 + 1) Ya? + Val d/2 + 2), fa em £3 
m'(m + 2)'(2m + 3)'” 
(A/2)(X/2)m(A/2 + 1) m ( 20° ) | ee 
—-— iF _- ——— ——_—- — xp {—2°/(2 1)} 
m\(m + 1)!(2m + 1)'” | ; 2m + 1/ (2m + 1) oF lists ame Eh 
(4/2)*(A/2 + 12 ( o _ | \ 
— — — xp {—2°/(2 2)} see 2, 
m!(m + 2)!(2m + 2) “78 m+ 1 — ¥/(am + 2)} | + J 
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The curves of Figs. 4 and 5 for p < 0.3 have been calculated with the help of (3.18).” 

3. Off-tune, unmodulated carrier. This is a generalization of the preceding case to 
include the effects of a constant deviation of the carrier from the /F-limiter-discriminator 
resonant frequency fy . It is now evident that 7; = waly , m2 = wa(to + ¢), and so m = 
M2 = w, . The correlation function (2.28) reduces to 


R(t) = v4 > > ie : C= cos (k + wal 


+ Ai om.ie-1; COS (kK — Lou!) 
+ ¢; = (a: oms1.z COS Kwyt — ; Hi. oms1.na2 COS (k + 2)wat 
— ; Hi. om+1.1e-2, C08 (k — 2)axt) (3.19) 


- 


+ (2bop)' ‘oa (1 = | Oe eae = Hy m+1,x) Sin (k + L)wyt 


we H bo. JA—3 ]) eae H, 2m+1,1k—-1 ) sin | k - l out 


+ 2b wip cos haut] SH k+1 a 5 (Es 20.0 — |) Pare » }}. 


The effects of the carrier with constant deviation are discussed in Part I. The contri- 
bution of the signal component is illustrated in Fig. 2, where it is seen (from the ordinate 
scale) that the signal power varies as the square of the deviation f, . The general determi- 
nation of the noise spectrum may be made in the usual way with the aid of (2.2) and 
(3.5). 

The example of a strong carrier is instructive. Then, with the Gaussian distribution 
as a particular example, we see from (2.42) and (2.47) that when there is no limiting 
the output spectrum is 

Wo(f)r-1 2 2m"? K’Brorbo[(Qs + Q)* exp {[—(Q, + 2)*} 
(3.20) 
+ (2, — 2)? exp {-(2, — 2)*}], Y= io 
6 


since here the Fourier development of [cos (nz — 7;)],,- is simply A, cos wt, A, = e.. 
The mean output power is (w#,«8)*b,/2, which is proportional to —¢,(0), and is inde- 
pendent of modulation, which in this case is the deviation of the carrier from resonance. 
This can be observed at once from (2.40) or directly by integration of (3.20). Figure 7 
shows the spectrum (3.20) for various ratios Q, . 

The spectrum after extreme limiting is found in a similar manner to be 


Wi(f)r-2 S ss Pd Ro sbop"[2? exp {—(Q, + 2)?} + @ exp {—(Q, — 2)*}], 
T 2b, 


(3.21) 
p >t. 
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Fic. 7. Low-frequency noise spectrum following discrimination, for a weak noise 
and an off-tune carrier and no limiting. 


Figure 8 illustrates the spectrum in this instance. The mean power output of the dis- 


criminator may be obtained from 


32 272 R; —} i 2 2 , Gs, 

(W, - Al — Ly Ol own |. {O. -L. 2)7} ete SE. wee SON) 1 ee 
Wo)rn2 372 k B 2b, w, bop j, 2 [exp {—(Q, + Q)°} + exp {—(2 )*} = 
3.22) 


) 


25 xfun)*(22-) 2 (0; + 
us 2b, Pp 


) 


showing that the noise depends on 2; , as can be seen from Fig. 8. 
The distribution in the neighborhood of f = 0 is easily determined with the help 


nol 


of Taylor’s series 


8 re) 20 1/2 = 
exp {—(Q, + 9)} = (2x) >> (+1)' Go” 920); 


o¢ = —_ 5 . — (exp {—3}), 
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Fic. 8. The same as Fig. 7, but for extreme limiting. 


a development which is valid for all Q@ and Q, . Equations (3.20) and (3.21) are accord- 


ingly modified to 


Wolfner oY 40/7 x7B8’e,b,{ 2 — 2003 + O7(1 — 503 + 205) + --- } exp {—93}, 
(3.24) 


<1 
and 


oa 64 2 —1 i “ 2 2 7 
Wo(f)r-2 oa « B’wsbop (H). Qo’ + --- } exp {-—Q3}. (3.25) 
0 
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A comparison of (3.24) and (3.25) bears out the general observation made following 
Eq. (2.47). 
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. 9. Typical low-frequency spectra for an off-tune carrier: (a) no limiting, 
(b) extreme limiting. 


4, On- and off-tune carriers modulated by a sine wave. In this case we have 


” =w, + D, cosw,t, and 


. Do do I 6 . 
" n=ot+ywsnat; p= Do = (2). . (3.26) 
Wa w,/ Q, 


Here yz is the modulation index, and D, is the maximum deviation of the modulation, 
while w, is its angular frequency. (For such modulations the inequalities Dp) K wo , 
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w, < w, must hold.)" Referring to the general correlation function (2.28) we see that 
the modulation enters by virtue of the coefficients 

[cos U(n2 — m)]av- 5 (im + m2) sin (qn, — m)Jav- 5 
and [m1 M2 cos [nz — m)lav- - 


A straightforward application of the techniques of Appendix III allows us to write 


these averages as 


feos lms — m)lay. = ) €,J,(lu)” cos (nw, + lw,)t, (3.27) 


n=O 


[((m + m2) sin Un. — m)Jav- = = én alu) J,-:(lp) sin [(n — lw, + le, Jt 


Pa . 
+ E J (Ip) + _ (Jnsi(lp) + J, (1) | in (nw, + lw,)t 


(3.28) 
dD, ; 
+ > Jnsi(lu) sin [(n + Io. + lwult?, 
= ) 
and 
. do d ); 
m2 cos Uns — m)lav- = rs etatn,| Pe J(u) + in J. | 
X cos [(n — 1)w, + lw, Je 
: doer, di : 
+ Ez + sae (Jnsi(du) + J, :(lp)) + ut (J n+2( Up) (3.29) 


+ J, (uw) | cos (nw, + lw,)t 


+ | Pe J nsi(te) + ae 1.4 | cos [(n + l)w, + lal) 
When the carrier is tuned to resonance, w, vanishes and (3.27) to (3.29) simplify con- 
siderably. However, a precise determination of the spectrum in either instance (w, = 0, 
or w, ~ 0) is a very tedious process. Fortunately, for many purposes the qualitative 
and semiquantitative conclusions obtained from the limiting behavior of the general 
case are sufficient; (see Part I). 

We note that for strong carriers (p’ > 1) the output spectrum may be found as 
before from (2.42) or (2.47). Equation (3.26) shows that the Fourier expansion of 
[cos (n> — m)]ay- is A, cos (nw, + w,)t, where A, = ¢,./,(u)*, so that for no limiting 
and our usual Gaussian noise distribution we find 


Wolf)rar & 20'? Barby >> FW (n0, + 2, + 2) exp {—(nQ, + Q, + 2)’} 
si (3.30) 
Wa 


, at 
+ (nQ, + Q — Q)° exp {-—(MQ, + Q — ay} fs 2=-, 


We 
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and the mean output power is simply («’6’w;)b,/2, [ef. (2.40)]. For superlimiting we 
obtain 
‘ fe) ee ee | 
Wo(f)re2 375 EB unbd( 5 \p 12 YS ,J2(u)sexp {—(nQ, + 2, + 2)} 
T | 


2bo =o 


+ exp {—(nQ, + 2, — 9)}?. 
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Fic. 10. Low-frequency noise spectra following discrimination, for a weak 
noise and a sinusoidally-modulated carrier on tune. 
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The mean output power associated with the noise in this instance is 
6 Ri 


(Welkus = | Wol(f)ra2 df ™ 6 (xBeo,)"bol 


is 
J, oD aws ? 3) ‘ 7S 


(3.32) 
+ z «(n° Qe + 2,4) Ja(u) f 


Figures 10 and 11 illustrate the variation in spectra for different modulation indices and 
modulation frequencies in the somewhat less general case of zero displacement from 
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Fic. 11. The same as Fig. 10, but for extreme limiting. 
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APPENDIX I: Second-order probability density, its characteristic function 
and associated moments 


Let X, (k = 1, 2,3, --- , 8) be s random variables, each of which in turn is the sum 
of a large number (N) of independent random variables z;, (i = 1, 2,3, --- , N), dis- 
tributed according to some law. Then the central limit theorem” states that the dis- 
tribution of the X’s is normal (i.e., multivariate Gaussian) in s dimensions, subject 
to certain restrictions [see e.g. S. O. Rice,”’ Eqs. (2.10-3)] on the distribution of the 
x;,’8 comprising X, , restrictions which are assumed to be satisfied here. The probability 


density of the X,’s is 


exp {— XM X 2} 
[(Qm)° ub - 





WAX,,°°*,X,;0 = 
(Al.1) 


kl 


172 ©XP E 5 > > ie a xX, |, 
k= | 


(2r)*" | i=1 | 


where M is the reciprocal matrix to u, with elements y |u{; uw’ is the cofactor of 
Me: ; | # | iS the determinant of u, and X is the transpose of X. The matrix p is symmetrical 


and has elements 
‘ 
te = 2, taku » (A1.2) 


in which z,, is the 7-th element in the sum comprising X, and 2;, is a similar element 
in X, . We assume here, without loss of generality, that po % = 0, i.e., that X, = 0. 


All odd moments are therefore zero. The Fourier transform of W,(X,, --- , X, ; é) is 
known as-the characteristic function and is” 


. lac , 
F(z, ,+:* ,2 30) = exp E = is os mits | (A1.3) 
= k=1 Jl=1 


Now, in particular, we denote here the various “in-phase” and ‘out-of-phase’ 


components V., , -:: , V.2 and their derivatives V., , --- , V2 [which by (2.6), (2.7) 
and the properties (2.4) are random variables] by X, , --- , Xs , (s = 8). The subscript 
lon V., V., or V,, V, refers to an initial time and the subscript 2 to a time ¢ later. 


The elements of our matrix uy may be computed as follows, from (2.6) and (2.7) in 
(A1.2). For example, in the instance of »,, we obtain 


f N 


tn = V., V se z (a, COS w/ to + b,, sin w/ to) 


(Al.4) 


\ 


N : . 
*< Z. (a. cos wt, + 6b, sin wile) fetat- aes = ; 2 (a, + b,)/2—- | w(f) df, 
1 J n=1 “0 


with the aid of the characteristic properties (2.4) on passing to the limit N —o. It is 


27See S. O. Rice, Bell Syst. Tech. J. 23, 282 (1944), see. 2.10. D. M. I. also contains a discussion 


of some of these questions—Appendix I. 





LY) 
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assumed throughout that w(f) is such that the various integrals involving w(f) con- 


verge. All 64 moments (u“’ = uw“, wx: = wie) are similarly found to be 


tig tx m = V*, = ws; = Vi, = pa = Vie = bo; 
uss = V2, = tee = Vo = wir = Vi = pes = Vin = be ; 
ae ts OE, ie ot Ne «ee a, oe ee Fa wh 
Mu = | | = poy = VV, pita VaVe 
= et ON te Cin © ha Bate ee 
(A1.5) 
tj = Vad = pea = Vali = AO; 
tia = VeaarVeo = —bes = —Vi2Var = rol); 
tis = VeiVeo = uss = ViiVi2 = —wso = —VaVes = —un = — V.iVie = ¢,(0); 
Mis | V.. = iw = Yes V., = —-f, = — V.iVes ~~ V..Ve = Ai(0); 
ss = Vial = try = V4 = —¢,(t); 
tse = Vel —per = —V.2V., = —),(2). 
Here we have 
P ere . 3” ; ; 
D, = | (w — w)"w(f) df, ¢,() = — | w(f) cos (w — wt df, 
“0 dt” “0 
(A1.6) 


ae a , 
A,() = — | w(f) sin (w — wy)t df. 
ot" Jo 


When the input spectrum w(f) is symmetrical, \,(¢) vanishes, as does also 


bo;., (7 = 0, 1, 2, ---). The probability density (A1.1) in such cases is given by 


+ a, 


; ( 
W(V.., °°» Vio) = (Qe) "* | aw * exp | - = 4a VS, + Vi. + Vi, + Viz) 
“| bh 


\ 


b Qu2(V Ver + Ver Vor) + Qu! (VerVer + VerVer — VeoVeo — Via V2) 


-|- Qu(Ver Ver + Vir Vez ea Vi2Ver st va¥.)}]. 
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rn . . kl 
Che determinant | » | and cofactors nu" are 


u | = [(bob. — $1)” + (bi — doe)” — (Ot + b5¢2 + b305)]° 


uw — = | w |'”*[bo(b: — 2) — boi], vw = |u|’ "[—dibs — Go(b2 — o2)], 

| (A1.8) 
. ee Te “ilgobe + bode], i — i.e |'”“dildobe + gob2 — $i], r. 
w= | pw |'[bo(bo — bo) — boi], =e = | w |'"[ba(bo — 0) + dido], | 


The characteristic function (A1.3) is represented by (2.13). When the spectrum w(f) is 
not symmetrical, the resulting expressions for W, and F, are correspondingly more 
involved, containing A,(t), (n = 0, 1, 2). Specifically, the characteristic function is now 


2 


1 a “ ‘ . » ' 2 
F(z, , °°: , 2,3) = exp (. 1 ie + 2. +23 +23) + doles + 25 + 27 + 232) 


+ Wolt) (2:22 + 2324) + 261 ( t) (2:26 — 225 + 2s%3 — 2427) 
(A1.9) 


— Wo(t)(Zs26 + 272s) + QWolt)(Z1%4 — Zo2%3) + 2ri(O(Zie%5 + Zeer — Zs2%q — 242s) 


— 2r.(t)(225 — fers 
Le 


APPENDIX II: Evaluation of K, and K. when R,—- 


The integrals K, and K, (Eq. 2.20) apparently are divergent when there is no 
limiting (u = 1/2 or \ = 1, and R, ~~), and it is not immediately evident that the 
results of integration (2.21) still apply in this case. The difficulty may be removed if 
we return to the general form (1.11) for arbitrary clipping. Our expressions for K, and 


K, accordingly become 


, 2 . 4 1 — ex ; iia Rod} ) 22 
KK. = = - a I (- es a Uo Jan iT] | jes 
(A2.1) 
tJ ((a° + y")’ » exp {—1ixz — tyé} 
x + y zP i _— Ys} 
and 
, 21 ° f ‘| — exp }—7zR,A} \ an@ 
k, = —— lim | |— ————_— )ar dx dy 
‘ T (Ro) Se | r || : 
(A2.2) 


2 21/2 
x ye + y) ) exp {—ixz — iyé}. 
x+y 
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Transforming to polar coordinates gives us for K, 
: a F (1 — exp {—7R dA} 
i= = lim | -—- — 

BEE 


WT (Ro 


* jan / dr J,(dr) 
(A2.3) 


4 | cos 6 exp {tZr cos (6 — w)} dé, 


where 


cos W¥ = => sin y = and Z=(24+8)'" >0. 


g 
. adh 173) > oreo 2 By 
(f+) (2° + &)'” 


We employ next Eq. (1.4), so that (A2.3) is found to be 








, 4z ,. r (1 — exp {—7RA} ® wee 
Kk, == Z _ I fi: B{—iedl J (rrA) J (Zr) ar| dn. (A2.4) 
The integral in r is a special case of Weber and Schafheitlin’s integral'® and has the 
values 
‘ — = 
| J\(rA)J\(rZ) dr = a2 2F (1/2, 3/2;2;0/Z’), Z2rX20) 
(A2.5) 
= =3'2F,(1/2, 3/2;2;Z°/N’), AS>Z>R>O!}. 
2r 
At Z = the functions are identical, viz., 2P'(0)/z, but logarithmically divergent. 


Since this divergence is only logarithmic, the integral in ) is still convergent, including 
the point Z = \, and this independently of R, . We observe now that (A2.5) may be 
extended to negative values of \ if we replace \ by | | and multiply the expressions 
(A2.5) by (—1). The integrands containing the real part of (1 — exp [—7Rpd]) then 


eancel, and we are left with 


, 4a. 1 f’ sn RA ,, , 2 72 
K, =_-_ Z jim E / — “ sf (1 2, 3 = 2; DN Z ) dr 
(A2.6) 
. | amar pes ce oe 
+Z] =" .F\(1/2, 3/2; 2; Z°/r) an| 
JZ nN 
Utilizing the fact that lim,,,..) (sin Rod)/A is a delta-function,”’* with the value 
76(\ — 0), we have finally 
; 2riz 2riz ™ 
¢ = Vi a (2” fy" (A2.7) 


the second integral in (A2.6) vanishes because of the property J2 f(x)6(x — 0) dx = 0, 
b > a > 0. A corresponding expression for K, is obtained on replacing z by £ in the nume- 
rator of (A2.7). If we set u = 1/2 in (2.21) we see that the result is precisely that given 


above, cf. (A2.7). 
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APPENDIX III: Integration for R(t), 
The Bessel expansions for the trigonometric exponentials in (2.27) are 


exp {—dopip2 cos (6. — 6,) + tAop, cos (0, — mm) + tAope cos (82 — 72)} 
= + z. 2. (—1)*e,0 (bop: p2)t "€r€nd 1(Aop:) Jd n( A ope) (A3.1) 


cos k(@. — 6,) cos (0; — ,) cosn(@. — m2). 


The integration over 6, and @, indicated in (2.27) is accomplished with the help of 


inde i LS] eae /.: F-cos ¥ 
( in yk-sin vm) ay oi ie al Z | wa pk-co ym) ay =0. (A32) 


2x /,  \eos wk-cos ym 2r Jo cos wk-sin wm 


We at length obtain the following five expressions for the average over 6, and @, : 


Oh COS (NN2 — ln,) +lek+1 elk a Q « 
CC! =- — : m [a ‘sh i % I 6!" My (A3.3) 
2e:€ 
Ol COS (NN. — ln) bak E chats A get csace ‘ 
Coin = . 6:6, — 5 6; 5, 9 9 6, , (A3.4) 
v(k l ce hk lf ek +2 k ° 
Com = = 6 (6; ° + 6;) cos 4, sin (nn. — In) 
deve 
+ 6.°°(6;°° — 6,) sin nm COS (nn. — In.) — 6, by ‘I+! sin nn. cos (l — 1)n, (A3.5) 
eC Pe : eke .'s elk-11/-1 
+ 6 6; cos nn» sin (l — 1)n, — 6, sin nn2( 6; COS 7; COS ln, 
— 6, sin 7, sin In,) + 6,°~"' cos nn (6;-” cos n; sin In, + 6,°~*' sin , cos Ln,)], 
yik I ek+17 ek+2 ch ° 
Cit«k = 6; (6, + 6,) cos m2 sin (nn. — In) 
fe, 
— 6, '(6,°” — 6,) sin m2 cos (nq, — In) 
+ 6; °'6 ‘sin In, cos (n — 1). — 6," ''6, cos ln, sin (n — 1) (A3.6) 
+ 6; ' sin 1n,(6,"~''~* cos m2 cos nn2 — 6, SiN ne Sin NN) 


— 6; "| cos In,(6,-* cos no sin nn, + 6," ~' sin n2 COs nn2)], 
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and 


(k) l f + t— t+ \k— 
C$",, = —— [eos m2 cos nq, cos 7, cos In,(6n°* + 8,57") (8r7* + 8)" 


46,€; 
+ COS m2 COS Ny Sin n, sin In,(6e*' + 6)*""')(8;** — e* 

+ sin 72 sin ny, cos , cos In,(6,*' — 85~')(8;** + er 

+ sin 7. sin ny sin , sin In,(6,** — &")(6,° — 5, ') (A3.7) 
+ COs no Sin nn. cos m, sin 1n,(d,** + &§-")(8;°' + 8:7") 

+ 0s no Sin nq sin m cos In, (6,*' + 8-")(6;)""! — 87°") 

+ sin 72 cos nq. cos 7, sin 1n,(6'""! — 8k*")(8;*' + 8:7") 


+ sin 72 COS Nm SiN 7, COS In, (6,," a evar —_ 6:°')] 


with the convention that 5;° = —6; . The correlation function may now be written 


ao 


Rod = 1% DY Laee(— rr? fp *Si(a,Ao) exp (—bopi/2) dor 


i 0 n=O 


x | Pe “J n(p2Ao) exp (—bops 2)1(@opiP2) Ape (A3 8) 


v(h 2 ’ . 7 y(k . . y(k) s° * itt 
X {O01 + 1 pi pol oh + tAodimpiCs, in +- tA od nope “9 - AgmmCS:tn} - 
The series in n and I are eliminated because of the Kronecker deltas in (1, k) and (n, k), 


which appear everywhere in pairs in the various C’s. If we separate the variables p, 
and p, by expanding /,,(@op,p2), integrate, and reduce the C’s, we finally obtain (2.28). 


APPENDIX IV: Some properties of Hy .om+n..(p; ») 





The determination of spectra requires a wide range of values of the order 
(k, 2m + n, q) of the amplitude functions Hy 2n+n,.(p; 4). These values are best found 
with the help of recurrence relations, which greatly facilitate the work of computation. 
The amplitude function H is described by (Eq. 2.29) 


rik + 2m+n+q+A— )) 2)p”* 


Ap sasn.ctia = Ra “% : 
giber 72% —“lmi(m + kt" 


(A4.1) 
xX Fi({k + 2m+n+q+A — 1)/2;¢4+ 1; —p), 


k, m,n, q > 0 and integral. Here q is equal to k plus an integer, and n is either zero or 
unity. Using this fact we distinguish, for convenience, two classes of functions, denoted 
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respectively by superscript zero and superscript unity. The amplitude functions are 
thus only of double order (in k and m). 

Recurrence relations between the H’s are readily found from those of the confluent 
hypergeometric function ,F, , a table of which we list below,” after multiplication by 
the appropriate factors [ef. (A4.1)]. We have 


f F Fs Fs I 
| ‘ a 3 : B 2a + p 
2 ¥B +p(6 a Bla + p 
3 , l B (8 a l A4.2) 
| 8 p B 
2 a 3 .=4 l—a¥} 
6. +p(p a B(6 1) Bl 8 +p 


The subscripts on F heading the columns refer to an addition or subtraction of unity 
ina or fof F = ,F,(a; 8B; +p), the rows list the factors which multiply the heading 
F’s, and the upper sign wherever (+) occurs refers to ,F;(a; B; +p), while the lower 


refers to ,F,(a; 8; —p). The sum of each row is zero. 


There are five distinct H’s, corresponding tog = k + 1,|k — 1), with n = 0, and 
q=k,k+2,\ k —2|, with n = 1. For each amplitude function, two recurrence formulas 
are needed to cover the complete manifold (k, m): one for k fixed, m = 0, 1, +--+ , the 
other for m constant, k = 0, 1, 2, --- . Since m appears only in a of , F(a; B; —p) |see 


A4.1) above] the first row of (A4.2) may be used directly in the former case, viz: 
B — a), F,(a — 1; B; —p) + (2a — B — p), F(a; B; —p)]/a 
= , F(a + 1; B; —p). 443 


In the latter, where k is varied, both a and 8 contain k additively, and so we seek a 


relation between F, F,-3— , and F,.3, . Now if we change 8 to B — | in 4. of (A4.2), 
and again alter a to a +1, use 5. and eliminate Fs, and F,, from these three relations, 
we obtain 
WW , y t 

—(8 — 1),F,(@ — 1;B8 — 1; —p) + (6 — 1 + D),F (a; B; —p)} 


= F(a T 1: B+ i; Pp). 


28See Appendix 4 in ref. 7. 








~ 
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lhe general recurrence formulas follow from (A4.3) and (A4.4), and are 





f a = + 2mt+nu+a-—l el 
”? 2 m+ 2 4 k)1'”" \ > j 

( k+ 2m+n+qt+A- ‘) ~ _ Aivansn.e Fa 

\ 2 I ((m + 1)(m+ k + 1]'” 

2m +n+k +X — pyHe's, aso a Os tucnk dea 

ane 
| Zz (é +2mt+nt+q+t+ A- \) _ 

m+k-+ 2)” | 2 (m+k+1)'” 


HH, 


We observe from (A4.1) that a = (k + 2m + n+ q+ — 1)/2 is either 


integer or half-integer. When a is a half-integer we use ~ 


_— 2°T(u + 1) 
(+p)" 


exp {+p/2}/,(+p/2). 


Phis Is a ¢ 
and ¢ *J,(2).2" When a is integral, ,F 
in p ana ¢ 
In order tha 


The initial functions 


o give us H ; are 
Hl P; 4 ayp . Fy 2; = —Pp) ayp' “( 4 *(To(p 2) . 2 I\(p ai (XY = 
ap '’*(1 — e”) , v= 
H Pp iF(A/2 + 1; 2; —p) = ayp' ¢ [,(p 2) — I\(p 2)] 2 
7a 
= ap ¢ 
Aa) . , »/2 , 
H . { Pp Fi(A/2 + 1; 3; —p) = ae I,\(p 2 


= a,(l—e” — pe”)/p 


Ref. 21, pp. 191. 
W. R. Bennett, J. Acoust. Soc. Amer. 15, 164 (1944). 


? In general, four functions must be known initially for n = 0 or n 
hat (A4.5) and (A4.6) can be used to give us a complete set of H’s. 
*° to be used in the recurrence relations when n = 0 (q = k + 1) 


169 


(A4.5) 


(A4.6) 


a positive 


(A4.7) 


onvenient form for computation because of the available tables of e-*Jo(x) 
may be expressed in closed form as a polynomial 


=]. 


1) (A4.8a) 


2) (A4.8b) 


(A4.9a) 


(A4.9b) 


(A4.10a) 


(A4.10b) 
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170 
and 

0 Aa 
it pied \/2 + 1) Fr oe sp iF ,(A/2 + 2; 3; —p) 

= ae ”’ [pl o(p/2) p + 1)1,(p/2))/2 *s st 5 Hose 
\s 
Ag pt 9 

where a) r(A/2)/b 2 


in particular, aq=7T 


A similar technique is 


k =0,H may be 
are 
Hioo = Q 1F,(d/2: 1: p 
AQ) . : 
Hi20 = 53 FPi(A/2 + 1;1; —p 
0 AQ) : o 
HS) = 5 pF (4/2 + 1; 2; - 
and 
Aa) + 
H = (A/2+ 1 Wt Pp F(d 
_ Pp é 12/9 
= 4 36 (3/2 — p)I.(p 


= 2a.p’e "(1 — p/2)/6° 


We continue in the same way 
tions are 


HEB = 0 F,00/2; 1; —p) 


pho = 


applied to obtain H, 
obtained from (A4.8-A4.11) above. The initial functions here 


"/2 anda a, = (2b,) 


0, q / 


n= 


pP)Io(p/2) + pli(p/2) 


2+ 2:2 —= 2) 
: : ay 
2) + (p — 1/2)1,(p/2)] | = Bsa: 
Cr 


to -optam FH, 5..;:1(n = 
= ¢,¢ P “I o(p 2) 
=Ccr -. 
= 2 (1 — plolp/2) + plilp/2) 
= 5} (1 — p)Lo(p/2) + pli(p/2) 
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A4+.1 1a) 


(A4.11b) 


A4.12) 


Ge When 


A4.13a 


A4.13b 


2 (A4.14a 


(A4.14b) 


A4.15a) 


(A4.15b) 


(A4.16a) 


(A4.16b) 


1,g = k). The initial fune- 


(A4.17a) 


(A4.17b) 


(A4.18a) 


(A4.18b) 
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1/2 p/2 


Hi? = d/2qp'F,\(d/2 + 132; —p) =e, ~ (Io(p/2) — I,(p/2)] (A4.19a) 
- = c.p'e” (A4.19b) 
H® = (r/2 +1) x p'?,F(d/2 + 2; 2; —p) 
p Sivs— ((8/2 — plolp/2) + (» — 1/2)L(p/2)] = 3'7HS (A4.20a) 
Qcup (1 — p/2)/2'”, (A4.20b) 


and here 
ce, = T(A/2)/?2°”: specifically, c, = (a4/2bo)'” and c. = bo’. (A421) 


There remain Hfinsi.nc2 (n = 1,g = k + 2) and Hy anes, je-2), (n = 1,g =| k — 2)). 
For the former we find the initial functions to be 


Hoi2 = (4/4) p-iF'(A/2 + 1; 3; —p) = ce” I,(p/2) (A4.22a) 
=¢(1—e” — pe ”)/p, (A4.22b) 
H\: r/2 + 1) “ pF \(A/2 + 233; —p) =e, —>— [plo(p/2) — (p + DI(p/2)) 
(A4.23a) 
= Cpe”, (A4.23b) 
Hy, 2+ | - pF (4/2 + 2;4; —p) = “pp” ((4 + p)I,(p/2) — plo(p/2)| 
(A4.24a) 
= cp (2 — Qe” — 2pe” — pe”), 
(A4.24b) 
an 
> 2. 9) 2 ie : 
Hs = v3 : * ~ Ss iF(A/2 + 3; 4; —p) 
C1 /%9-”/1(Dn? + p)Io(p/2) — (2p° + 3p + 4)1i(p/2)] (A4.25a) 
9 
Cop 4 (A4.25b) 


We have finally, for si sti 


Aan Les Fy(X/2 + 1; 1; —p) = s~ [(1 — p)lo(p/2) + pli(p/2)| 
(A4.26a) 
C . 
==—75(1 — pre; (A4.26b) 











nN 
7 
> 
= 
7] A 
= &% 0 vo 
(Q - )9° - L 
_[19/ (9% — m)] - 
‘(Q — r)eim| | r (| oi —) dee LS. Sa a 5 L = ({)n poyd, (‘¢) 
0 Wr ; lL OPHEO0,, > 
io<7 ’ 
<< 2° ‘(| x | —) dxo *m— 
7 "(| x | —) dxa % 
~ {x . f>s+o 
> . \e us r soo © 
= e : I 2 z ‘ ; 
= — |°o - — 17 { ‘7 — Of = ippnbunjoay (*Z) 
a. F UIs — 3) is soo a) x us ae ? l Wd 
ze : Ps ' 1S 
— _ t (In F ’ 
> routs °/: f+ %> { wf — Of SOM as (fa 
= al f ‘ y F 
t U— ) dxo 
(Z/,«—) dxo = (¢/.7—) dxa : im—| (¢/ x—) dxa 1 20?/ (0 — m)—} dxo °yy = (f)m upissnps) (1) 
ed (, &8\G i ' mee 
(T — <2) 
4 — Ty = (2) (7)*u (7) (7) ° (74M = X) Woy aYNA)ouy unwoade 
2 } } A Iyhpouy ; 
N 
~ 





-~ 
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(1) r is i] 
Ho30 = (A/2 + 1) oa 2 iF, (A/2 -+- 2; i; — = Ls oo 3p + 3/2)1o(p/2) 
(A4.27a) 
+ (2p — p’)I,(p/2)] 
C2 . 9 2 
= “7a (2—4p+>p), (A4.27b) 
H® = (r/2 + 1) NOP F\(A/2 + 2; 2; —p) 
ee e?((3/2 — p)Io(p/2) + (p — 1/2)1,(p/2)] (A4.28a) 
2Cop e?(1 — p/2 (A4.28b) 
and 
H331 /2 + 2)(A/2 + 1) a pF \(A/2 + 3; 2; —p) 
: E — [(p? — 9p/2 + 15/4)Io(p/2) — (p? — Zp/2 + 3/4)I(p/2)]  (A4.29a) 
(A4.29b) 


cope”, » 
~t7s (pp — Gp + 9). 


APPENDIX V: Evaluation of Series (Noise alone) ; 
Remarks on Optical and Rectangular Spectra 


We note here some possible methods of summing the series (3.7) and (3.9) in the 
‘i . ° + . vr . 31 
instance of noise alone. For no limiting (A = 1) Kummer’s transformation” puts the 
series (3.7) in a form that is more rapidly convergent. Here we write 


A,, = 2m + 1) exp {8/(m+ 1)};) a, = (1 2)? exp {—B/(m + 1)}/(m!)*(m + yr": 


28 = w/w, 


and 


1/2 
2(m + 1) exp {8/(m + ifs — (mtiye (m + 1/2) |. 


(m + 1)” 
The series becomes 


ae + dy (l —,,)@,n = 2+ Th = dns; (A5.1) 


m=O 


31—. P. Adams, Smithsonian mathematical formulas and tables of elliptic functions, Publication No. 


2672, Smithsonian Institution, 1939. 
- 








? 
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Because lim,,.... \,, — 1, the series ‘= (1 — X,,)a,, converges more rapidly than pe a, 
alone. 

On the other hand, for superlimiting (A = 2), we may use Euler’s summation 
formula®’ to sum the remainder after m + 1 = M terms. Our series (3.9) becomes 


' M-1 
exp {—8/(m + 1)} exp ;}—B8/(m + 1)} 72 , 
—_ a5 = =55 (r/B) °° 8 ([B/M)"’") 
(m+ 1)°° 2 (m+ 1)" +r 


m=0 


+ M~°” exp {-8 uis[3 + a — of | (A5.2) 


1 [105 , 1058 — 218" 4 ) 
r aur | 8 * 4M 27 wit 


| 

p 

) 
The function 0 (x) = 2r°'”” {i exp (—y’) dy is the error-function.”' Observe that when 
8 — the series approaches the value (#/8)'”’. 


a * ok sa * K 


When there is only random noise entering the limiter-discriminator elements, the 
correlation function of the output [ef. (3.2) or (3.8)] in the two extremes of no limiting 
T. are 


and superlimiting is proportional to A(t) = [ri(t) — ro(t)re(t)], where ro , 7: , 
defined by (2.14) and (2.15). It is instructive to determine A(é) for the case of ‘‘optical’”’ 
and rectangular spectra, as well as for the Gaussian distribution used extensively in 
the present paper. The results appear in the table. Here 6(¢ — 0) is the familiar delta- 
function of Dirac, and the spectral ‘“width-factors” w, for the different distributions 
are related by 

(A5.3) 


1/2 
(We)Gauss = —173 (Wo)rect- = TF (Wo)opt- 5 
if the input noise power (b,) and the maximum spectral intensity W, are the same in 
these cases. 
When the input spectrum is “optical,” the output noise spectrum is found from 


(2.2) and (3.1) to be 


W(flooe- = ¥,00 2 T'(A/2)? ce | cos wl oI',(d/2, d/2; 2; 75) 6(t — 0) dt 
WT Jo 

(A5.4) 
Qe T(2—/2) Fe 
There is no d-c component, since A(~) = 0. The spectrum has a uniform intensity 
for all frequencies, and hence the power output is infinite. This catastrophe is never 
realized physically, because it is not possible to produce an exact optical distribution. 
Owing to parasitic effects in the circuit elements the falling-off of w(f) at the higher and 
lower frequencies, away from resonance, is always greater than the ideal optical be- 
havior, and so A(t) is finite for all values of t. The output spectrum, widely spread when 
parasitic effects are minimal, nevertheless encloses a finite area and hence finite power. 
The saturation of the discriminator in practical cases imposes a further deviation from 
the optical shape. 


Reference 31, Section 1.86. We follow a suggestion of Rice* here. 
* 
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CURRENTS ON THE SURFACE OF AN INFINITE CYLINDER 
EXCITED BY AN AXIAL SLOT* 


BY 
C. H. PAPAS anp RONOLD KING 


Cruft Laboratory, Harvard University 


1. Summary. The distribution of surface-current density is determined on an 
axially infinite, perfectly conducting cylinder. A constant distribution of electric field 
across a narrow axial slot is assumed given. 

2. Introduction. The analysis of the electromagnetic field radiated by a perfectly 
conducting circular cylinder of finite length fed by an axial slot becomes quite involved 


— 


Fic. A, Infinite cylinder with axial slot. Constant # across mouth of slot. 


especially when the end surfaces are rigorously treated. However, if the slot is situated 
half way between the ends of the cylinder, the peripheral surface current distribution and 
the radiation pattern in the plane through the slot and perpendicular to the axis of the 
cylinder is approximately the same as for an infinite cylinder with an infinitely long axial 
slot. Therefore, in so far as this peripheral surface current and radiation pattern is con- 
cerned, it is only necessary to solve the following problem. 

The problem investigated is the distribution of currents on the surface of an axially 
infinite cylinder of radius a fed from an axial slot across which is maintained a constant 
distribution of electric field defined by 

E = 6,E,(a, 0, Z) = 0,E,G(@)e ***, (1) 
where G(@) is defined as follows: 
1 for —0,/2 < 0 < &@/2 
G(@) = 
\O for | @| > 0/2 


and @, is the unit vector in the 6-direction of the cylindrical co-ordinates r, @, z. 

teceived April 15, 1948. The research reported in this document was made possible through support 
extended Cruft Laboratory, Harvard University, jointly by the Navy Department (Office of Naval 
Research) and the Signal Corps, U. S. Army, under Contract N5-ORI-76, T.O.1. 
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Since the cylinder has infinite conductivity, all the components of the /-field for 


r-=aand| @! > 6/2 are zero except the radial component EF, . Therefore, for r = a, E, 
9-- 80 ae. .. 
=+4 VQ 
\ e 0 a+ 2 
x , 
r=a 





Fic. B. Cross-section of cylinder with slot. 

is constant in the gap and zero outside it. The applied field, #, , maintains currents in the 
@-direction. The field outside the cylinder has only the components /, , E, , and H, so 
that it is TM (transverse magnetic) and may be regarded as a superposition of cylindrical 
waves traveling radially outward from the surface of the cylinder. At great distances 
from the cylinder, E, is negligible compared to E, , so that the field is predominantly 
TEM (transverse electromagnetic). As 7 approaches a, E, of the radiated field must 
approach the EZ, of the applied field, and H, becomes numerically equal to the 6-compo- 
nent of the surface-current density J, . 

3. Fourier expansion of the applied field. To obtain an analytical expression for 
the applied field as given by 1), it 1s expanded ina complex Fourier series; that is, 


E,G(@) « = Z. ic (2) 


where a, are the Fourier coefficients given by 


a, = E,(nr)™* sin (n6,/2) 3) 


sy substituting (3) in (2), the Fourier expansion for the applied field is 
E, = E,G(6) « * = Bow” Ss n' sin (n@,/2) ec” (4) 
This infinite series represents /, for all values of 6 between —7 and z. 


4. General expression for field surrounding cylinder. The most general expressions 
for E, and H, for r > a are the following:''” 


i, = —— > HH *(kryb. e"*"' (5) 
and 
H, = > Hi (eb, ** (6) 


where k = 27/xX, \ being the free space wavelength, and where H,."’ (A) are the Hankel 
functions of the first kind and order n, H‘"’(kr) are the derivatives [d/d(kr)]H<” (kr), 
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This equation determines b, for all values of n. Explicitly, 


2 1949 
or and b,, are constant coefficients. EZ, is not involved in the determination of the unknown 
- coefficients. Expressions (5) and (6) hold for all values r > a and for —7 < 6 < =. 
| 5. Matching at surface of cylinder. At the surface of the cylinder, r = a, expressions 
| 4) and (5) must be identically equal for all values of 6 between —z and 7. Therefore, 
Ie, (n7) : sin (n@,/2) = — (inw/k)H”’ (ka)b, . (7) 
b, = {ikE, sin (n0,/2)}/{pornH<”’ (ka) }. (8) 
By substituting (8) into (6), 7, becomes 
B™ the) s.0- ies 
2 mn (9) 


H, = (ikE,/pwr) D> n™ sin (nO/2 ae ( 
a... H,, °'(ka) 
at the surface of the evlinder. 
6. Surface-current density. Since J, is numerically equal to H, , the surface-current 
density around the periphery of the cylinder is 


es -. | 
Je = (tkE,/pwr) n~ sin (n6o/2) —~ > — € (10) 
>> HS” (ka) 


The limits of summation in (10) may be changed by using the formulas 


H(z) = e""* HS” (zx) 
(11) 


HO) "(¢) = e°**H,” (2) 


= 1/(1207) mhos, the following expression 


Substituting (11) in (10), noting that k/(uw) = 


is obtained for J, 


i [ 70, H;” (ka) - ~ sin (n6,/2) H\” (ka) | 
he 21 Sie a cos 10 (12 
1207 | 2 HS’ (ka) 2» n Hi” (ka) 12) 


Let a be the real part and 8 the imaginary part of the bracketed expression in (12) so that 


this he ‘omes 


‘ 2\-1m iw ° ’ 
Jo = (1202*)*Eye** ‘(a + iB) (13) 
For any given values of the parameters 6) and ka, a and 8 are functions of the angle @. 


The absolute value of J» is given by 


Jo | = (1202*)"'E,(a’* + B°)'” (14) 
7. Radiation pattern. The radiation pattern is the distribution of | KE, about the 
cylinder along a circle whose radius r is so large that it is in the far zone defined by 
kr >> 1. The formula for E, is obtained by substituting (8) in (5). The result,’ which is 
valid for kr > ka, is 
' bd sin (n0,/2) HS?’(kr) ine (15) 
te = al 8 
asa n H}"’' (ka) 


E, = x 'E,.ec ** 
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Since in (15) a very good approximation to the infinite series is obtained by summing 
from n = —M ton = +M where M = 40, it is permissible to use the following asympto- 
tic form when kr > 1 and kr > M. 


H"' (kr) ~ i(2 ‘kr)' 2,' kr 2n+1)/4) 4 (16) 
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Substituting (16) into (15) and changing the limits of summation, /, (for far zone) 
becomes * 
Ey, = «Be **'(2/akry ee [0,( 2H,’ (ka)} 
(17) 


+2 >> sin (n6,/2)e'""”* {nH (ka)}~' cos n0 
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The factor in front of the first square bracket is a complex constant independent of @. 
Denoting the complex quantity in the square brackets by 6 + 7, and taking the absolute 
value of both sides of (17), the following expression is obtained for £, which defines the 


radiation pattern 
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8. Computations. For suitably chosen values of 6, and ka, a and 8 were computed 
using tables of Hankel functions for orders 1 QOton 20. For n greater than 20, 


} 


use was made of the asvmptotic form 


Hi.” (ka) ka) (19) 
H (ka) vi P i 


holds for (ka)? «< 2n + 1. Forty terms were used in the calculation. 


whi h 
9. Results. The current density on the periphery of cylinders of various sizes 
ka 1/10. 1/2. 2.3. and 4) and for two sizes of gap (@ = 1/10 and 1 20) were calculated 
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using (13). The results are in Figs. 1-6. Figures 1, 2 and 3 show plots of a, 6 











and (a@° + 6°)’ as functions of the angle 6 for 0. = 1/10. Note that a and 8 are propor- 
tional, re spe tively, to the real and imaginary parts of the surface-current density J, ; 
a + B°)"*is a pe aes ul aes e J, |. Figures 4, 5, and 6 show corresponding 
plots of a, 8, and (a” + 6’)'” for 6 = 1/20. 
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For ka = 1/10, the current density about the periphery of the cylinder is almost 


constant. For ka = 1/2, 2, 3, and 4 there is a standing-wave pattern about the periphery. 
In all the cases considered, the current near the gap is larger than the current near 
6 180°. By narrowing the gap from @ = 1/10 to @ = 1/20, the overall level of current 
density was decreased, since reducing the size of the gap is essentially tantamount to 


decreasing the strength of the source. 
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The radiation patterns for the two gap sizes 0, = 0.1 and 4 0.05 (Figs. 7 and 8) 
show that the radius of the cylinder has a marked effect. For ka = 0.1, the pattern is 


uniform. With increasing values of ka the patterns become nonuniform. This is to be 
expected since when ka is small the current distribution has a uniform amplitude about 
the periphery of the cylinder, whereas when ka is increased, standing waves obtain. 

10. Acknowledgment. The authors thank Dr. C. T. Tai of Cruft Laboratory for 
his aid in the preparation of this paper. The computations were made by Mrs. R. Stokey. 
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THE EQUIVALENT LOADING METHOD AND 
THE EQUIVALENT BEAM METHOD’ 


BY 
PEI-PING CHEN 


Brown University 


Introduction. Mven under the usual engineering assumptions, the analysis of beams 
with non-uniform cross section subjected to distributed lateral loads, axial forces, and 
elastic support action is often difficult. It is the purpose of this paper to present a 
technique by which one can reduce more complicated beam problems to standard ones. 
This technique consists of two successive methods which we call the equivalent loading 
method and the equivalent beam method. 

In the equivalent loading method an arbitrary distribution of lateral loading on a 
portion of a beam is replaced by concentrated loads and couples so chosen that they 
are entirely equivalent to the given loading so far as the remainder of the beam is con- 
cerned. Thus a simple equivalent loading of concentrated forces and couples may be 
used to evaluate rigorously the stress and deflection that the given loading over a 
portion of a beam produces in the remainder of the beam. 

In an analogous fashion, the equivalent beam method is used to replace a portion 
of a beam by an equivalent, beam portion in such a way that the stresses and deflections 
in the remainder of the beam are not affected by the change. In this method the original 
portion may be of non-uniform cross section and may have axial forces and elastic 
support, but must be free from lateral load; the equivalent beam portion can always 
be chosen to consist of segments having uniform cross sections, free from axial forces 
and elastic restraint. 

It is assumed everywhere in the paper that the axial forces are such that elastic 


instability does not occur. 


PART I: THE EQUIVALENT LOADING METHOD 
Structurally balanced loading systems 


1. Fundamental fixed-ended beam problem. The equivalent loading method, which 
will be introduced in this paper to analyze complicated beam problems, depends on the 
analysis of the fundamental problem of the bending of a fixed-ended beam. The reason 
for this will be shown in Sec. 7. 

Let us consider a fixed-ended beam of non-uniform cross section subjected to arbitrary 
distribution of both lateral load and axial force (see Fig. 1); furthermore, for generality, 
we will assume the beam to be resting on an elastic foundation. The deflection function 


y(x) is governed by the differential equation, 


Lly) = —{ET(a)y"\" — [Pla)y’}' — kay = F(x) GR 
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together with the boundary conditions 
y(0) = y(l) = y'(0) = »y'() = O 12 


where £J(x) is the rigidity, always positive and finite, 
P(x) is the axial compression, 


k(v) is the foundation reaction intensity due to unit deflection (k(x) = 0, 


if the beam is elastically supported), 
and F(x) is the lateral load distribution (continuous or discontinuous). 


m— X {| F(x) 





—* Pe) 


Aq. “TET OY ‘B 


TXTTI7T 77 TIFT 


4 ee | 


Fia. 1. 











Here, as in all of the following discussion, the magnitude of the axial force is assumed 
to be such that buckling will not occur. The known functions, P(x) and EI(zx), are also 
assumed to be piecewise differentiable and twice differentiable, respectively, in the 
open interval (0, 2). 

In the derivation of Eq. (1.1), the following well-known formulas 


M(x) = —EI(z)y''(z), S(z) = M(x) — P(az)y'(z) 1.3 


are used for the bending moments and shear forces. The sign convention used is shown 
in Fig. 2. 





— 7 , e 
S G 


The solution of Eq. (1.1) together with the boundary conditions (1.2) may be ex- 
pressed in terms of Green’s function G(z, &): 


al 


y(xz) = | G(a, &)F(&) dé. 1.4 


Here G(x, &) is the Green’s function for the operator L[y] with the boundary condition 
1.2); stated in physical terms, it is the deflection at x due to unit lateral load at & with 
the effects of both axial force as well as the elastic foundation taken into account (see 


Fig. 3). 


2. Structurally balanced loading system. 
Definition. A loading system f(x) is said to be structurally balanced with respect to 
a fixed-ended beam AB, if and only if it produces no reaction forces and couples at 





= 





= EE 
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both ends. It follows from Kqs. (1.2) and (1.3) that these reaction forces and couples 
vanish if and only if 

y’ (0) 3 y”’ (l) - y’” (0) _ y’"(D = 0, (2.1) 


4* 

















— 
j / 
/ SIVSITITITTT TTT ITT LY, MY 
/ ba OF ~G (x,€) ) y, 
nine 
Fia. 3. 
whence, trom Kg. (1.4), 
G,(0, &)f(§) dé = 0, 
| GAL, &)f(€) dé = 0, 
(2.2) 
| G.2s(0, D§@ dé = 0, 
| G,.:(l, E)f(—) dé = 0. 
It is shown in Sec. 4 that the functions 
G..(0, &), Gi, £), G,,.(0, &), G.. (1. ) (2.3) 


must be linearly independent. 
The above equations (2.2) constitute necessary and sufficient conditions for struct- 
irally balanced loading systems. These conditions might be further simplified by the 


tid of a linear transformation (see Example I). 

EXAMPLE I: Yo find conditions which must be satisfied by all structurally bal- 
unced loading systems acting on a fixed-ended uniform beam under the following 
onditions: 


EI(x) = constant, P(x) = k(x) = 0. (2.4) 


SOLUTION: The Green’s function is found to be: 


wte 


a’(l — £)°(3tl — al + 22) /6EIT’, e< 


G(x, é) = 4 


fe 
A 
t 


lel — 2)*(Bal — €l + 2tx)/GEIL, 


¢ 








LS6 


from which we have 
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G,, (0, £) = &(l — £)’/EIT, 
G,, (l, ) = (1 — &)/EIP, 
2.6) 
G,..(0, &) = (l — £)°(2 — D/EII’, 
and tt, £) €(3l + 2¢)/EIT’. 


The above set of functions may 
independent functions, 
by a linear transformation 

£ = HIG,O, &) + 


Thus, according to Eq. (2.2), 


balanced loading system are the f 


be further simplified to the following set of linearly 


wtr 
ctr 
Jv¥ 


LIG,,.(l, =). ete. 


oj ~ 


- EIG,,Al, §) + 
o ‘ 


necessary and sufficient conditions for a structurally 


ollowing equations: 


3. Auxiliary functions of a beam. 
Definition. Four functions ¢,(£), g2(&), ¢3(€), ¢.(&) are said to form a set of auxiliary 


functions of a fixed-ended beam 
such that the relations 


| ¢(&)f 


« 
) 


if and only if they are linearly independent and are 


(&) dé = ¢e= 1.2.3.4 (3.1) 


hold for every loading system f(x) which is structurally balanced with respect to that 


beam AB. For example, the funct 


ions (2.3) form a set of auxiliary functions. 


In order to find the differential equation for auxiliary functions, we rewrite Eq. (1.1) 


in a more general form: 
Lily] = a(x)y"™ 


For any auxiliary function g(x) ar 


f(x), we have 


+ B(x 


2 9) 


yi" ae y(x)y”’ a a 6(x)y’ + exjy = F(x). (.. 


id any arbitrary structurally balanced loading system 


J o(é)falé)y" + BG)y’” + vy” + 5G)y’ + ey} dé 
é al 
= [Ay’” + By” + Cy’ + Dy} + | y(&) L*[p(€)| dé 
E an “0 


5 
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by integration by parts where A = a(é)g(t), B = — {la(é)e()}’ + B(e(8); ete., and 
L*|g] = falé)e}™ — [Ble}'" + lrvel” — [6e}’ + ee. This last operator 
L* is usually called the adjoint operator of L. Since y = y/ = y” = y’” = 0 at both ends 

see Eqs. (1.2) and (2.1)), we get 
[ y(é) L*[e()] dé = 0 (3.3) 

for which a sufficient condition is 
L*|o(&)| = 0. (3.4) 


lo show that Eq. (3.4) is also necessary for an auxiliary function ¢(é), we assume 
that L*{e(é)]} > O at some point ¢ By continuity, there exists a neighborhood of &, 


0 a: wll L 
*-Y) 


~ 
Fic. 4 





say (a, b), such that L*[y(é)| > O for € in (a, b). Let us construct a deflection function 


- 


(x) defined by (see Fig. (4)): 


y= 6 when x is not in (a, b) 
(3.9) 
>. when x is in (a, b) 
so that Eqs. (1.2) and (2.1) are satisfied, but not Eq. (3.3). Then there exists a struct- 
urally balanced loading system F(a) obtainable from Eq. (3.2) and Eq. (3.5) and an 


auxiliary function o(&) such that 


£)| dé + 0 


| g(t) I (é) dé | y(é)L* lo 


which contradicts hq. (3.3). Hence the necessity of Eq. (3.4) is proved. 

Since the operator L defined by Eq. (1.1) is self-adjoint, any set of linearly inde- 
pendent complementary functions” of that equation may be used as a set of auxiliary 
functions 

4. Necessary and sufficient conditions for structurally balanced loading systems in 
terms of auxiliary functions. Necessary and sufficient conditions for structurally 
balanced loading systems have already been obtained in terms of Green’s function (see 


Iq. (2.2)). Now we try to express them in terms of auxiliary functions. 


rom the definition of auxiliary functions, it is seen that 
| oe) fe) dé = 0, i= 1,2,3,4 (4.1) 


are necessary conditions for a structurally balanced loading system f(x). Here the 


gy 


o,(£) are the linearly independent complementary functions of Eq. (1.1). 


L. Inee, Ordinary differential equations, Dover Publications, New York, 1944, p. 402. 
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To prove the sufficiency of Eq. (4.1), that is, to show that a loading system f(z) 
is structurally balanced if it satisfies Eq. (4.1), we need to have the following additional 
structural stability condition: 


2,(0 gi(l ¢1(0) ¢i(L) 
g(0) ¢2(l) 03(0 5(l) 
Ale] = # 0 £.2 
¢3(0 v3( Ll) v3(0) ¢3(L) 
¢,(0) g,(l v4(0) g3(L) | 


For if any loading system f(x) satisfies the conditions (4.1), then by repeating the 
the argument of Sec. 3, we obtain 


1A ,(&)y’’’(E) + Be)y’"(8} = 6. 1=1,2,3,4 1.3 


due to the facts 
(i) that y(0) = y() = y'(0) = y'(l) = 0 for a fixed-ended beam, 
and (ii) that Z*[¢,(£)] = 0 obtained from the given condition (4.1). 
Since the determinant of the coefficients of y’’’(0), y’’’(2), y’’(0) and y’’(1) in Eq. (4.3), 


A,(0) A,(1) B,(0) B,()) 
A,(0) A,(I) B.(0) 3,(L) 
= {EI(O)EI()} Ale), 
A,(0) A,(l) B,(0) B,()) 
A,(0) A,(1) B,(0) BD 


does not vanish by the stability condition (4.2), we get Eq. (2.1). Thus f(x) is a struct- 
urally balanced loading system. 

It is interesting to point out from the above argument that the functions (2.3) are 
also linearly independent and they exist only when the stability condition (4.2) holds. 
EXAMPLE IT: To find the auxiliary functions for a fixed-ended uniform beam column 
under the following condition: 


EI(x) = constant, P(x) = constant > 0, k(x) = 0. 4.4) 
SOLUTION: The differential equation (1.1) becomes 
—Khiy* — Py" = F(z). 
Hence the auxiliary functions are 
l, g, sin &/J, cos &/j, (4.5 


respectively, where 7? = EI/P. And the stability condition (4.2) becomes 


ia e.. £3 ,. 5 l l .) 
Alg| = Y hoe 2j (sin 2j = 2j cos +) ~ 0 (4.6) 


which must hold for a stable fixed-ended beam column. 


; 
| 
i 
| 
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Structurally equivalent loading systems 


5. Structurally equivalent loading systems. 

Definition. Two loading systems, f(x) and g(x), are said to be structurally equivalent 
with respect to a fixed-ended beam AB if and only if they give the same reaction system 
at both ends (in other words, the difference, f(x) — g(x), forms a structurally balanced 
loading system with respect to the given beam). 

Generalizing this definition, two loading systems are said to be structurally equivalent 
with respect to a portion AB of a beam if they are structurally equivalent with respect 
to the fixed-ended beam obtained by clamping the points A and B. To be sure that 
two loading systems, f(x) and g(x), are structurally equivalent, we need* 


| g(t) f(é) dé = i] g(t) g(é) dé, += 1,2,3,4 (5.1) 


where the g; are auxiliary functions. It is understood that the stability condition (4.2) 
f J 


Thedttdeel . 


is satisfied. 
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6. Reversed reaction system. Fig. 5a shows a given loading system f(x), acting on 
a fixed-ended beam AB. The loading system produces a reaction system as shown in 
Fig. 5b. Then, by definition, the reversed reaction system (Fig. 5c) also constitutes a 
loading system structurally equivalent to f(x). Thus it is seen that there exists at least 
one loading system which is structurally equivalent to a given loading system—namely, 
the reversed reaction system. However, in some special cases, other useful types of 
structurally equivalent loading systems might be more easily obtained (see Fig. 6b). 

[t will now be shown how the auxiliary functions may be used to determine the 





*The following facts are noted: 

(1) If the loading g(r) consists of a concentrated load F at £, the right-hand side of Eq. (5.1) be- 
comes ¢i(é)F. 

(2) If the loading g(x) consists of a concentrated couple C at ~, the right-hand side of Eq. (5.1) 
becomes ¢’;(£)C; this can easily be shown by using two concentrated loads —C/e and +C/e 


t £and — + e, respectively, and then letting « — 0. 
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reversed reaction system. If ¢,() is a set of auxiliary functions and if Alg] 4 0, then 


Now if the ¢, are chosen so that 


oO v.(O) 00 g3(l v,(l) = v1(0) (), 
().2) 
(0) = g2(l) = ¢(0) = oi) = Aly] = 1. 
then the reversed reaction system may be readily found to be 
n o'(l) 
[. n ( 
n 
n yy (l elO a) 
L(O oll 6.3) 
n | ¢g (0) oh(l 
F, ¢ Fy 
n 0 | oh (| 
0 
N 0 Q | 


The above choice of g; is made possible by using the following linear transformation 


process on any given set of auxiliary functions y, , w , vw; and y, with Aly] # 0: 


LO v(l y.(E) LO v,( Vi(0 
¢ v.(0 y.(l W(é W.(0 y.(l W!(O 
v.,(0 W2(1 W(E L(O y,(l) W(0 6.4 
v(0 yi(l Vi(0 Wilé 
; / 
: v.(0 y( WO iE) 
os(é Aly 
¥ , . 
y,(0 W(1 W3(0 Wi(é) 
1.0 vi) VO) He 
Since Aly] + 0, all denominators in Eq. (6.4) are different from zero if a proper choices 


2 2 and 4is used. Thus the existence of a reversed reaction system is assured 


of indices 1, 2 


as is phvsically obvious for a stable beam). 
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EXAMPLE ITT: Given a tapered beam problem with the following conditions: 
E I(x) # constant, P(x) = k(x) = 0, (6.5) 


it is required (i) to find the auxiliary functions, 
ii) to find the reversed reaction system due to a unit load at & and 
(iii) to investigate some sufficient conditions for 
a) structurally balanced loading systems, and 
b) structurally equivalent loading systems 
in « special case in which the reciprocal of F1(&) is a polynomial in & 


¢: 


SOLUTION: (i) The Eq. (1.1) becomes 





EI (x)y")" = F(z). 


Hence the auxiliary funetions are 


y¥ - Po £ ov, = &(£), y, = T(&) (6.6) 
" | a at : : 
wher S(t) = j | | H(é) dé dé, 
on :  - 4 : _ 
T'(é) rE EH (é) dé dé, (6.7) 
and H(t) = 1/EI(é) 
with the value of A[¥] = —A/l # O where d is defined later in Eq. (6.9). 


ii) In order to find the reversed force system, we use the following auxiliary fune- 
tions (according to Eq. (6.4)): 


LS(é) 
1 | Co =¢€ / 3 — © oo s , 
»>— ¢ 
, : (6.8) 
l - ais 
o; ‘ (p — q)T(&) — (q — r)S¢&) 
A 
where the quantities p, g, r and \ detined by 
‘ —_ 
p = | H(é) dé, q = ] tH (é) dé 
(6.9) 


oT a 
r = | f° H(&) dé, A = i 


are all positive and finite. In this example the given loading system is a unit load at 
£, i 
(0) when x # &, 
f(z) = (6.10) 
l when x=. 
Hence we have n; = ¢,(£). Thus from Eq. (6.3), we obtain the reversed reaction system 


due to a unit load at. é: 
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F,=1-—F,, F, = [qS(&) — pT(&)]/r 
6.11) 
and OF =f — Fpl — Cp 9 Cp = De. 
iiia) If H(£), the reciprocal of HI(£), can be expressed by a polynomial of degree 
n in &, viz., 


H() = Dh, (E/D", (6.12) 


then S(é) and 7'(£) are exactly polynomials of degrees (n + 2) and (nm + 3) respectively. 
Hence under the conditions (6.5) and (6.12), we have the following sufficient conditions 
that a loading system f(x) be structurally balanced: 

| E™f(£) dé = O, m=0,1,2,--- ,n+ 3. (6.13) 


Ss 


For instance, if f(£) can be expressed by the sum of Legendre’s polynomials of orde" 
(n + 4) or higher, it constitutes a balanced loading system with respect to this fixed- 
ended tapered beam. 

(iiib) Under the same conditions (6.5) and (6.12) as mentioned above, we may also 
have the following sufficient conditions for the structural equivalence of the loading 


systems f(x) and g(2): 


| E"f(§) dé = | E"g(&) dé, m= 0,1,2, --- ,n+3. 6.14) 

t(§)=1 
[ITT] ud 
T¥rvvevevy 







12 = — f /> e— 2/5 = 2 




















Fic. 6a Fic. 6b. 


For example, if 1(&) is a polynomial of second degree as given by 


H(&) = ho + hi (El) + ho(E/D’, 6.15 


then the two loading systems shown in Figs. 6a and 6b are structurally equivalent 
because they satisfy Eq. (6.14) with m = 0, 1, 2, 3, 4, and 5 (compare with the reversed 
reaction system shown in Fig. 6c for ho = 1, hy = 0 and h, = 2 only). A detailed dis- 


cussion of various other types of structurally equivalent loading svstems* is beyond 


the scope of this paper. 


43 37 

ao 8 

— AL p2 -13 2 
{racl ON 200 


Fic. 6e. 














‘See, for example, P. P. Chen, Dyadic analysis of space rigid frameworks, J. Franklin Inst., Nov., 1944 
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7. The equivalent loading method. In ig. 7a is shown a loading system f(x) applied 
, portion AB of a beam PQ. Suppose that g(x) is a loading system structurally equiv- 


f(x) 


= (x) 


Fia. 7b. 





Fig. 7a. 


alent to f(x) with respect to this portion AB. By definition, if the loading systems 
—g(x) are applied simultaneously to the fixed-ended beam AB (Fig. 7b), they 
will produce no reaction forces and couples at the two ends A and B. Moreover, if 
and —g(x) are applied simultaneously to the portion AB of the original beam 


f(x 
also produce no reaction forces and couples at any support. 


PQ (Fig. 7c), they will 

























F(X) 
~~ : § (x) 
sail B P Mrrrt! Q 
ea ee 
A B 


-8 (x) 
Fia. 7e 
This follows from the facet that we may consider the portion AB removed from the 
beam PQ, then loaded with the loading systems f(x) and —g(x), and then replaced in 
PQ without having any condition of static equilibrium or geometrical con- 


the b 

tinuity violated. Now consider the beam problem shown in Fig. 7d, obtained by re- 
placing f(x) in Fig. 7a by g(x). Since Fig. 7a can be obtained by superposition of Figs. 
7e and 7d, and since the loadings of Fig. 7e produce no reaction, it follows that the 


loading produces the same reaction forces and couples at all supports as f(x) does. 


Thus we have 


THEOREM I: A loading system f(x), when applied to a portion AB of a beam, may 
' d by a loading system g(x), structurally equivalent with respect to the portion AB, 


/ ’ 
be replar 
without affecting the regions of the beam other than AB. 

Thus the equivalent loading may be used to determine shear, moment, slope and 


deflection outside of AB. Within AB, these quantities can be found by superposition 
of the effects of the loadings in Fig. 7e and Fig. 7d. 
PART II: THE EQUIVALENT BEAM METHOD 
Stiffness matrix and reaction matrix 
8. Fundamental simply-supported beam problem. Just as the equivalent loading 
method depends on the analysis of a fixed-ended beam, the equivalent beam method 
depends on the analysis of a simply-supported beam. The reason for this can be seen 


in Sec. 12. 
Let us consider a horizontal simply-supported beam AB (Fig. 8) of length | loaded 
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by two couples C, and C, . It is required to find the slopes and shears at both ends. 
In this case the deflection, y(x), is given by 


Liy| = 0, (8.1) 





together with the boundary conditions 





S(o) S(2) 


yit})) = y(L) = () (8.2) 
and EI@)y'(0) = —C,, Biber (b= Ce . (8.3) 
The general solution may again be expressed by 


y(x) = | G(a, &)F(&) dé (8.4) 


where the Green’s function G(x, £), the deflection at x due to a unit load at & acting on 
this simply-supported beam AB, is assumed to exist, and the load system /’(&) is con- 
sidered as the couples C, and C,; . Thus Eq. (8.4) gives: 


y(x) = G(x, 0)C, + Gi(a, DC» . 


i 


By differentiation, we have the required end slopes: 


’ y'(0) = G,(0, 0)C, + G.(0, DCz , 
8.5 
y'(l) = G,,(l, OC, + G,(l, Cw , 
or, written in matrix form, 
G,(0, Q) G,,(l, 0) 
w(O), 7D) = [Cx Col ; 8.6 
G,-(0, l) G,;(l, l) 
In the above equation, the square matrix 
G.(0,0)  G,;(I, 0) a 8B 
= say, 
G.(0, G,-(l, LD) B ¥ 


which may be called the flexibility matrix, is a symmetrical one due to the fact that the 
Green’s function is symmetrical with respect to x and &. Its existence, as will be dis- 
cussed later, depends on the fulfilment of the stability condition (9.6) for a simply- 
supported beam. However, its reciprocal matrix, 
a i a 8B 
= :, say, (8.7 


6 Y B ¥ 
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which may be called the stiffness matrix,’ is more useful in the moment distribution 
method. It is also a symmetrical matrix, but its existence depends on the fulfilment of 
stability condition (4.2) for a fixed-ended beam instead of (9.6) which is the stability 
condition for a simply-supported beam. 
EXAMPLE IV: For a simply-supported uniform beam under the condition (2.4), 
the Green’s function is 
fa(l — &)(2k — & — 2°)/6EII, x : 
G(x, ) = (8.8) 
t(l — x)(Qle — x — &°)/6ETI, Gx B. 


A 
rer 


Here the flexibility and stiffness matrices are, respectively, 


a 8 2 =] ] a 8 2 I OKI 
‘a 6EI’ = a (8.9) 
b ¥Y —| 2 = B Y I 2 
The end shears may be obtained from the static condition 
S(O) = SY) = —(C, + Cz) /l. (8.10) 


9. Stiffmess and reaction matrices in terms of auxiliary functions. The general 
solution of Eq. (8.1) may also be expressed in terms of auxiliary functions ¢,;(2): 


1 
y(z) = : a.0(2). (9.1) 
i 1 


From Eq. (9.1) and the boundary conditions (8.2) and (8.3), the coefficients a, may be 
determined. Again if the yg, functions given by Eq. (6.4) are used, we have, 


a, = a, = 0, 


y'(0) = a; , 


y’(l) = ase3(l) + a, , (9.2) 
y’’(0) = —C,/EIT(O) = aze5’(0) + ayei’(0), 
and y’’(l) = C,/ET(O) = ages’(l) + ayei’()), 


provided that the stability condition (4.2) holds. Eliminating a; and a, , we have 


a B 
gg Col = iD: 7! (9.3) 
8 Y 


where the elements in the stiffness matrix are given by 

a = [—¢f'(0) + of(lei’(0)|EI(0), 

8 = —gi/(O)ETO) = les’) — os(De (DIET, (9.4) 
and y = of (DED. 


6S. U. Benscoter, Matrix analysis of continuous beams, Trans. A. 8. C. KE. 112, 1109 (1947). 
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The non-vanishing of the determinant of the stiffness matrix, 
v3’(0) ¢3'(1) 


ay—p = — EI(NEI() ¥ 0, (9.5) 


g'(0)—h"(D) 


where the g, are given by Kq. (6.4), or, more generally, the non-vanishing of the fol- 
lowing determinant, 


¢(0) yg, (1) v’’(0) of'(1) 
¢2(0) ¢2(l) ¢2’(0) g2'(l) | 
Diy] = ~ 0, (9.6) 
| (0) 3(l) i’ (0) vi'(D) 
g4(0) o,(L) ¢’(0) 4’(l) | 


where the ¢; are linearly independent auxiliary functions, gives the stability condition 
for a simply-supported beam problem.” 
To find the end shears, we start from Eq. (1.3), 


S(z) = —{EI(ax)y’"(x)}’ — P(a)y’(x) 
9.7) 
=— Do a;{EI(x)g}’(x)}’ — P(x)y'(x) 
-3,4 
where the coefficients a; may be obtained from Eq. (9.2): 
] — v3(L) 
[az , a4] = [y’(0), y’(D) t (9.8) 
0 ] 
Putting x = 0 and I respectively in Kq. (9.7), we have the end shears 
€1) €\2 
S(0), SCD] = [y’(0), y’(D] 9.9) 
€ E20 
where 
e; = —P(0) — {EIle3’}/(0) + o3()) {EI¢%'}'(0), 
exo = —{EIg3’}/(D) + oD Elos’}'(0, 
é, = —jEIg4'}’0), 
éoo = —P(l) — {EIgs'}"(0. 
If the relation (9.6) holds, then the end shears can also be expressed by 
Pi Pie 
(S(O), SCD] = [C4 , Ca] 9.10) 
P21 P22 


‘For example, for a stable simply-supported beam column under the condition (4.4), the stability 
condition (9.6) gives Dig] = —(l/j*) sin (l/j) ¥ 0. It is interesting to point out that if a portion AB 
of a beam is such that D[g] = 0, then this portion acts as a hinged link. 
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where the square matrix 


Pi Pi a B €11 €\2 
= (9.11) 
P21 p B 7 €o) € 
may be called the reaction matriz. 
However if the beam is under the conditions 
kz) = 0, P(x) = constant, (9.12) 


then we also obtain Eq. (8.10) even if the relation (9.6) does not hold. Thus the re- 
action matrix has the following simple form: 


Pi Pi2 1 ] ] 
(9.13) 


P21 P22 ] ] 


with the distinctive properties that 
(i) it depends on the length / only 
and (ii) all elements p,; , pic , Po: aNd po. are equal. 
10. Simply. supported beams with two or more segments under the conditions 
k(x) 0 and P(x) = constant. In Fig. 9 is shown a simply-supported beam of two 


x, | A 
la é, “ 4 


eee 
cS 





1> 











segments AG and GB, of lengths J, and /, respectively. If the flexibility matrices are 


ay B, Qe Bo 
and , respectively, it can be shown that the elements a, 3 and y 
Dy 71 Be Y2 a 8 
of the resultant flexibility matrix of the whole beam are: 
8 Y 
a = a, — 2B, + ny + a), 
B= Bi + NoBo — NyNl(¥; + ae), (10.1) 
and Y = Y2 — 2n,B. + ny + a). . 


By repeated application of Eq. (10.1), the flexibility matrix of a beam of any finite 
number of segments under the conditions (9.12) can be calculated. 


Various types of equivalent beams 


11. Equivalent beams under the conditions (9.12). 
Definition. Two horizontal simply-supported beams AB and A’B’ (Figs. 10a and 10b) 
are said to be equivalent, if and only if they have the same stiffness and reaction matrices. 





198 PEI-PING CHEN (Vol. VII, No. 2 

Generalizing this definition, two beam portions are said to be equivalent if they are 
equivalent when considered as simply-supported beams. 

If these beam portions are not elastically supported, and if P(z) = constant, then 
they need only be of equal length in order to have the same reaction matrix (see Eq. 
9.13). In this paragraph we consider only this situation, unless otherwise specified. 

Let us investigate the following types of equivalent beams: (Type I) Let the beam 

a 8 
AB (Fig. 10a) have a given flexibility matrix , and let us find an equivalent 
BY 


beam A’B’ (Fig. 10b) composed of two uniform segments A’G’ and G’B’ with lengths 


4) 


f3 ¥! 






















hint a o> oH El, Eds 

i 

A’ -_—- B 
ea A 
-—2 hie £2 4 
Fig. 10a. Fic. 10b.  ~ 
El, EI, ET; 
a Bl 
| 4/3 | 4/3 ° 4/3 | 
Fic. 10c. 
I, and /,(=1 — 1,) and rigidities FJ, and EJ, respectively. Equating the flexibility 


matrix of beam A’B’ (as given by Eq. (10.1)) to that of beam AB, we get 


7 + 28 Ba 28 


lL, = ——— 1 l, = l 
ya a- ¥ 
(11.1) 
and EI, =s— oo EI, = = | ones 
' 3a + (a + 28)1,/1,’ ” 3y + (y + 28)1,/l. 
provided that 0 Ss 1, = l, 
that is, a ra — 28 S 7. (11.2) 
< = 
a B 
(Type Il) Assume that the beam AB has a given flexibility matrix , and find 
8 Y 


an equivalent beam A’’B” (Fig. 10c) composed of three parts of equal length but 
different rigidities EI, , EI, and EI, . We get 


EI, = 21/(1la + 2y + 148), 
EI, = —21/(7a + 7y + 408), (11.3) 
and EI, = 21/(1ly + 2a + 146). 
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Hence under the conditions k(a) = 0 and P(x) = constant, it is always possible to find 
an equivalent beam composed of only three uniform segments, if the use of zero or 
negative values of EI is permitted. 

However, if (9.12) is not satisfied, the elements p:; , pi2 , P2i aNd p22 in the reaction 
in general, not equal to each other. Hence, we cannot find an equivalent 
for the original portion AB (Fig. lla) unless some additional supports 
lib). A detailed discussion of equivalent beams for this case 


matrix are, 
portion A’B’ 
are provided (see Fig. 
is not given here. 


EXAMPLE V: To find the equivalent beams for a tapered beam (Fig. 12a) under 


the condition (6.5). 


he ee 
A a A'| = B' 
<—___—) —____» ———— 


Fia. 11b. 





























Fic. lla. 
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,” an 
Ane —8 | SRT Tar 
me E - ] l 
sane’ Jcmmmaaaes x 3/5) ae 2/54 
Fic. 12a Fig. 12b. 
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= Bf b/3~ 4/3 


Fig. 12c. 











SOLUTION: Use the auxiliary functions (6.8). 
From Eqs. (9.4), (8.7) and (9.12), we have 


a=r/xX, B = (q—n)/A, y = (p — 2q+71r)/r 


a=p— 2q+r, B=r-— q, y=r (11.4) 
and Pir = Pio = Pa = Po» = —I1/l. 
\s a numerical example, we assume that H(£) of the beam (Fig. 12a) is given by 
H(é) = 1 + 2(&/D)’. (11.5) 
From Eqs. (6.9) and (11.4), we have 
p = §1/3, q=l1, r=y = 111/15, a = 21/5, B= —4l/15. 


Then by Eqs. (11.1) and (11.3), the equivalent beams may be obtained as shown in 


Figs. 12b and 12c. 
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12. The equivalent beam method and an illustrative example. In Fig. 13a is shown 
a loaded beam with a portion AB which is free from lateral load. In Fig. 13b is shown 
the same loaded beam except that the portion AB is replaced by its equivalent one, 


A’B’. By reasoning similar to that used in Sec. 7, we have 


‘i | Y 
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Fic. 13 Fig. 13b 
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] 4s Fig. 14h. 


























THEOREM ITI. In any beam problem, a portion AB, if free from lateral load, may 
be replaced by an equivalent portion without affecting the regions of the beam other than AB 
As a final illustration of the method presented, we will show the application to the 


beam shown in Fig. I 4a. 


EXAMPLE VI. To simplify the beam problem shown in Fig. 14a by using both equiv- 
alent loading and equivalent beam methods (use | = 1 for convenience). 


SOLUTION: First, using the equivalent loading method to remove the distributed 
loads (see Figs. 6a and 6c), we get Fig. 14b. 

Secondly, applying the equivalent beam method to the portion AB in Fig. 14b to 
remove the varying cross sections (see Figs. 12a and 12b), we get Fig. 14c. Then the 
analysis of the beam in Fig. 14e will give precise results for the beam in Fig. 14a every- 
where other than the portion AB. Thus we may at least obtain the precise results at 
the ends A and B of the affected portion. These precise results might be then used in 
the determination of exact local value of shear, moment, slope and deflection within 
the portion AB. 

3y using these methods, we obtain the linearized M/EIJ diagrams which greatly 


facilitate the subsequent computations. 


A NEW THEORY OF PLASTIC FLOW* 


BY 
D. TRIFAN 


Brown University 


1. Introduction. This paper is concerned with a new mathematical theory of plastic 
flow for an isotropic, incompressible strain-hardening material exhibiting a gradual 
transition from the elastic to the plastic state. The discussion is restricted to states of 
stress and strain which can be reached by a single loading, followed by at most one 
partial or complete unloading. 

The new stress-strain relations are introduced in Sec. 2, and the criteria for loading 
and unloading are established. The boundary value problem for the strain rates corre- 
sponding to given surface velocities is formulated in Sec. 3 and uniqueness of solution 
is proved. Sections 4 and 5 are concerned with the minimum principle associated with 
this boundary value problem. Sections 6 and 7 bring a comparison between the new 
theory and the more familiar theory of plastic deformation.’ Finally, the conditions 
under which the two theories furnish identical results are investigated in See. 8. 

2. Stress-strain relation. The customary form of the stress-strain relation of the 
flow theory for strain-hardening materials [2, Eq. 27] gives the rate of change of the 
strain tensor as a function of the stress tensor, the rate of change of the stress tensor, 
and the stress invariants. The treatment of some problems of practical interest is greatly 
facilitated, however, by the use of stress-strain relations giving the rate of change of 
the stress tensor explicitly in terms of the strain tensor, the rate of change of the strain 
(ensor, and the strain invariants. Equations of this form cannot be obtained, as a rule, 
from the customary stress-strain relations of the flow theory. We thus are led to consider 


the following stress-strain relation for an incompressible material:” 


a? 2G.e%, — P(E)E*e,; . (2.1) 
Here the stress deviation s,;,; is defined as follows in terms of the stress tensor o;; : 
| 
Ci5 — = SisOun 
3 
As. usual, repeated indices indicate summation in accordance with the summation 


. 
i= 


convention of tensor calculus; 6;; is the Kronecker delta, ¢;; is the strain tensor, E = 
2e,,¢,;/3 the strain intensity, G) the shear modulus in the elastic range, and p = p(£) 
an arbitrary function depending on the material. Variables with asterisks indicate the 
rates of change of these variables with respect to a parameter ¢ which may be thought 
of as time.” Stress-strain relations of the type (2.1) may be called stress theories of plastic 


*Received Aug. 13, 1948. The results presented in this paper were obtained in the course of research 
conducted under contract sponsored jointly by the Office of Naval Research and the Bureau of Ships. 
1This terminology was introduced by A. A. Ilyushin [1]. (Numbers in square brackets refer to the 
bibliography at the end of the paper). 
2The author is indebted to W. Prager for the form of this relation which was introduced in [3]. 
3Since Eq. (2.1) is homogeneous in the rates of change, no viscosity effects are introduced, and 
the equation remains unchanged if ¢ is replaced by any monotonically increasing function ¢, i.e. by 
any other parameter which increases with increasing ¢. 
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flow in contrast to strain theories of plastic flow such as are treated in (2, Sec. 2]. Equation 
(2.1) is to apply for “loading’”’ only. The precise meaning of this term for states of 
combined stress will be given later. 

The arbitrary function p(£) in Eq. (2.1) can be determined for a particular material 
from the stress-strain diagram in simple tension. Denoting the tensile stress by o and 
the corresponding unit extension by e, the results of such a test for the type of material 
being considered here is given in Fig. 1. The stress-strain relation (2.1) reduces to the 
single equation 


Oo 





Fia. 1 
: Lde\. i+ 
P(E) = (4, — -—//E, (2.2) 
3 de 
while the strain intensity becomes E = ¢’. Assuming that o(¢) can be represented with 
sufficient accuracy by the power series 





c=aet+ae+a;e+-:--, (2.3) 

where a, = 3G, , and a; , a;, are constants, it follows from Eqs. (2.2) and (2.3) 
that p(£) is given by 

p(E) = b, + bE + b,k* + ---, (2.4) 

where b, = —(2n + 1) aon:1/3, n = 1, 2, 3, --- . The stress-strain relation (2.1) for 
loading then becomes 

st, = 2Goe%, — [b, + b.H + bk” + ---JE*e;; . (2.5) 


The stress-strain relation for unloading will be assumed to be given by the differ- 
entiated form of Hooke’s law: 

st, = 2G 7; (2.6) 

Under combined states of stress and strain, loading may occur in certain portions 

of a body and unloading in others, so that it becomes necessary to establish some sort 

of loading-unloading criterion. For that purpose we will assume that loading takes 

place wherever e¢,;;do;; > 0, and unloading wherever ¢;,do,;; < 0. If we multiply Kq. 

(2.1) by the strain tensor, we find that 

3 


2 


bo 
“J 


E*(G, — Ep(E)) > 0 


os 
Oy 55; 


for loading. On the other hand, all along the stress-strain curve of Fig. 1 the slope 
satisfies 

da a 
0<— < 3G,, 
dle 


vs pala 
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and therefore, by Eq. (2.2), 


Go 
E° 





0 < p(k) < (2.8) 
It follows from Eq. (2.7) and (2.8) that #* > 0 for loading. A similar procedure with 
Eq. (2.6) shows that the unloading criterion is E* < 0. 

The stress-strain relations for loading and unloading can be combined into one 


expression, namely 
s¥, = 2Goe%, — 3p(E)e,(E* + | E* |) (2.9) 


which reduces to Eq. (2.1) when #* > 0 and to (2.6) when E* < 0. 

3. Uniqueness of solution. With the aid of two lemmas which are to be established 
in this section the following uniqueness theorem will be proved. 

Theorem. For a given state of strain ¢;; throughout the body, and given velocities 
u* on the surface such that the surface integral of the normal velocity component 
vanishes, the rate of change of strain ¢*; is uniquely determined throughout the body. 

Lemma 1. Let €;; be some fixed state of strain throughout the body, and «*/ and 
e*/’ two arbitrary sets of strain rates which satisfy the condition of incompressibility. 
Furthermore, let o*/ and o*/’ be the corresponding stress rates obtained from the stress- 


strain relation (2.9). Then 
(et! — eF/’)(ot/ — ot!) > 0 


holds throughout the body independent of whether there is local loading or unloading. 
Proof. Setting d¥; = «*/ — «*/’ and noting that d¥,(o%/ — otf’) = d¥,(s¥/ — s*/’), 


we find from Eq. (2.9), that 
d*,(st! — st¥/’) = 2G,d%,d¥, — }p(E)dte,,{E*’ — E*’’ + | E*’ | — | EF” |}. 


Because 


{E*’ — E*”’ + | E*’| — | E*”’ |} | < | BY — Et’ | + || Et’ | — | Be’ | | 
<2)\fh" — g° 
and 
’ 4 
ye = i"! = — ¢ ay 
= ial 
we have 
* (o%! *IT 90) A* j* 4 7 * 2 
d*(s*/ — s*/’) > 2G,d*,d*, — 3 p(E)(d*e; ;) 
since p(E£) > 0 by Eq. (2.8). By the Schwarzian inequality 


: 3 a 
(dt.c,,)° < d?dte,.4,. = 9 didhE 


so that 
d*,(s*’ — st¥/’) > 2d*,d*,{G, — Ep(E)} > 0. 


ait 


This proves the lemma (see Eq. (2.8)). 
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Lemma 2. Let V be any “‘regular region” [4, p. 112] with surface S which is the sum 
, V, , and oF! and e¥/’ 


a 


of a finite number of “normal regions” [4, p. 85], V, , V2, 
functions satisfying the following conditions: 
1) in each normal region 


a) of! = of! , 

b) of! is continuous and differentiable, 
c) of, = do*/ /dzx, is continuous, 

d) o%/ ; = 0 (equation of equilibrium); 


2) the rates of change of the stresses transmitted across the surface of separation 
between two adjacent normal regions are the same whether evaluated from ¢*/ on one 
side or o*/ on the other; 

3) in each normal region 

a) e¥/’ is continuous, 


b) ¢¥/’ = 1/2(ut4 + u¥’/); and 
4) the velocities u*” are continuous throughout V + S. 
Then 
ote?!’ dv = | ot/ut’'n; d8, 
oo “SS 


where n; is the external normal to the surface. 


Proof. Using the conditions (1b), (3a), (3b), (1a), (1d), we find 


n . 


er n 
[ one d= [| ote d= [| (otfut”),, de. 
“VE s Fa 


1 


Use of conditions (1c), (2), (4), and the divergence theorem then gives us 


~ 


3 | (ot/u?") , dv = | ot/ut''n,; dS. 

It is now possible to prove the uniqueness theorem stated earlier. Let us assume that 
two solutions e¥; and e*/ exist each deriving from a system of velocities in accordance 
with 


eF = du, + uF; 


and 

ef? = d(u?’, + u?’,) 
where u* = u*’ on the boundary. Since o%; and o*/, the corresponding rates of change 
of stress obtained from Eq. (2.9), represent solutions, they must satisfy the equations 


of equilibrium, of; = 0 and o*/; = 0. It is assumed that no body forces exist. All 


the conditions of Lemma 2 being satisfied, we can conclude that 


| (et — ef)(o% — off) dv = | (ut — ut)(o%, — off)n, dS. 

v¥V “S§ 
But since on the surface u*¥ = u*’ the surface integral vanishes and hence the volume 
integral also. On the other hand, by Lemma 1 the volume integral is greater than zero, 
unless e¥; = e¥/. Thus, the solution e¥; must be unique. 


47 
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4. A minimum principle for loading throughout the body. We shall confine ourselves 
in this section to the case in which loading takes place throughout the entire body. We 
denote existing strains by e;; and the solution, i.e., the strain rates which actually occur, 


by «*,. These strain rates must satisfy the following: 
1) the equation of equilibrium o*;,; = 0, 
2) €& 0 
3) e? Liu* + u¥ ;) and 


1) €,,e%, > O, 


where u* equals the prescribed boundary values on the surface S. Let e¥/ be an arbitrary 
set of strain rates which satisfy the conditions: 


Ey eee Q 
2) e® Liu*’, + u*’;) and 
eee 0 
where u*’ c* on the surface. With these restrictions on the natural strain rates e¥, 


and the admissible strain rates e*/ the following minimum principle can be established. 
Minimum Principle. 


| eto dv < | etfot! dv 
Jy Jy 
when o*, and o?/ are obtained from e¥; and ¢*/, respectively, by means of the stress- 
strain relation (2.9). Equality holds if and only if e#; = €?/. 
In proving this principle we use the following lemma. 
Lemma 3. For the stress rates and strain rates introduced above 


, / 
tot; = 67; . 


Proof. From e#jo%/ = e%s*/ and e*/o%, = e¥/s*, , the lemma is obtained by applying 
Eq. (2.9) for E* > 0 and E*’ > 0. The equality also holds when E* < 0 and E*’ < 0. 
We note that Lemma 3 would not hold if, for instance, e¥; were to constitute loading 
and e*’ unloading. 

teturning to the minimum principle it follows from Lemma | that 


| e*.a* dv — | e*/a*. dv — | eXo%! dv + [ e*/a*! dv > 0, (4.1) 
J} Pe ig Jy éy 
where equality holds if and only if e#; = ¢¥/. According to Lemma 2, 
| etfot, dv = | ut’otn, aS; 
Jy Js 


and 


e 


| e*.0%, dv = [ utatn,; dS. 
7“V 


‘7 
“Ss 
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Since- u*’ = u* on S it follows that 


| et/a?, dv = | e*a*. dv. \4.2) 
J} Pa 

3y Lemma 3 and (4.2), we thus have 
| etot/ dv = | eto, dv. (4.3) 


i 


éy Jy 


Substitution of relations (4.2) and (4.3) into (4.1) completes the proof. 

5. Extension of the minimum principle. For the case where loading occurs in a 
sub-domain of the body and unloading occurs everywhere else,* some of the arguments 
presented in Sec. 4 are no longer valid because Lemma 3 cannot be applied throughout 
the body. Since the actual surface of separation between the two domains of loading 
and unloading for e*; is not known, the surface which separates the regions of loading 
and unloading for e*/ is in general distinct from the actual surface. Accordingly, the 
domain V is divided by these surfaces into what we assume to be four normal sub- 
domains, namely V,. , V_, , V._ , and V__ ; the first of the two subscripts refers to 
the actual state, and the second to the artificial state, the plus sign denoting loading 
and the negative sign unloading. 

Both (4.1) and (4.2) are still valid in this more general case, so that 


| eta?! dv < | et aF! dv. (5.1) 
“Vv Vv 
Since Lemma 3 applies over the domains V,, and V__ , we have 


J \ “Fy \ \ 
or 
| cor, @ = | e*/o* du + | Ldv+ | L di (5.2) 
} J “) J 
where 
L = Ft! — ef) ak, = €F,87/ fa? 


3v use of Eq. (2.9), we find that 


3 epee : : a a 
L= g P(E) E*’ E* | — £* | E* |}. (5.3) 
Using Eqs. (5.2), (5.3), (4.2) and the fact that E* < 0, E*’ > 0 in V_, and £* > 0, 
E*’ < O0in V._ . we write (5.1) in the form 
+ 3 f ”) Je fe! 3 / ”) Ee Je? */ &? 
| e*o*. dv — | D(E)E*E*’ dv + | P(E)E*E*’ dv < | et at; a. 
Ce | ae | + 7 “V4 “ys 


In the domains V_, and V,_ , we have pE*E*’ < 0 since p(E) > O (see (2.8)). 


‘It is assumed here that the domains of loading and unloading separate the body into two normal 
regions. This condition is not necessary and is introduced for the purpose of simplification only. The 


generalization to any finite number of normal regions is immediate 
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We conclude that the inequality 
| tot, dv < | etfot! de 
Jy Jv 
holds even in this more general case, whenever the artificial state is made to constitute 
loading in those regions where the actual state constitutes loading. In this case there 
is no V,_ domain.” 

6. Stress-strain relations for a theory of plastic deformation. The customary stress- 
strain relations [6, p. 75] for incompressible materials possessing the property of gradual 
transition from the elastic to the plastic state are 

$3; = 2G" €,; ’ (6.1) 
where G’ = G’(E) is the secant shear modulus. The secant shear modulus, being a function 
of the strain intensity F alone, can be evaluated for any material in the manner in 


which the quantity p(/) was obtained for the flow theory. 
For simple tension Eq. (6.1) reduces to the single expression 


o/e = 3G’. 
Using the same series expansion (2.3) as in the plastic flow theory, we obtain 
3G’ = 3G, +a,E +a;,E°+---. (6.2) 
Substituting this into Eq. (6.1) we find 
8; = 2Goe;; — [e, + oF + ¢,k" + ---]Ee,;, (6.3) 
where c, = —2a.,.,/3 andn = 1,2,3,---. 


7. Comparison between the two plastic theories. Consider a solid composed of a 
strain-hardening material which is deformed from a state of zero stress and strain by 
an equilibrium system of surface stresses such that loading occurs throughout the body. 
We assume that the relation between stress and strain in simple tension for the material, 
as given by the series expansion 


o=aetae +ase+-:--, (2.3) 


is known. The complete systems of equations for the two theories are then as follows. 
I. Plastic Flow Theory. 

a) Equations of equilibrium 

of, ; = 0, 
where body forces are neglected; 

b) condition of incompressibility 
ey, = 0; 

¢) stress-strain relation 

s® = 2 ye*, — [b, + DE + bE + ---JE*e,; , 


where b, = —(2n + 1)don+1/3; 


For a discussion of minimum principles based on strain theory of plastic flow see [5]. 
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strain rate-velocity relation 
fy = 3(ud,; + uf..), 
where the velocity vector u* is assumed small; and 
boundary conditions 
ir = [iF . 
where U* is the prescribed velocity of the surface. 


{ the solid, while I(e) 


Equations I(a)-I(d) must be satisfied at each interior point o 


must hold on the surface. 
II. Plastic Deformation Theory. 


ri) o, = 0 
D) € = 
S;; = 2G C, + Col esl” + -++|Ke;; , where ¢,, — 20 5 
a € = 5 (u; OU and 
r u-=U 


and 


+ 


Phe possibility of solving these non-linear partial differential equations for ¢ 
to 


in closed form is remote for most practically important problems. The procedure 


uM, 
in this paper assumes that the solutions may be represented by power series ex- 


be used 


pansions of the form 


0 5; = Zz, en ‘i 


2p+1 2prl 
é;; = > €é;; , 


i.1 


“= ui ‘i 


i BE Saale 
where the unknown coefficients of the odd powers of the conveniently chosen paramete! 
t, which may be thought of as time, depend on the space variables only. Substituted in 


the basic equations of the flow theory these expansions lead to the following relations 


+... =0; 


2G.{e;;) + Bet? + Sep tf + --:} (7.2 
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b(E,,t° + 2E,s¢° + ---) 


4 
whan 
Wl bo 


where E.. i, a 
d) 6; + Bei? 0 + Sef ti + e+) = Affusl) + Busse + Sujiti + ---} 
t fuss + Bust + Supt +--+}; 
and 
e) us? + Bure + Set eee = UY 4 BU + BU + ee. 


The unknown coefficients of the linear terms in ¢ of the expansions (7.1) are the 
solutions of the set of differential equations obtained for { = 0 in Eqs. (7.2), namely 


i) o 0 

D) eis 0 

( S QWGn€); (7.3) 
d) « Afus; + uj; ch; and 

i US 


These equations represent the basic equations for an incompressible elastic body, for 
which the existence [7, 8] and uniqueness [9, pp. 92, 93] are known. We assume that 
the linear equations (7.3) can then be integrated either in closed form or by some 
numerical method to yield of; and u{"’. Substituting Eqs. (7.3) into (7.2), factoring 


out ¢, and letting ¢ tend towards zero, we obtain 


a) o 0; 

b) « 0 

( 2G.€;;, — . be: : (7.4) 
9 

d) « = Afusth + u;*s} and 

e) u = ['° 


These equations, too, have the form of the basic equations for an incompressible elastic 


body. Indeed, they are equivalent to 
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a) "@:;- 3 =; Oe) rg Be 

b Phd = 0 

ec) 8 = 2G; (7.5) 
d) e = dus; + u;.3}; and 

e) ua; = [ 


where we introduce a fictitious body force expressed in terms of the solution of the 
first system of Eqs. (7.3). It is readily seen that all other unknown coefficients in the 
expansions (7.1) can be calculated from sets of equations similar in structure to Eqs. 
(7.8). 

Carrying out an analogous procedure in the case of the basic equations of the plastic 
deformation theory we find that the calculation of the unknown coefficients of the linear 
terms in ¢ of expansions (7.1) reduces, as for the flow theory, to the solving of Eqs. (7.3). 
The differential equations for o;;’, €;;’, ete. are those of Eqs. (7.4) with the exception 
that the factor 4b,/9 in Eq. (7.4¢) is replaced by the factor 2c, /3; however, from the 
definitions, 

4 2 4 


1 =3% = —9%- 


_ 


We thus find that the two theories give solutions which are identical as far as the terms 
linear and cubic in ¢ are concerned. 

The next higher order approximation will produce a difference between the solutions 
as given by the flow and deformation theories. For the flow theory, we have 


5 er 4 ; ; 16 , s 2 (1 —— 
si; = Ges; + T a3E,,¢€;;. + T: a3;E \3€;;, + a7 OsEne,; ; (7.6) 
oO o ad 


and for the deformation theory 


oO > 


5 y (5) 4 ’ 8 ’ 1 2 (1) 
8:; = 2G66:; + 9 a,E,,¢;, + 9 a;E y3€;; + a QF 51633 « ( 
¢ 24 


‘ 


~I 


“J 


Thus, from this stage on, the two theories will differ as a rule. 

8. Conditions for similarity. In order that the fifth order terms be the same for 
both theories, so that the solutions agree up to the terms of fifth order, we have the 
condition 

Bia? = Bye. (8.1) 


This equation is satisfied whenever 


3 


es; = f(z, y, Ze; » (8.2) 


where f(z, y, 2) is an arbitrary function. However, the strains ¢,;’ as given by (8.2) will 
satisfy the compatibility equations [9, p. 27] only if the function f(z, y, z) reduces to 
a constant. Accordingly, 

i Ban (8.3) 
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If we define ¢;;/E'’ as the normalized tensor of ¢;; , condition (8.3) states that the 
normalized tensors ¢:;’ and ¢:*’ must be identical. Substituting this value of «€{°’ into 
the stress-strain relations for s;;’ and adding a further restriction that 


ej; = Bej;’, (8.4) 


we find that now the two theories are identical for seventh order terms. Thus the first 


) 


four terms of the expansion for s;; , i.e., 
: ,(1) ,(3) 3 ,(5) 5 ,(7) 7 
8; = 8; t+ 38, € + .8;;€ + 8;; t 
are identical for both theories whenever 


eg, =e? (t+ AP + Be] + &7t'. (8.5) 


Higher order agreement may be obtained by similar restrictions. 

It is of interest to determine conditions under which the actual solutions as given 
by the theories of plastic flow and plastic deformation, and not as approximated by a 
finite number of terms of expansions (7.1), are identical. Differentiating the plastic 
deformation equations II with respect to time and comparing with the plastic flow 
equations I, we find that there is no difference between the theories if 


2ke*, — E*e,, = 0. (8.6) 


This‘condition is fulfilled if 


e, = 6; (2, y, Dg, 


where g(t) is a monotonically increasing function of ¢t. The two theories thus are seen 
to furnish identical predictions if the ratios between the various strain components at 
a generic point remain constant throughout the loading process. 
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—NOTES— 


A NOTE ON A SOLUTION TO POSSIO’S INTEGRAL EQUATION FOR 
AN OSCILLATING AIRFOIL IN SUBSONIC FLOW* 


By JOHN W. MILES (University of California at Los Angeles) 


Summary. Following a procedure analogous to that used by Schwarz for solving 
the integral equation for an oscillating airfoil in incompressible flow, a method of attack 
on the same problem for subsonic, compressible flow is outlined. The circulation and ker- 
nel to the circulation integral equations are expanded in series with the coefficients M™ 
and M*" In M, and a corresponding expansion for the Green’s function is thereby im- 
plicitly determined. The solution up to (but not including) terms of O(1/*, M‘ In M) is 
stated. A method of iteration for improving a given solution is also indicated. 

The circulation integral equation. The problem to be considered is that of a thin 
airfoil located near the x axis between x = +1 in an airstream having a subsonic velocity 
U in the positive x direction. The theory is linearized in the sense that the downwash 
w(x, t) (positive down) may be specified along the x axis between z = +1 and is small 
compared to U. The function which it is desired to obtain from w(z, t) is the circulation, 
y(z, t) which is so defined that the pressure jump, II(2z, ¢) at any point on the airfoil is 
given by the Kutta-Joukowsky law 

I(x, t) = pUy(a, t). (1) 
For a harmonic time variation of the form exp (iwé), Kiissner’ has shown that y(z, ¢) and 
w(z, t) are connected by the integral equation 


a+1 


J : : — . ; 
w(x, t) = a | v(é, t) exp [ik(é — x)]R[M, wk(x — §)\(x — 8)* dé, (2) 
a7 J} 
where M is Mach’s modulus, u is the Lorentz contraction factor (1 — M’)~'”’, k is the 


frequency parameter wbh/U (b, the semi-chord, is the unit of length used throughout the 


anslysis), and R(M, y) is given by 


R(M, y) = = My | exp Gu)H?(M |u|) |u| du, (3) 


H(z) being Hankel’s cylindrical function of the second kind.* The notation has been 
converted to that generally used in this country; Kiissner uses 6 = M, v = @, w = tk. 
The incompressible problem. Soehngen® has shown that the only solution to 


o2)= = | s@@-era (4) 


*Received March 29, 1948. 
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ea 


g(a) = (2m)7 | (x — &)~'f(®) dé 


v-—a 


und deren Anwendung in der Tragfligeltheorie, Math. Z. 45 (1939). 
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for which f(1) is finite (in the present problem, the Kutta condition), is given by 





al as 9/1— xy 2 pl aa": _ e 
fee) = — 2(4=4)" f° o@(t+ i) @ — a" ae. (5) 
Using this result, Schwarz‘ has given the incompressible (4M = 0) solution to (1) and 


the Kutta condition 


y(1, t) = 0. (6) 
His solution is 
ee a 
y(cos 6, t) = rs | G.(k, 6, v)w(k, ¢, t) dy, (7) 
x = cos 8, & = cosg, (8) 


1 — cos (6 + ¢) 


2(cos ¢ — cos 6)" sin @ 
1 — ce (Oo (cos ¢ cos 8)” sin 


G(k, 6, ¢) = tk In | 


+ [(1 — cosy) + T(k)(1 + cos ¢)] tan (2), 


—tHy'(k) + Hi’ (k) (10) 


i Media oe a 
iH,” (k) + Hy” (k) 


where H;” and H<*’ are Hankel functions in Watson’s notation.’ 
Compressible case. In seeking a solution to (2) for M # 0, it is expedient to introduce 
the modified functions 


w*(x) = pw exp (—ixx)w(2), (11) 
y*(x) = exp (—7txx)y(2), (12) 
ar _M nae 2 1% ¢ 
c= (. = )k = M*k*, (13) 

k* = pk (14) 


and rewrite the integral Eq. (2) as 


ol 


w*(z, t) = | K(M, k(x — dIy*t, 0 ab, (15) 
K(M, y) = 5 y~ exp (—éy)R(M, 9). (16) 


‘L. Schwarz, Berechnung der Druckverteilung einer harmonisch sich verformenden Tragflache in ebener 
Strémung, Luftfahrtforschung 17, 379-386 (1940). 
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Expanding the Hankel function in (3), the kernel in (16) may be written 
K(M, y) = > exp (—7zy) | : exp Guy, - u~® 
= (—)"M?"*7y?” | ( | ) l — (1) | 
+ —|In|M +C- —— - -~)}+ =I? du 
p> 2°"m\(m + 1)! 2 2(m + 1) 2. n 2 


where C is Euler’s constant. It will generally be convenient to integrate each of the terms 
retained in the integrand repeatedly by parts until the only integral remaining in the 


(17) 


kernel is [¢.. wu’ exp (tu) du. 
A solution (15) can be carried out by writing 
K(M, y) = >> M*"K,,(y) + In(M) > MK, , (18) 
y*(x) = D> M*"y2,(x) + In (M) >> M?"3,(x) (19) 
n=O n=l 


and proceeding, in a manner entirely analogous to that utilized by Schwarz, towards a 


result of the form 


7(6) = . | exp [7x(cos @ — cos g)|G(k*, 6, g)w(k, ¢, t) de, (20) 
G(k*, 0, ¢) = >> M”’G,,(k*, 0, ¢) + In(M) D> M*Gi,(k*, 8, ¢). (21) 


The algebraic manipulations required to solve for the G,, and G3, are rather lengthy, 

and it is probably not feasible to include more than a few terms. However, the solution 
. ° ° mn “¢ (N) ° 

can be improved by iteration. Thus, if y*'"’ denotes the solution to 


el 


wt = | K™y*™ dé, (22) 
J} 
the remainder [y* — y*'*’] is given by the solution to 
al al 
wr = | [K“% — K}y*™ dt = | K[y* — y*‘"’] dé, (23) 
e 1 s 1 


. ° » . + - . (N 
where the left side to (23) may be regarded as a new velocity distribution w*"’. 
The author has carried out the solution for N = 1*. In this case G, is given by (9) 
noting that k must be replaced by k*), while G, and G} are given by 


( m\2 
sae re oe fe) oe .45) 
G.(k*, 6, ¢) = k* tan (S\a + cos | i(a + 1) sl 5) 


/ 


\Vl+T . a _, 
+ ar i)( +1) + i _ (itt) tan (2) sin’ ¢ 


— k** sin | o( + a + cos yg) + in*(a - i) sin’ e| 


*Details of the solution are available from the author. 
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k* /i+T 
aa E \ 9 ) sin 6 cos 6 
ee k* 1 — cos (@+o 
k*° sin o(cos 6 — cos ae — (cos 6 — cos ¢) In = r |, (24) 
a P. : | — cos (6 — o). 
t 
= \r had * k 
‘ pet @ 1 6\) (1 1 ae 
Gi(k*, 0, = ( s ) tan (2) ( + a + cosy) — k* sin’ ¢ iy 
- 7 . —— ~ 0 
em? /1+7'\ . 1+T' : | I 
7 | 5} } sin Hl ( - Ja + cosg) + 7k* sin g |, (25) c 
c 
ou l ir | k 
a= ; + :) + 5 In (‘), (26) t 


The function [(1 + 7'(k*))/2] is the C(k*) used by Theodorsen.” 

Other solutions. Other solutions of (2) have been given by Possio’’’, Schade’, ’ 
Dietze” and Eichler’’. The only extensive numerical results available are those of Possio 
for the lift and moment on a flat plate airfoil due to plunging and pitching motions. The 
collocation technique used by Possio is probably not practica! for control surface motion. i 

Due to the large expenditure of labor which will be required to calculate subsonic 
flutter derivatives over (usefully) complete ranges of M and k, it appears that a prelimi- 4 
nary study of all available methods should be conducted, probably under the auspices of 
one of the government agencies. While it appears to the author that the work of Schade" t 
is probably best suited to a large scale computing program, some calculations should prob- 
ably be made using each of the methods mentioned. 


Further discussion of the German work has been given by Biot, et al.”’ 


T. Theodorsen, General theory of aerodynamic instability and the mechanism of flutter, NACA T.R 
196 (1935). 

°C. Possio, L’azione aerodinamica sul profilo oscillante in un fluido compressibile di velocita iposonora, { 
L’Aerotechnica 18, 441-458 (1938). 


™R. A. Frazer, Possio’s derivative theory for an infinite aerofoil moving at subsonic speeds, A.R.C. 4932, ; 
0.205, 1941. ;' 
‘Schade, The numerical solution of Possio’s integral equation for an oscillating aerofoil in a two- 
dimensional subsonic stream, Aero-Versuchsanstalt, Géttingen; Reports B 44/J/27, 31.8.44; B44/J/43 1 
27.11.44; B 44/J/44, 1.12.44, Trans. by S. W. Skan, Aero Dept. NPL, 1946 (Brit.). 1 


'Dietze, Die Luftkraft des harmonisch schwingenden Fliigels im kompressiblen Medium bei Unter- 
schaligeschwindigkeit (Ebenes Problem), Deutsche Versuchsanstalt fiir Luftfahrt, 1943 (AAF Tr. F-Ts- 
506-RE, Wright Field, 1946). 

0M. Eichler, Auflésung der Integralgleichung von Possio, Forschungs-Bericht Nr 1681, Junkers 
Flugzeug und Motorenwerke A-G, Dessau, 1942 (photostat available through N.A.C.A.). 

uM. A. Biot, S. N. Karp, 8. S. Shu, H. Weil, Aerodynamics of the oscillating airfoil in compressible 
flow, AAF Tech. Report F-TR-1167-ND. 
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THE NORMAL EQUATIONS OF THE METHOD OF LEAST 
SQUARES AND THEIR SOLUTION* 


3y M. HERZBERGER (Eastman Kodak Company) 


The literature of the numerical solution of linear equations is extensive.** The fact 
that war work forced many mathematicians into contact with numerical work has 
brought with it a better insight into the mathematical problems involved. In 1947, von 
Neumann and Goldstine’ discussed an analysis of the accuracy of the best-known methods 
of solving linear equations and an analysis of the steps involved. The present paper 
provides the applied mathematician with a method by which he can improve the ac- 
curacy of his solution, step by step. At the same time, an attempt is made to guide the 
computer by transforming the problem into one of analytic geometry. 

The problem of solving n linear equations with n unknowns is equivalent to finding 
the components of a vector y with respect to n independent vectors a; : 


wer 


a; = y. (1) 
Multiplication by a; gives the equations 
Eia;-a; = y°a;, (2) 

having a symmetric determinant. 

The normal equations of the method of least squares can be written in the same form. 
The geometrical problem is as follows: 

Given in n > k-dimensional space the k vectors a, , --- , a and a vector y, we want 
to determine a vector d = £,a; such that 


(y — d)* = min., (3) 


i.e., we want the best representation of y as a linear manifold of a; . 
The solution is obviously the vector d, which is the projection of y into the manifold 
a; so that 
(y — d)-a; = 0, (4) 


a result which also can be obtained by differentiating Eq. (3) with respect to a; . In- 


serting into Eq. (4) the value of d, we obtain 
€,a,-a; = ya; , (5) 


the normal equations of the problem of least squares. These equations are identical 
with Eqs. (2), the only difference being that the a; are vectors in a space of higher 
dimensions. 

Equations (5) are the normal equations of the problem of least squares. Their solu- 
tion, in general, presents a difficulty if, and only if, the determinant of a; is small, that 





means, if the vectors a; are ‘“‘nearly linearly dependent.”’ 
In an earlier article’ the replacement of the vectors a; by unit vectors was suggested: 


e, saat a, (a;)*”- f im y/(y’)'”. (6) 


*Received July 26, 1948. Communication No. 1207 from the Kodak Research Laboratories. 
**Comprehensive summaries have been given by Hotelling*, Dwyer,’ and Bodewig.‘ 
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bo 
—_ 
@ 


In this case, 


(7) 
e;-e; = cos (@; , e;), 


i.e., the matrix of Eqs. (5) has all its diagonal elements equal to unity and the other 
elements smaller than unity. 

We now re-order e, , --- , &, with the aim of choosing e, and e, such that [e,e,]’ has 
the largest possible value, choosing e; such that [e,e2e,]” has the largest possible value, 
and so on.* 

For practical reasons, this criterion may be replaced by the following, which, although 
not rigorous, is easier to apply. We choose e, , ez such that (e,e,)” has the smallest value, 
then add e; such that (e2e;)” + (e,e;)” has the smallest value, and so on. These values 
form, on account of Hadamard’s famous theorem, a majorant of the values [e,e, --- e;]’. 
The re-ordering is done to investigate whether a projection into a space of a smaller 
number of variables gives a sufficiently close approximation. We must keep in mind 
that in practical problems the data are only known to a certain number of digits. There- 
fore, it is only necessary to determine a solution so as to reproduce the accuracy of the 
data. 

To solve the normal equations, two procedures seem possible: one is to change the 
right-hand side, i.e., to replace the vectors y successively by smaller vectors; the other 
is to replace the vectors a; on the left-hand side by a series of other vectors which permit 
the solution of the normal equations. We shall make suggestions with respect to both 
methods. The determinant of the normal equations is the discriminant of the quadratic 
form, 


(> a.é.) = > a,a,fé; ? (8) 


and as such is positive and definite. Since it is derived from a finite number of digits in 
machine computation, however, the number of significant figures in the calculated value 
of the determinant may be small. In either case, the homogeneous equations (right- 
hand side equal to zero) have only the solution: £, = --- = & = 0. 

Let us assume that we have given a vector b which might be an approximation to y. 
Such a vector b permits us to obtain a new and smaller value for the right-hand side 


of Eqs. (5). If we introduce 
b = Ba, (9) 
into Eqs. (5), we obtain 
f,a,a; = (y — b)a; = ya; — 6,a,a; . (10) 
The sum of the squares for the new vector is given by 


(y — B.a:)” = y’ — Biya, + B,B,a.a; . (11) 


*The symbol [e¢2 ... ¢,]? is an abbreviation for the determinant of the products e,e; for 1,7 = 1, 
- , p, taken from Grassmann’s symbolic notation. It is the square of the volume of the parallelepiped 


formed by the vectors ¢, ... , @&. 
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The reader may notice that the computation of the coefficients in (10) and (11) 
presupposes only a knowledge of the 6, the a,a; , and the ba; , i.e., the data of the k + 1 
by k matrix. It is not necessary to go back to the original vectors in n-dimensional space. 

If a number of approximation vectors b, , b. , --- , b,, are given, then it is possible 
to compute the best multiple of them to be taken from the right-hand side of the equation. 
The best approximation is given by 


B = mb; (12) 
with 
njb;b; = byy. 
The reader may notice again that the expressions in Eq. (12) can be computed from the 
coefficients of b; with respect to a; without having recourse to n-dimensional space. 

Of much greater importance for the solution of the normal equation is the trans- 
formation of the left-hand side. 

There are two different transformations of the left-hand side of the equation which 
have played a role in the literature of the problem. We shall analyze these transforma- 
tions in the remainder of the paper. The first method will be called “the improved 
Gauss-Doolittle’’® method. The second one is the square-root method found independ- 
ently by Schur,® Banachiewicz,’ and Dwyer.* We shall study both methods theoretically 
for the general case in which the vectors are not reduced, but we recommend the second 
method particularly, for the reduction to unit vectors. 

Strangely enough, the best way to discuss these methods seems not to be found in 
the literature.* Both methods can be described as using orthogonal vectors. Let us 
discuss first the Gauss-Doolittle method. We find a system of orthogonal vectors 0, , 

- , o such that 0, = a, , that o, lies in the plane of a, and a, , that 0; is a linear 
combination of a, , a , ag and such that in all cases 


[a,a,a;]" = [0,0,0,]°. (13) 
This leads to 
a= 0, 
a = 6.0; + O2 , 
(14) 


A; = 5:30; + 52302 + 0; , 


The matrix A of the 6,; is a triangular matrix which we can calculate as follows. 
From Eqs. (14) we find that 


2 2 2 
a, = 01, aa; = 4,;0; , 

2 2 2 2 2 
a. = 6:20: + 2, A.A; = $126,;0; + 42;52 , (15) 


2 2 
a; = 5,30; + 52302 + 053 ’ a,a; = 5435; ;0; + 52352;02 + 53;03 . 





*As far as we know, Kolmogoroff is the only other author who has used vector notation for the 
least-squares method. His paper, published shortly before our first communication,’ reached us only 


recently. 
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If we abbreviate the diagonal matriz, 


1o; 0 0 0 0 
10 0 0 0 Oo 
| |= O (16) 
, 
‘ 
iO 0 0 0 O; 
these equations can be written in the matrix form 
A = A’OA. (17) 


Since a, and a; are given, the computation of 6;; can now be done in the following 


way. Equations (15) give the relations 


0; = a, 5,; = (a,a;)/o% , 

- 
O; = & — 51201 , 52; = (@:8; — 51261;0;)/02 , (18) 
0; = @; — 53301 — 52302 , 53; = (a3€; — 5136,;0; — 52352;02)/03 . 


Having obtained o; and the triangular matrix, 6;, , we can solve the least-squares 


equations in two simple steps. The equations Af; = n; are equivalent to A’OAE; = 1. 
We write 
Ag; = f, (19) 
and 
OA’S, =. (20) 


These equations are easily solved because of the triangularity of the matrix A. 


Written in full, Eqs. (20) are 


0; 613f1 + O2de3f2 + O3f3 = m3 , 
or, solved 

[1 = m/0; 
f2 = (nz — 01642f1)/0 (22) 
ite a ( 2e rm 2¢ ~ a 
63 = (M3 — 0101301 — 0202362)/03 .« 

The computation of &, is then obtained by the equations 

&, — oe ] 
(23) 


€,-1 = Cex = Ox, 4+18k . 
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It may be noted that the ¢; are the coefficients of the vector x = £,a; with respect 

to o, , for 
i= £,a; = £;6,;0, = 6,0; ° (24) 


The discussion of the square-root method is similar, except that, instead of choosing 
orthogonal vectors which fulfill Eq. (13), we choose orthogonal unit vectors o, in the 
direction of 0; . This leads again to a triangular matrix of coefficients: 


a, = 7119; , 
A. = 1201 + 2202 , (25) 
@; = ¥1301 + Y2302 + 3303 , 


with the distinction that 
A=MYT, (26) 
since 
QB; = Visi - (27) 


Hence, by this procedure the original matrix becomes the product of two triangular 
matrices. 
The computation of y,; proceeds similarly to that of 6;; . Equations (25) give 


V1 = ¥ YuVis = 4,8; , 
(28) 
Vi2 + Yo = a Y12V1; + Yo2x¥2; = A; . 
Thus, we find that 
: oa, Niu = 44; , 
: Y22 = (a; = Ya): Y2iY22 = 228; — Vii7Viz ; (29) 
Ys3 = (a3 — Vi3 _ Y23)' 4 Y3i¥33 = @38; — Yoi¥ee — VYiYu - 


The computation of square roots on a calculating machine can be reduced to simple 
division if an approximate value is known, because, for any value z and a small increase 


dx. we have 
a/(x + dx) = x — dz. (30) 


Having found the value of y,;; , we can proceed with the solution of Eqs. (2), again 
in two steps: 


I’Té; = Aé; = oe (31) 


which splits into 


and (32) 
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Again we find that because of Eqs. (25) and (32), 


x = £,a; = &,0;, (33) 


so that the ¢ are components of x with respect to 0; . 

The advantage of the second method in connection with making the vectors a and 
y unit vectors is that all the numbers calculated, with the exception of the final &; , 
are projections of vectors of unit length into a Cartesian coordinate system; they, 
therefore, are all smaller than one. The accuracy of all the operations can hence be pre- 
served to the last significant figure, since division does not lead to a loss of significant 
figures. 

A check consists in forming 

CT = 4 (34) 
and making sure that the elements of A are reproduced to the last decimal. 

These remarks hold up to and including the computation of the ¢. If the values of 
£; , when inserted into Eqs. (5), leave small residuals, a correction can be obtained by 
repeating the solution of Eqs. (5) with the residuals at the right-hand side. A small 
number of digits will be sufficient to carry through the computation in this case. 

The author suggests the following procedure in solving a system of linear equations 
having a small determinant and a large number of variables. The first step is to change 
the vectors a; and b; into unit vectors and to rearrange the vectors in the manner de- 
scribed previously. We then solve part of the equations according to the square-root 
method, until y,, has sufficient significant digits, using the number of decimals which 
the calculating machine permits. This gives a solution for p < k unknowns and re- 
duces the first p numbers on the right-hand side to zero or nearly zero. These values 
are put into the remaining k — p equations, part of which are solved, and the same 
procedure is followed until values of all the unknowns are determined. If the sum of 
the squares of the residuals is small enough, we stop. If not, we repeat the procedure. 

Going through it for a second time, we look for the best multiple of the last solution 
to diminish further the sum of the squares. The procedure is repeated until the sum of 
the squares becomes compatible with the error of the original data. 


Acknowledgment. I should like to express here my gratitude to Mr. R. Morris, who 
worked with me on many theoretical details of the paper, and to Dr. F. Kottler, of 
these Laboratories, for his helpful criticisms. 
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ON REISSNER’S THEORY OF BENDING OF ELASTIC PLATES* 
By A. E. GREEN (Durham, England) 


1. Introduction. The classical theory of bending of elastic plates has recently been 
extended and improved by Reissner.’’*’* His theory takes into account the transverse- 
shear deformations of the plate and the equations of the theory are obtained by an 
application of Castigliano’s theorem of minimum energy. The object of the present 
note is to show that Reissner’s equations can be obtained directly from the stress equa- 
tions of equilibrium and the stress-strain relations. Moreover, by consistent use of 
complex variable notation, the form of the results is simplified. The equations are first 
obtained for an isotropic material and are then extended to an aeolotropic material 
which is transversely isotropic in planes parallel to the faces of the plates. 

2. Fundamental equations for isotropic plates. Consider Cartesian coordinates 
xz, y, Z and let z = x + ty denote the complex variable with Z = x — iy the complex 
conjugate of z. Stresses connected with the coordinate Z are denoted by Tes , Tys » Os 
since there is no need to confuse the z in this notation with the complex variable. At- 
tention is directed to stresses in plates bounded by the planes Z = +h. 

When body forces are absent, Stevenson‘ has shown that the stress equations of 
equilibrium can be expressed in the form 





OP 00 ov 

aa 2” 0, (1a) 
av , WW , do, _ 

a*s’*a’” (1b) 


where a bar placed over a quantity denotes the complex conjugate of that quantity 
and where 
6=e,+4,, © = o, — o, + 2r.,, V = 7, + it, - (2) 
If u, v, w denote the Cartesian components of displacement and if D = u + iv, the 
complex form of the stress-strain relations is 


(1 — 2n)0 = 2 UP + F) +2n oul (3a) 


*Received Aug. 13, 1948. 

1—, Reissner, J. Math. Phys. 23, 184-191 (1944). 

2E. Reissner, J. Appl. Mech. 12, A68-A77 (1945). 
3—. Reissner, Q. Appl. Math. 5, 55-68 (1947). 

4A. C. Stevenson, Phil. Mag. (7) 33, 639-661 (1942). 
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“ 0,.2-c%UP — ow a 

2n)o 2u) 7 Vi(F + F) + (1 — n) az? (3b) 
92h 

& = 16u —, (3¢e) 
oz 

re) OF \ 
Y=72 (2 2 t wy), 3d 
- Oz az * , 


where 


OF , om O° 0 
4— V;= 4 — = 3 + 3e) 


az” Oz Oz Ox ay” 


D= 


and where Poisson’s ratio 7 and Young’s modulus F are related to Lame’s constants by 


r K ,; 
= i = } 
’ 2(\ + pw)’ . 21 + 7) 
3. Change of axes. Stevenson has pointed out that the stress combinations 0, ®, V 
are particularly suitable to the problem of transformation of stress and his main results 


are recorded here. Thus, if 


oO = 6, + ¢,.; @’ = o, — 6, + 2Tn 5 . se oe 5a) 

where 
Yahi a =e. 5b) 

then 
6’ = 6, ® = eS, WwW =e, D’ = e"**D. 5e) 


4. Formulation of the problem. I[t is assumed that the faces of the plate are free from 
applied shear stresses so that 


v¥=0 (Z = +h), 6a 
while the normal traction on the faces is such that 
o, = + 3p (Z = +h), 6b 
where p is a given function of x and y. 


Stress-resultants and stress-couples are defined by 


ah ah ph 


M.=| Zo.dZ, M,=| Zo,dZ, H=]| Zr,,dZ, 


~] 
os 


or, in complex form,” 


*These definitions are different from those given by Stevenson [4]. 
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A=M,+M,=[| Zaz, 


h 
r= M,—M,+ 2H = [ Z® dZ, (7b) 
Jah 
Y, = V.+i1V,= | WV dZ. 
J—h 


By using the boundary conditions (6a) and (6b) and Eqs. (la) and (1b) it is found 
that 


a 
Oz Oz 
(8) 
) 0 4) 0 
oat 4, Set yeh 


dz Oz 
lhe boundary conditions (6b) may be satisfied by taking 


s,=-€ {3(2) ws (2) (9a) 


tU\h h 
and then, from (1b), 
3 f (2) 
v= 41 - VY, , (9b) 
Ah | wy), 
which gives zero shear stresses at the faces of the plate Z = +h. 


Weighted displacements are now defined by * 


* : ; 
Dt = 4—, Ft = — | ZF dZ, 
(10) 


\ 


Pe Se is f _ FT 
wt = 7h | at (5) pu dZ. 


Relations between these weighted displacements and the stress-resultants and stress- 
couples can be found by multiplying Eqs. (3a) to (3d) by suitable functions of Z and 
integrating with respect to Z through the thickness of the plate. Before carrying out 
this process it is necessary to eliminate déw/@Z from Eqs. (3a) and (3b) in order to avoid 
introducing further unknown quantities. Thus Eqs. (3a) and (3b) are replaced by 


(1 — n)O — 2no, = Qu(l + n)Vi(F + FP), (11) 


and the remaining equation which involves dw/d@Z is not used. Hence 


1Ouh 0 . ; 
VY, = —— — (w* + 2F*), (12a) 
3 Oz 
*The choice of a weighting factor Z for the displacement D is natural. Equation (3d) then shows 
that w must be weighted with a quadratic factor in Z which vanishes at Z = +h, in order not to intro- 
duce further unknown quantities involving values of F at the faces of the plate. 








226 NOTES [Vol. VII, No. 2 


_ 82uh® aF* 
~ @o wf’ 





—3 h'np fe Ault + »)V2(F* + F*). 


(1 — mA 


Equations (8) are now satisfied if 


Shiu 


‘orm ve" = 


ee w* + 2F*), 


Vw + F* + Pe) = 22. 
5uh 


5. Solution of the differential equations. If 
F* = $+ ty, 
then, from (13b), 


2 


w+ mes fey @) + are} + ar VP, 


provided that the applied pressure p can be expressed in the form 
p= VP. 
Also, from (13a), 
2 
= a Viv = 0, 
Shiu —2 10uh 
A + =a Vid = = (w* + 26), 


and, from (12c), 


1l—-yAt+e 4 itp = Sa (1 + )Vie 


With the help of (15), (17b) and (17c), a simple calculation now gives 





We Ft 3(1 — 3nViP 

>= “3a, " {009 + 20@) + w(@) + sa} + Sin + 20h” 

wt = - +I {200 + 29@ + we) +) +E {20 + va} 
16u J 5yu 





=a, wa n)ViP 
4h*u 10hu , 
2h*n 


~ dam {or (2) + aa} - OP + (1 + VIP 


9] 


(12b) 


(12c) 


(13a) 


(13b) 


(14) 


(15) 


(16) 


(17a) 


(17b) 


(17 ¢) 


(18a) 


(18b) 


(18¢) 
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r = 2h T o+iv, (18d) 

az” 
vy, = we 2 sy (we + 26 + Div), (18e) 
DP m4 ~ + iy). (18f) 


The functions Q(z), w(z) are arbitrary functions of z. These equations for the weighted 
displacements and for the stress-resultants and stress-couples are in a compact and 
convenient form for applications. With them are used formulas for the transformation 
of stress-resultants and stress-couples, which are easily found by integrating Eqs. (5c) 
with respect to Z, after multiplication by a suitable function of Z. 

6. Results for certain aeolotropic materials. Only small changes are needed in the 
theory when the material of the plate is aeolotropic but is transversely isotropic in 
planes parallel to the faces of the plate. The stress-strain relations (3a) to (3d) are 


replaced by 





8 = (ce, + en) VF + F) + De oo (19a) 
= os VXF + F) +e (19b) 
GO, = Cig V1 C33 97” : 
2 nu 
® = &(e,, — Cy2) : A (19¢) 
~ o 2(e ) 
y= de 2 (2 9F + w (19d) 


Since the details of the analysis are similar to those for an isotropic material only the 
final results are given here. Thus 





C93 Pane 
¢ = as ( Q(z) + 22(2) + w(2) + sa 
(20a) 
4, 8esViP Beas 
10A(c1 133 sai Cis) 4h*(cy1C33 - Cis)’ 
™ : - ios h? ~ 
wt a — (206) + 22(2) + w(z) + sa} += {oe + wa} 
8(C31C33 — Cia) Dae 
(20b) 
_ ___ BeasP 8 {1 _ tu opp 
Qh? (ce, 1Caa — Cia) + 5h {2 C3;C33 ~~ Cis Vil ‘ 
3 fe ( — 
A= = clea + ea) = Beal + va} 
3 C1033 — Cis 
(20c ) 


a 2h*crs(Cr1 ex? Ci2)p i‘. Ca3(Cr1 + ¢i2) — 13 Vv? 2p 


5(Ci1C33 — Cis) €11C33 — Cis 











228 NOTES (Vol. VII, No. z 


16h°(e1, — C12) 8 
r= eu C12) ; ~@+ @, (20d) 
o ie 
Ohes, ¢ , 
Vv, = as = (w* + 2h + 2iy), (20e) 
. Oz 
D* = 4 9 (¢ + iy) (20f) 
> ae oa ; 
Oz 


The function y satisfies the equation 


y—- U2. = Wl = ©, (20g) 


SOME REMARKS ON THE FLAT PLATE BOUNDARY LAYER’ 
3y J. A. LEWIS anp G. F. CARRIER (Brown University) 


1. Introduction. In a recent paper’, the Blasius solution for the boundary layer flow 
past a flat plate was supplemented by an investigation of the character of the flow field 
near the leading edge of the plate. In that paper, the leading edge solution (of the Stokes 
type) was matched numerically to the Blasius solution. At first glance, therefore, it 
would seem desirable to find the solution of the Oseen type and either verify or improve 
this match. In the present paper we shall compute the exact solution of the Oseen 
equations associated with the flow past a semi-infinite plate and shall present arguments 
which verify our belief that this solution is physically not acceptable. However, we shall 
introduce a modification of the Oseen linearization such that the exact solution of the 
modified equations completes the flow pattern with a reasonable degree of accuracy. 

We shall also indicate an iteration procedure from which the exact flow pattern 
(i.e., the exact solution of the non-linear equations) can be obtained as the limit of a 
rapidly converging sequence of functions. Since the actual calculation associated with 
each step of this iteration would be very tedious, and since the solution given by Blasius’ 
and the present result give the most interesting information, we shall not complete the 
integrations. 

Finally, we shall indicate the ‘modified Oseen solution” for the flow past a flat 
plate of finite length. Again, the interest does not seem to justify the necessary algebra 
so that no numerical results are obtained. 

2. The Oseen linearization. The equations which govern the flow of a viscous, 
incompressible fluid 


ou; op Ou; i 
je p <P — (1) 
Ox; Ox; Ox; Oz; 
Ou; 9° 
“i = 0) (2) 
Ox; 


*Received Nov. 5, 1948. 
1G. F. Carrier and C. C. Lin, On the nature of the boundary layer near the leading edge of a flat plate, 
Q. Appl. Math. 6, 63-68 (1948). 
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can be linearized in a manner first suggested by Oseen’. For our problem, the arguments 
leading to such a linearization might be the following. The stream function and velocity 
functions defining the flow must have a branch point at the leading edge of the plate 
(i.e., the origin) and, furthermore, the velocities vanish at this point. This implies that 
the third term in Eq. (1) will provide much larger contributions near the origin than 
will the first term. Moreover, the velocity at large distances from the plate is essentially 
the constant free stream velocity, uo . Thus, one replaces Eq. (1) by 


Ou; , Op Ou; 7 

“ Ox, < a its ax, az, (3) 
and anticipates that the solution of Eq. (3) will demonstrate the proper qualitative 
behavior very near the plate and very far from the plate. However, if the very low 
speed field is separated from the “far field’? by a region in which the non-linear terms 
predominate, we shall be faced with a phenomenon not unlike the Stokes phenomenon 
encountered when investigating asymptotic solutions of ordinary differential equations. 
In fact, we cannot anticipate that the solution of Eqs. (2) and (3), which holds for the 
far field, will be the same solution which holds in the slow flow region even though in 
each region the asymptotic form of the flow functions must be solutions of these equa- 
tions. Nevertheless, it is profitable to find the solution of these Oseen equations and 
then utilize the extra information at our disposal to modify them and render the theory 


more applicable to our problem. Therefore, let us define u, = uo[l + ¥,], vw = —Uoyz ; 
L = Ul, /v, ¥ Unt. /v, Where X, , X2, U; , Up are the physical coordinates and the asso- 
ciated velocity components. With these definitions, Eqs. (2) and (3) become 
fe] 
a(a a“ 2\y = 0 (4) 
Ox 


where A is the Laplacian operator. The boundary conditions require that y. , ¥, — 0 
far from the plate and that y, = 0, ¥, = —1 on the plate, (i.e., on y = 0, x > 0). We 
note here that y, ¥, , ¥,,, are continuous across the plate but that y,, will be discontinuous. 
Indeed, let us define f(x) to be the discontinuity in y,, across the plate and note that 
f(x) vanishes when x < 0 and that f is twice the velocity gradient (du,/dx, in dimen- 
sionless form) along the plate. We may most readily formulate the problem by the use 


of Fourier transforms. In fact, we define 


VE, n) = | Y(x, ye '**™ dx dy, (5) 


multiply Eq. (4) by exp [—7(é + ny)], and integrate over the infinite region. Under 
our boundary conditions, integration by parts with respect to y leads to 


© 


a & > a3 ; 2 4) ~ — ) ° —itz > 
[ | (S-r(5-z2- 7) ve, yer de dy = | — infde de, 
d CO. / - . o—a 


where f(x) is defined in the foregoing. Defining f(£) in the obvious manner and completing 
the integration, we obtain 
+ nv + n° + py = — inf. (7) 


2H. Lamb, Hydrodynamics, Dover Publications, New York, 1945, p. 610. 
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At this stage, we can avoid future difficulties by observing that Eq. (7) is the limit 
(as k — 0) of the equation 


+ P+ WYP + E+ DE — ky = — inf. (8) 
We also note that the boundary conditions y, = —1 on the plate may be written 
as the limit, as a — 0, of ¥,(z, 0) = —e”°** on x > 0. Therefore, we now are prepared 


to solve Eq. (8) for ¥ and f according to the boundary conditions, with the obvious 
intent of carrying out the indicated limiting processes. Let us now define the function 


o(é, y) to be 


; a 
ag, y)=s- | VE ne dn 


oF Jz 


(9) 
= | v(x, yje~*** dz. 


The latter identity follows from the conventional inversion formula. Solving Eq. (8) 
for ¥ and performing the integration in Eq. (9) by the use of simple contour integrals, 


we obtain 


flexp {— |y | @ + k*)'*} — exp {—|yl|E+O"E-—m'""}] ag 


t 


ol, y) = (1 — KE — 2h) 





and, in particular, 


a ae , 
OK) = eH E + iby + EF OT - 


However, our boundary conditions are associated with g, , and it is convenient to define 





o(é, 0) = Uo(—) + UE) 
where 
lip(E) == | ” he e dx = ‘ = a (12) 
and 
0 ° 
u(~) = / u(xje*** dz; (13) 


that is, % is the transform of y, along the plate and u, , of y, along the line y = 0 
x < 0. The latter function is, of course, not known. Using these substitutions in Eq. (11) 





‘ f© 
q //( _— => 7 > -s r 7 5° 14 
u,(E) + L, g ta) QE as ik)'“*[(é + ik)'”” -}- (é + i)*”*] ( 4) 
We can now determine f and %, by the Wiener-Hopf method.* We note that, by 
definition, %,(£) is an analytic function of the complex variable in some upper half- 
plane which includes the real axis, provided we anticipate that y, should decay ex- 


*This technique has been extensively used, of late, on diffraction problems by H. Levine, J. 
Schwinger, A. E. Heins, and others. 
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ponentially in (—z) on the negative z-axis. Similarly, f(£) is analytic in some lower 
half-plane including the real axis. Therefore, when we rewrite Eq. (14) in the form 


mE + i)? + E+ 7) + e+ H+ E+ 97) 
— [a +b)? + (@ + 1)")} = HOE - iw"? (15) 


~ "lla + BY + @ + 1)", 


—& — ta 





the left-hand side is analytic in the upper half-plane mentioned before and the right 
side is analytic in the lower half-plane. Thus, by the arguments of analytic continuation, 
each side is equal to a constant, which is readily seen to be zero. Consequently, in the 
limit k — 0, a— 0 


a i\'? 

je = 22)”, (16) 
and we find by the conventional inversion formula that 

f(x) = 2x)”. (17) 


But f(x) is twice the velocity gradient along the plate. Thus we obtain .555/z'” for 
the velocity gradient, as compared with .332/z'” which, of course, is the Blasius result’. 
Thus, we must conclude that the region anticipated earlier in this section exists and 
that the solution is not valid. 

3. The modified Oseen solution. It is most convenient to introduce our proposed 
modification of the Oseen linearization with a minimum of a priori argument and, 
having found the desired solution, to establish justification arguments using the result 
to strengthen these arguments. We therefore merely suppose, for the moment, that the 
inertia term of Eq. (1) should be replaced, not by udu;/dx, , but by u*du;/dx, , where 
u* is some weighted average of the horizontal velocity component. It is evident that 
different methods of averaging would give different u* values so we may write u* = cUo 
where we expect that 0 < c < 1. Reserving the privilege of choosing c after the equa- 
tions are solved, we replace up by cup in Eq. (3), and Eq. (4) becomes 


a(a _ o\y ~—" (18) 
It is evident that the solution can be found as before, and in fact, we now obtain 
S(c, 2) = 2(¢/rx)'” (19) 
where f is again the jump in y,, across the plate. 
It can also be shown after some manipulation that 


v=s-| a6 we ae (20) 


ean be written 


ax 


v=] [Gr + 2)/2er*P7(1 — &-) de, (21) 
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and that, in particular 


u’ Y, = | — {[(r + x)/2em"*]'7*(1 — eh" "”)} dx (23) 
OY . 
where”? = 2° + y’. 

We now observe that the velocity gradient at the plate which is associated with 
the solution agrees with the Blasius result when c is chosen to be about .35. We also 
observe that the character of the far field flow is unaffected by changes in ¢ provided 
c > 0. That is, the leading terms in the dimensionless velocities u’ and v’ are essentially 
independent of ¢ for large distances from the plate. Thus the following observations 
are valid: the modified Oseen equation (18) is as valid as Eq. (4) for predicting the 
far field behavior and the Stokes (slow) flow field for any positive c. The choice c = .35 
allows us to use a single solution of Eq. (18) to predict these asymptotic flows and also 
to approximate the flow in the intermediate region. Since any Oseen attack is intended 
to do no more than this, it appears that this modification constitutes an improvement 
(for this particular problem) over the conventional attack. It must be noted, however, 
that we needed to know the value of y,, at the plate in order to choose our correction. 

There are various ways in which one might formulate a criterion for choosing c 
when the velocity gradient at the obstacle is not known. We might, for example, demand 
that the average (along the negative real axis) of the neglected inertia term be made 


to vanish, that is, that 


This, in fact, vields c = 1/2. A second way would be to weight this average. For example, 
if we wished to make the inertia term have mean value zero in a wedge shaped region 


ahead of the plate, we could require that 


Unfortunately, however, this integral does not exist, even in the Abelian sense. However, 
using the argument that when w’ is very small the value of c does not matter, we might 


weight the inertia term with w’ itself. That is, we might require that 


f ‘ ; ou’ 
| u(x, O)}1 + 2 c| — dz = 0. 
" OZ 
This produces the result c = 1/3 which is, of course,.an excellent approximation to 


the value formerly found. Since this averaging makes sense physically and yields the 
proper result in at least one case, we suggest that, in general, the Oseen linearization 
associated with the viscous flow past an obstacle can be improved by using u* instead 
of uo , where u* is found by the foregoing averaging process. However, in view of the 
fact that this yields u* CUo up/3 Without regard to the form of wu’ (and therefore 
without dependence on the symmetric obstacle shape), one must immediately conclude 
that u* ~ u,/3 for any symmetric obstacle. The flow so obtained can then be expected 
ear the obstacle, far from the obstacle, and 


to be valid in the Stokes slow flow region 1 
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in some strip along « <0, yy < yo. This is certainly as large a region of validity as 
one could expect from such an asymptotic solution. 
4. The exact solution. If we now write down the exact equation for the stream fune- 


tion Ww (as defined before) we obtain 


a(a —C = \y = F(y), (24) 
OX 


where /(wW) is that portion of the inertia term not contained in ey, . We can now find 
a Green’s function for this equation (treating F as a non-homogeneous term) and write 


aa 


y Yo(xr, y) Ft I G(x, %o 5 Ys YE (Plo, Yo)) dxo dYyo . (25) 


where Wo is the modified Oseen solution just obtained. We note that if a function y*; 
which was very nearly the exact solution of Eq. (24), were placed in the integral term, 
he function w obtained from the integration and addition of y. would, in general, be 
still more accurate. This, of course, is a standard procedure. However, we have at our 
lisposal the fact that for | are tan (y/| 2 |) | < (say) 7/18, the Blasius solution is very 
weurate, and for the remaining field the present solution seems adequate. Therefore, 
f one divided the integral of Eq. (25) into two parts, chose ¥* as suggested here, and 

rried out the integration, he could expect a marked improvement in the result. Actually, 
the results already in hand seem good enough so that the amount of labor required is 
essentially prohibitive. Thus we shall not pursue the question further except to write 
lown the appropriate Green’s function. When y* is odd in y, the integral of Eq. (25) 

| always obey the conditions implying y = ¥, = 0 on y = 0, x > O, if the Green’s 


+) 
8) Is 


kK, (x Xo) + (Y%1 — Yo) |? da, dy 


where 0 <arg2'* < aie < arg (2)'” < 2x; and K, is the modified Hankel 


Inetiohn of orce ZCTO 
5. The plate of finite length. The technique by which the foregoing problem was 
olved cannot be directly utilized when the plate is of finite length. However, if we let 
ate extend over the interval 0 < « < A, we may write (for the jump in ¥,,) 


2. a,x" for O<r<aA 


0 for 2 < 8, and oA 
When this is done, ¥,(a, 0) in 0 < «2 < A can be written 


W(x, 0) > alx” exp (kx). 
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It ean also readilv be show 


aft 
a 


The series will fail to 
1.1] 


consideration does not 


con 
to exceed 


verge, 


NOTES 


n that one may successively find the 


ii) | on Q < A 


MH 


of course, at A. However, 


the value obtained for the semi-infinite plate 


the unmodified Oseen procedure 


AC 
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a’ (and hence a,) so that 


n 


a) can readily be shown 


and thus the finite plate 


more valid. In fact, one 


might use the modified procedure to improve the finite plate solution although, for this 


problem, we do not know 
In any event, 


uny detail 


1 
Ts) 
i 


velocity gradient at the plate as we did for the long plate. 
the numerical evaluation of the velocity field will not be pursued here in 


AN ITERATIVE METHOD FOR FINDING THE SMALLEST 


sy JOHN 


EIGENVALUE OF A MATRIX 
I 


> 
> 


LOHMAN 


Brown Unive 


The eigenvalues and eigenvectors of a matrix A are often found by repeated ‘pi 
multiplication of an arbitrary column vector v by the matrix A. The products A"v con 
verge to the eigenvect« corresponding to the eigenvalue \, largest in absolute magni- 
tude: convergence is | \, is much larger than any other eigenvalue of A. 

If one already h ) imate values of the eigenvalue of A, this iterative method 
can be modif O give ore rapid convergence. Aitken’ suggests subtracting calar 
matrix al (I denoté T t matrix) from A; a is chosen so as to give the maximum 
separation of Ii “ carries out the iteration using a polynomial P(A) in A, in 
place of A itse lf; riz » chosen that the igenvalues hot immediately sougnat are 
nearly zero as possibl 

Iteratio ( ( mallest eigenvalue of natrix A commonly perform i using 
the inve em | | note describes moaned pro edu i( viding to mot rap | 
converge! 

Let the or i A be..22. .4, 1, Z id assume we know them only to within 
LO% Simple ersion Oo ® ld | with eigenvalues Lo 2 5: the large { ivwen- 
value of A correspond yr to the smallest eigenvalue i 1.8 times as large as thr 
second-larges \ssum¢ elieve the smallest eigei f A to be .2 instead of .22; 
subtract t scalal 21 from A, and then invert ting matrix (A 2] 
has eigen 50, 5, 1.25, .56. But we ean still improve tae ‘paration by subtracting 
inothe } \\ noose the mean yf the se yi v und the small wen 
Values 1, a i at \ yuld be yppropri ( som C1? l Liues Were n FALLIVE for 
maximum separation of the largest and the ond larg renvalu Our final modifie 
matrix, (A 2] 2.781, has eigenvalues 47.22, 2.22, 1.53, —2.22. The ratio of 
largest to si 1 larg o’ value ha ) l {from 1.8 to 21.3. 

Since nodified has all but one of its eig wlues tightly clustered about 
zero, 10 1 ) suitab form from which to determine iese other eigenvalues 

I 2 ) I rk was supported under Navy con 

A.C. 2 P I s Mdinburgh 57, 269-270 (1937 

WM. I i, fath. 5, 329-345 (1947 
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AN ALTERNATE PROOF OF THE CONSTANCY OF CIRCULATION* 
By WALLACE 1) H AYES (Brown { niversily 


In a non-viscous, barotropie fluid with only conservative body forces the circulation 
about any circuit that moves with the fluid is constant. The direct proof of this theorem 
due to Kelvin is the one customarily given but other proofs due to Cauchy and Helmholtz 
are also available. By the application of Stokes’ theorem the statement may be changed 
to the constancy of the flux of vorticity through a surface bounded by a contour, both 
of which move with the fluid. The theorem in this latter statement is here proven 
directly, using vector language. 

Let A be a vector which satisfies the equation 


) C 
OA 8 OB +g. VA = (A-V)q — A(Y-9), (1) 


q is the velocity vector and D/ Dt the total time derivative. 
\ differential surface element dS satisfies the kinematic condition 


Dds 
Di (dS-V)q + dS(V-q) — dS X (V X q). (2) 
This condition mav be derived by considering dS to be the vector pro luct (a *& b) of 


lifferential distance elements Iving in the surface. The total derivative may be 


Dia Xb 
a qx<b+a (b q 
Dt 
a X b) ‘ q 
Ls tO on (2). Taking the sealar product of Eq. (1) with dS and of Eq. (2) 
A and ad vield 
D(A-d$S (3 
Dt \e 
on ght cancelling because of vector identities. If the integral of this 
be taken over a surface which moves with the fluid and its bounding contour, 
D fy 
A-dS 0 | 
Dt I 


tor through the surface is constant. It may be noted that V-A 
he same continuity equation as does the density. 


B satisfies the condition 
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The result (4) together 


where y is an arbitrarv scalar, its curl satisfies condition (1 


with Stokes’ theorem giv 


2 ) B-dl = 0, 6) 
i.e., that the line integral of B about a circuit which moves with the fluid is constant 
This may also be proven directly. 

or q satisfies (5) as a ly namical condition and the vorticity vector 
taney of flux of vorticit 


7 


The velocity vect 


V X qsatisfies (1). Thus (4 


and (6). the constanev of circulation 


< 


gives the theorem for the cons 


A GENERAL APPROXIMATION METHOD IN THE THEORY OF PLATES OF 
SMALL DEFLECTION 
REMARK ON MY PAPER 
YTARTERLY OF APPLIED MATHEMATICS, 6, 31-52 (1948 


By M. Z. KRZYWOBLOCKI (University of 1 


V. Southwell and Dr. A. Weinstein 


The author wishes to apologize to Professor R 
for not mentioning their methods in his paper. As Professor Southwell wrote in one of 
his papers, it is unlikely, in present circumstances, that any review of numerous m ds 
which appeared lately in the theory of plates could be made complete and any judgement 
regarding the merits of different methods is premature 

‘mentioned in the Introduction >| 


So let u Lad just two more methods to thos 
and on the bottom of p. 42 of the author’s paper. Weinstein’s variational method allows 
one to solve probl governed by the differential eq tion 
AAw = gq, i.e. with ransve \. Weinstein, Mémorial des sciences 
mathématiques, No. 88, 1937; N. Aronszajn and A. Weinstein, Am. Jour. Math. 64, 
623-43 (1942); A. Weinstein and D. H. Rock, Quart. App. Math., 2, 262-6 (1944); A 


Weinstein and J. A. Jenkins, Tr 


ature 1s applicable to all cases as may 


arbitrary shapes 


, } 
rse load only see, 


R. S. Canada, See. IIL, 60-67 (1946). Southwell’s 


be seen trom 


relaxation method, from l 

examples in Southwell’s “Relaxation Methods Applied to Engineering Problem nd 
from the following papers: L. Fox and R. V. Southwell, Trans. Roy. Soc. (A) 239, 419-460 
1945): D. G. Christopherson, L. Fox, J. R. Green, F. 8. Shaw, and R. V. Southwell, 


tid. 239. 461-487 (1945 
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Tables of Be SSé | fun haon of fractional orde Prepared DY the ( omputation Labor wory 


of the National Applied Mathematics Laboratories, National Bureau of Standa 
I New York, 1948. XLII + 413 pp. $7.50 


Volume I. Columbia 


University Press, 


The first part (pp. 1-272 yntains tables for J,(a rv = +1/4, +1/3, +2/3, +3/4 an ingin 
from 0 to 25, at intervals of 0.001 for small values of 2 ip to 0.9 for 3/4, 2/3, to 0.8 fo 
= 2 1/4. to 0.6 for 1/4, 1/3, and to 0.5 for Zia, a/4 I intervals of 0.01 Lr ere 


y = —1/3 
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valu The second part (pp. 273-383) gives the functions A,(r) and B,(x) defined by J,(7) = A,(z) 
cos rv/2 r/4 B(x) sin (x — wv/2 r/4) for the same values of v and for x = [25 (.1) 50 (1) 
500 (10) 5,000 (100) 10,000 (200) 30,000]. All functions are tabulated to ten decimal places 

ibles facilitate interpolation and give the zeros of J,(2). 


A number of aux ul 
W. PRAGER 


Table of the Bessel funet ONS of the first kind of orders zero and one By the Staff of the 
Computation Laboratory (The Annals of the Computation Laboratory of Harvard 
Volume III). Harvard University Press, Cambridge, 1947. $10.00. 


University 
ission of the Bessel functions, the computations, the interpolations, and a list 


Following a brief 
| Jessel functions Jo(r) and J;(r) are tabulated to eighteen digits. The argu- 


the Bessel coefficients, the Be 
001 for 0 < x < 25, and in steps of .01 for 25 100. 
G. F. CARRIER 


lis 


rt < 


Bv the Staff of the 


; 


ions of the first kind of orders two and three. 


Tabi f the Bessel funct 
itory (The Annals of the Computation Laboratory of Harvard 


( omputation Labor 
- TV). Harvard University Press, Cambridge, 1947. $10.00. 


University Volume 
and J3(x) are tabulated to eighteen digits. The argument varies in steps 


J o(x are 
25, and in steps of .O1 in the range 25 < «© < 100. 


G. F. CARRIER 


f ON he rang 
Tables of the Bessel functions of the first kind of orders four, five and six. By the Staff of 
he Computation Laboratory (The Annals of the Computation Laboratory of 
Harvard University —Volume V). Harvard University Press, Cambridge, 1947. 
$10.00 
Che Bessel fur ‘ J <x) and Je(x) are tabulated to ten digits. The argument varies in steps 
H } rang ) < < 25, and in steps of .OL in the range 25 < xr < 100. 
G. F. CARRIER 


Lables of the Besse lf inctions of the Sirst kind of ord rs seven, cight and nine. By the Staff of 
the Computation Laboratory (The Annals of the Computation Laboratory of 
Harvard University -Volume VI). Harvard University Press, Cambridge, 1947. 
$10.00 

are tabulated to ten digits. The argument varies in steps 
25 <r < 100. 


» 


Ji(xr), Js(x) and J9(x 


5, and in steps of .01 in the range 
G. F. Carrter 


) 


Tab cs of the B s! f 
of the Computation Laboratory (The Annals of the Computation Laboratory of 


essel functions of the first kind of orders ten, eleven and twelve. By the Staff 
Volume VII). Harvard University Press, Cambridge, 1947. 


Li 


Harvard University 


$10.00 
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The Bessel functions / Jii(z) and J,0(x) are tabulated to ten digits. The argument varies in 


steps of .0O1 in the range 0 < 25, and in steps of .O1 in the range 25 < x < 100. 


G. F. CARRIER 


Table rs of the Bi gsel fun l20NS of thie Jt st / and of orde rs thirtee hn, fourtee MH and fifteen 13\ thie 
Staff of the Computation Laboratory (The Annals of the Computation Laboratory 
of Harvard Universitv— Volume VITI). Harvard University Press, Cambridge, 1947 


$10.00. 


The Besse unctions J ] r) and J il tabulated to ten digits The argument varies in 
steps of OO1 1 the range 0 25, and in ste] of .O1 in the range 25 < < 100. 


Table r of the Bessel functio of the first kind of orders sirleen thre ugh fwenty-seve) Iv thie 
Staff of th Computation Laboratory Che Annals of the Computation Labora Or 
of Harvard Universitv—Volume IX). Harvard University Press, Cambridge, 1948 


$10.00 


The Bessel fu sons -cl ] J uulated to ten digits. The argum 
steps of .O1 r OQ OO 
cz, FF, 4 
Tables of the Bessel f s of the first ki) f ord twenty-cight through thirty- 13) 
4} ; 4 ; ms t+] ri 
the Staft i the lL omp iTIO! Laboratory i fie Annals of the ( omptutation I. MOraloO 


of Har 1 U1 ¢ { Volume X Ha rd Universitv Press. Cambride 18 


$10.00 


The | / lig The argu 
steps of .O 0 ) 
( I. ¢ 
Eleven a pia fh / / aul significa / 
By Enzo Caml York, Dover Publications, 1948 154 pp. $3.95 
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The coefficients A,, of this Neumann series are given by the formulae 
Land 


{ A. (k,n) = J{k) 


{ A.(k, n) = J,..(k) + (—1)"J,..(k) (m > 0) 


which follow at once from the well known addition theorem for J;,(x). 
These auxiliary tables were used by the author to build up the values in Table I for n even and > 2. 
and J; were taken from Meissel’s well known table). The values for n odd were then computed from 


/ 
the recurrence formula 
9\-1 / 
Jora(2) = x(4r + 2) [Jo,(2) + Jor+2(2) 
An unusual feature of Tables I and IT is a “cheek”? column showing for each value of x the errors 
le in substituting the tabulated values into the identities 
| y Js kei. 
sm 2 2 (—1)"J..4;(2) 
r/2 S* (Qn + 1)J (xr) 
— . 
Phe ire a wctorily small. Unfortunately these checks were not applied to the proof sheets. 
Otherwise the following tvpographical error (noted by the author) might have been detected in time for 
tion: p. 47. J;(6.983), for .35385 88467 2 read .35385 88367 2. 
\nother ident ised by the reviewer in spot checking the table, is the Kapteyn series 
| 2) i+ 2 > J (nz). 
Phere ig no d n of interpolation nor are any differences tabulated. . 
| compact a relatively inexpensive table of J/ will prove very useful to almost all workers 
u inetions. As explained in the introduction, the table emphasizes the Bessel 
the highe rders n. Its choice of the variables n, 2 will make it especially useful in dealing 
umanr 
S* a,J,(x) 
has 
H r for me a tions, such as certain Fourier-Bessel series, the limitation 2 < 10.5 may prove 
ictive. hi he table gives very few waves of the wavy surface P Thus, within the 
ibulati ind J; have onlv 3 zeros, Jo and J; have 2 and J,, J _ Js have one each. J; and 
| maxil hile Je and J,» only one point of inflexion each. All other funetions are mono- 
| 

I} ffset p typescript is verv clear and well set out 

D. H. Leumer 
Introduction to the differential equations of physics. By L. Hopf. Translated by Walter 

Nef. Dover Publications, Inc., New York, 1948. v + 154 pp. $1.95. 

Chis book is a troduction to the differential equations of mathematical physics. The exposition is 
le nd coneist ne might like the treatment to go further than it does. The text should be very 
1 i] to seniors and beginning graduate students in physies and chemistry 

Che contents he book are divided into seven chapters. The first four chapters are devoted to 

discussion of the o1 ry differential equations of the mechanics of particles and the simpler partial 


diff tial equatio physies as applied to potential theory, heat conduction, potential flow, and elec- 
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tromagnetic theory. Chapter V is devoted to solutions of the partial differential equations by eigenfunc- 
tions with a discussion of the product method, Fourier series methods, expansions, regular and asympto- 
tic, the Fourier integral. Chapter VI discusses the solution of partial differential equations by change of 
variables methods; this section includes a brief discussion of conformal transformations. Chapter VIT is a 
discussion of the use of singularities in obtaining solutions of equations. The chapter contains some dis- 


eussion of the use of Green’s Theorem, and the transition from differential equations to integral equations 


R. TrRvuewi 


Introduction to applied mathematics. By Francis D. Murnaghan. John Wiley & Sons, Inc., 
New York and Chapman & Hall, Ltd., London, 1948. ix + 389 pp. $5.00. 


The text b gin with a devel pment of the algebra and calculus of vector analysis based on a (Car- 
tesian reference system. The material is extended to include vectors in an n-dimensional space and then 
Fourier series are introduced from the function space point of view. In chapter four the vector calculus 
is carried over into the curvilinear coordinate system. A chapter then follows on the conventional method 
for solving the more elementary problems involving the Laplace equation. Following a chapter on 


Legendre and | functions, concerned largely with power series developments, zeros and orthog- 
onality, is a discussion of eigenfunctions and eigenfunction expansions associated with boundary valu 
problems. The final chapters deal with the Fredholm equation, the calculus of variations, and the opera- 
tional calculus. 

On the whole, the material is sufficiently elementary and presented in such detail as to be of int 
as a text for a senior « i1troductory graduate course in applied mathematics. 


G. F. CARRIER 


Introduction to omple triables and applications. By Ruel V. Churchill. MeGraw-Hill 
Book Company, Inc., New York, Toronto and London, 1948. vi + 216 pp. $3.50 
The algebra of complex numbers and the introduction of analytic functions are followed by a discus- 

sion of the elementary functions from both the analytic and geometric point of view. The integration of 

complex functions is developed and power series are discussed with a very brief section on Laurent series. 

The residue theorem is discussed with emphasis (via examples) on the evaluation of definite integrals. 

Conformal mapping is discussed and applied to the conventional heat transfer, electrostatic, and fluid 


mechanics problems. The Schwarz-Christoffel transformation is treated in detail and a brief section on 


analytic continuation and Riemann surfaces complete the text. The appendix contains a short ‘‘dic- 


tionary”’ of simple conformal mappings. The book should be of interest to those concerned with the use 
of function theory in engineering problems and similar applications. 


G. F. CARRIER 


Principles of jet propulsion and gas turbines. By M. J. Zucrow. John Wiley & Sons, Inc 
New York and Chapman & Hall, Ltd., London. xiv + 563 pp. $6.50. 


A discussion of the properties of gases and the fundamental laws which they obey is accompanied by 
detailed material on airplane performance, propellers, gas turbines, compressors and jet engines. High 
temperature alloys are also discussed. A large proportion of the text is devoted to related topics rather 
than to the principles of jet propulsion. The emphasis is on the presentation of engineering information 
and the text is perhaps most appropriate for an engineering course at the senior or early graduate level. 


G. F. Carrier 

















SUGGESTIONS CONCERNING THE PREPARATION OF MANUSCRIPTS 
FOR THE QUARTERLY OF APPLIED MATHEMATICS 


Mathematical Work: Only very simple symbols and formulas should be typewritten. All others should be 
carefully written by hand in ink. Ample space for marking should be allowed above and below all equations. 


Greek letters used in formulas should be designated by name in the margin. 

The difference between capital and lower-case letters should be clearly shown; care should be taken to 
ivoid confusion between zero (0) and the letter O, between the numeral one (1), the letter / and the prime (’), 
between alpha and a, kappa and k, mu and uw, nu and », eta and n. 

All subscripts and exponents should be clearly marked, and dots, bars, tildes, etc. over letters should be 


avoided. 
Square roots should be written with the exponent $ rather than with the sign /. _ 
Complicated exponents and subscripts should be avoided. Any complicated expression that reoccurs fre- 


quently should be represented by a special syinbol. 
or exponentials with lengthy or complicated exponents the symbol exp should be used, particularly if such 


exponentials appear in the body of the text. Thus, 
exp [(a’ + b’)'”’] 1s preferable to e‘ 
Fractions in the body of the text and fractions occurring in the numerators or denominators of fractions 
should be written with the solidus. Thus, 


a*+b*)2/2 


cos (r2/2b) . —_ 3 
cos (wa/2) is preferable to ——— - 


cos = 
2b 


In many instances the use of negative exponents permits saving of space. Thus, 
ait ce ° sin u 
J u~ sin u du is preferable to | a du. 


Whereas the intended grouping of symbols in handwritten formulas can be made clear by slight variations 
his procedure is not acceptable in printed formulas. To avoid misunderstanding, the order of symbols 


in spacing, th 
should therefore be carefully considered. Thus, 
(a + bz) cos ¢ is preferable to cos ¢ (a + bz). 
In handwritten formulas the size of parentheses, brackets and braces can vary more widely than in print. 
Particular attention should therefore be paid to the proper use of parentheses, brackets and braces. Thus, 
{{a + (b + cz)"] cos ky}? is preferable to ((a + (b + czx)") cos ky)’. 
For further suggestions concerning cuts, bibliography, footnotes, and abbreviations please see the inside back 


cover of the January issue of 1949. 





Al 4 K Wy ] l on the Dover mailing list? Do you receive regular announcements of such 


books as Fr. A. Willers’ PRACTICAL ANALYsIs (A 422-page presentation of 
the numerical, graphical and many of the instrumental methods of obtaining a specific numerical answer in 


low priced and important scientific 


applied mathematics. Includes a special section on modern calculating machines by Tracy W. Simpson. Now 
ATIONS OF NUCLEAR Ppuysics (Compilation of thirteen classic papers on nuclear 


$6.00): Robert T. Bever’s FOUND 
Cockroft, Fermi, Curie, et al. Large 122- 


physies, unabridged and in the original language of their authors 
page bibliography of approximately 5000 papers on nuclear physics. Only $2. 95) ; H. 8. Carslaw’s INTRODUCTION 
TO THE THEORY OF FOURIER’S SERIES AND INTEGRALS (Formerly $7.00, now $3.95); Felix Klein’s ELEMENTARY 
MATHEMATICS FROM AN ADVANCED STANDPOINT (Volume I, in English pate covers arithmetic, algebra 
also translated, deals with geometry and is priced at only $2.95); 


and analysis, and costs $3.75. Volume II, 
umbi’s TABLES OF BESSEL FUNCTIONS (Eleven and fifteen-place tables of the first kind to all significant 


orders. Only $3.95); T. M. MacRobert’s sPHERICAL HARMONICS (second revised edition. Now $4.50); A. I. 
Khinchin’s MATHEMATICAL FOUNDATIONS OF STATISTICAL MECHANICS (First. english translation by G. Gamow. 
Only $2.95); H. M. Barlow’s MICRO-WAVES AND WAVE GUIDES (Only $1.95); Howard W. Emmons’ GAs DYNAMICS 
TABLES FOR AIR (Semi-stiff binding. Only $1.75). 

We'll be glad to add your name to our list, if you’ll write and tell us which field or fields you are most interested 


in—pure mathematics, mathematical physics, etc. Incidentally, you may order any of the above books direct 


from us. They’ll be shipped to you postpaid, fully covered by our ten-day, money-back guarantee. 


DOVER PUBLICATIONS, Inc., 1780 Broadway, New York 19, N. Y. Dept. QAM. 


Enzo Cz 
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method 
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NOTEs: 
J. W. Miles: A Note on a solution to Possio’s integral equation for an 
oscillating airfoil in subsonic flow 
M. Herzberger: The normal equations of the method of least squares and 
their solution 
A. E. Green: On Reissner’s theory of bending of elastic plates . 
J. A. Lewis and G. F. Carrier: Some remarks on the flat plate boundary 
layer 
J. B. Lohman: An iterative method for finding the smallest eigenvalue of 
a matrix 
W. D. Hayes: An alternate proof of the constancy of circulation 
M. Z. Krzywoblocki: A general approximation method in the theory of 
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HEAT CONDUCTION. With Engineering and Geological Applications 
INGERSOLL, University of Wisconsin; Orro J. Zope, Bell Telephone 
and ALFRED C. INGERSOLL, University of Wisconsin. 290 pages, $4.00. 


3y LEONARD R 
Laboratories, Inc.; 


An elementary text in Fourier Series and the mathematical theory of heat conduction, with 
many interesting and practical applications. 

By AprRIAN ALBERT, The University of Chicago. 164 pages, $3.00 
Presents an exposition of the analytic geometry of three-dimensional space. The aim of the 


book is to provide a modern and simpler treatment of the subject matter which permits easy 
generalization and fits the subject into its proper place in modern mathematics. 


THE REAL PROJECTIVE PLANE 
By H.S. M. Coxeter, University of Toronto. 200 pages, $3.00 


An introductory college textbook in projective geometry, including a thorough treatment of 
conics and a rigorous presentation of the synthetic approach to coordinates. 


INTRODUCTION TO STATISTICAL MECHANICS 
3y Ronatp W. Gurney, The Johns Hopkins University. 268 pages, $5.00 


Against a background of atomic physics, this book studies the behavior of large groups of 
particles, showing how the behavior of matter in bulk follows directly from the properties of 


the individual atoms. 
Send for copies on approval 


McGRAW-HILL BOOK COMPANY, Inc. 
330 West 42nd Street New York 18, N. Y 








